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7 _ FOUNDATION STABILITY FOR A SUBMARINE LIQUID SULPHUR PIPELINE 


By 
By Basil W. ‘Wilson, F. ASCE 


‘SYNOPSIS 


This p paper examines the vertical and tateral stability of a7- mile long hot — 
- sulphur pipeline, entrenched in the sediments between Grand Isle, Louisiana, 
and the offshore sulphur mine inthe Gulf of Mexico. Hot water under pressure 
a ina jacket pipe maintains fluid sulphur flow in an inner line at an average a 
temperature of 300°F. These hot lines are insulated within an outer comes 
‘Pipe. Under normal operation, the expansive tendencies of the hot lines 
blocked by end anchorages, prescribed for this purpose; the lesser tendency 
to expansion of the casing pipe is similarly prevented. The resulting com- 
“pression i in the composite pipeline « could involve buckling and | critical exposure 7 
of the pipeline in the absence of adequate soil support. The adequacy of this 
support and the nature of the thermal stresses in the pipeline are subjects of 7 
enquiry, having regard to the special nature of the foundation. The latter in- 
volves fine sands and silts over 50% of the pipeline length and ve very soft co- 7 ; 
7 _hesive sediment over the remainder. It is shown that the pipeline will be 
‘stable — buckling, but may suffer limited subsidence from thermal 
2 “pipeline will probably sink below trench level ina pocket of very weak sedi- — 
ment at two-thirds the distance from the shore. It should stabilize, = 
without suffering unduly severe strain from the sag. A finding of the study is © 
: ‘that initial prestressing of the casing pipe during the laying operation will be 
felled, by the resistance forces” of the sediments; from this 1 it is 


_ Note.—Discussion open until January 1, 1962. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE. This paper is part a 

of the copyrighted Journal of the Soil Mechanics and Foundations Division, Proceedings — 
he: of the American Society of Civil Engineers, Vol. 87, No. SM 4, August, 1961. 
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1cluded_ that soil resistance also renders unnecessary the use of costly 
7 anchorages for ni needless control of pipeline expansions. 


NATURE OF THE PROBLEM 

- The decision in 1956 to mine the large sulphur deposits discovered in 
_ 1949 offshore from Grand Isle, La., in the Gulf of Mexico, , led to the design of 

an over- water Frasch extraction plant, : unique in being located seven miles 
from shore in water of 44 ft depth. 2 Recovery here of liquid sulphur at = 

temperature of some problems Of in ‘a sometimes 


con 


K= COEFFICIENT OF THERMAL CONDUCTIVITY 


@= TEMPERATURE (°F) WALL, 14"0.0. API 

- \\, CASING PIPE (K,=26; 


0.81 "AIR SPACE 


FIBERGLASS INSULATION 
0.30 "WALL 0.D. 
SHOT WATER 
=2630 to 420; 6,2 300°F) 
wor WATER | WALL, 6'0.0. 
300°F 


=| SULPHUR 
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FIG. 1.—SCHEMATIC REPRESENTATION OF CROSS-SECTION OF PIPELINE | 

SIZES, MATERIALS, AND HEAT CONDUCTIVITIES 

= ile environment that, i it was decided, - could best be overcome with an 
underwater hot-sulphur | pipeline from the offshore s structures to the Grand 

Isle base.2,3,4 But design? and construction of this pipeline presented prob- 
lems concerning the foundation for which there were no immediate answers. 


a 


“The Grand Isle Mine; Freeport Sulphur Company’s Offshore Venture,” by C. O. _ 
- Z. W. Bartlett, and R. Feierabend, Soc. of Mineral Engrs., AIME, Reprint No. 
59H97, November, 1958. 


_ 3 “The Grand Isle Sulphur Mine Heating Plant and Underwater ’ Sulphur Line Design,” 
by C. M. Cockrell, ASME, paper no. 59-A-183, August, 1959. 
: 4 “Construction Features of a Seven Mile Hot Sulphur Line under the Gulf of Mexico, ” 0,” 
” E. J. ee Proceedings, ASCE, Vol. 86, April, 1 1960, p. 47. 
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‘FIG. : 2. .— LOCATION OF SULPHUR PIPE R PIPE LINE OFF GRAND ISLE, ‘LOUISIANA, 


GULF OF MEXICO 


(a) LOCATION MAP| 
4 
— 
| 
— 


1961 


TABLE 1. —DIMENSIONS AND PROPERTIES 


‘Moment of 


area, in | Neutral Axis, 


Air 


625 


with Rollers 


Hot ter line 
‘Annulus of Hot 

Sulphur Line 
Molten Sulphur 

Return Water Line 


Outer — 
(Return Water) 


Return Water 
Return Water 


The foundation was known to ci comprise extremely st soft sediments, deposited 
on the continental shelf in geologic time by the Mississippi River deltaic 
_ System. 5,6 The questions at issue concerned the anticipated behavior of the 
pipeline in such a medium, having regard to the high temperatures the steel 
pipes would acquire from passage of the hot fluids and the possibility that the 
; pipeline might be subject to exposure from the devastations of a Gulf hurri- 
cane, Reasonable assurance was needed that the pipeline would be both longi-_ ; 
1 -tudinally and laterally stable when entrenched in the sediments and that coastal “a 
_ erosion in the area of Grand Isle would not lay bare the line within the prob- 


life expectancy of the mine. The depth of burial of the and 


Sediments of the Continental Shelf ” by H. N. 


_ Proceedings, 8th Texas Conf. on Soil Mechanics and Foundation Engrg., Univ. of of Texas, 


Univ. of Texas, Tex, 1957. 


— 
4 
hoe, 
= 
| 
— 
f 
— 
4 | 7.625 0.300 6,90 46.20 &g 
= 
a 
" 
— 
— 


Material Coefficient | Film | Coefficient 


—| Thermal | Coefficient | of Linear 
Yield | ConductyK,} Expansion 
‘Strength, | in BTU per perft| ,in | 
“Fperbr 
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in | ft °F per! hr 


2.8 


ca 
2630 to 420 
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w arm Water 


271 Pipeline with fluids 


- additional ‘specific unknowns. It is i purpose of ‘this pz paper to describe the 
_ exploratory studies that yielded information on on the vertical and ‘thermal — 
stability of the pipeline; limitations of length prevent discussion here of the 
Notation.—The letter symbols adopted for use in this paper are defined 


Principal Features of the Pipeline.—A schematic representation of a cross- 

section of the embedded pipeline, as it was conceived in 1958, is shown in 7 E 
Fig. 1. ‘Fluid ‘sulphur was to be conveyed in the innermost line of 6 in. (out- © ; 

side) aioe. 2 arranged annularly within an insulated hot-water line of 7 5/8 

Heat from the surrounding water would maintain the sulphur in : 
the aes onan at about 300°F, and the insulation and annular airspace © 
round this assemblage were expected to prevent serious heat loss to the outer 
casing pipe (of 14 in. diameter) that was to function as external support and q 
| protection for the system. Not shown in Fig. 1 is a return-flow water-line 7 . 
that, in 1958, , was was envisaged as a 4 1/2 in. outside diameter pipe strapped to 


per ft | 
75-109 (Timcoat) o2 |] ... | 699 
(range) | ~ fe = 
4 
ty 


| the p 

line, 2,3,4 that however, does not feature in the coaaidasetines of this paper. 
~The location of the offshore mine (Fig. 2) made the effective length of the | 
pipeline approximately 35,000 ft. Over this distance the anticipated ate 
water and sulphur temperatures at the mine were 320°F, and their outlet 
temperatures 280°F at the shore. In enter to secure an equitable distribution 


pipe, _ when meomtes on site at a probable temperature of 80°F, should be pre- 
tensioned by end-pulls of about 200,000 Ib and anchored at both mine and 
_ shore-ends. The hotwater and sulphur lines, supported in the casing pipe at 
intervals by peripheral rollers, were then to be secured to the casing pipe 
& the mine-end and expanded at atemperature of about 190°F before fastening» 


their ends to the shore anchorage. The whole pipeline was to be laid in an 


TABLE. 2,—AVERAGE PROPERTIES OF PRINCIPAL. SEDIMENT TYPES ES 


Classification om _Con- ‘Atterberg Limits 
and Grain tent 
in % 


Greyortan, | 
loose, fine, a 
w/small shells 


and silt inclusions 
mm, | q 
90%<0.2 m.m, 
Soft, grey, mainly 
Montmorillonite, 

flliteand 
Kaolinite | 


Pertinent information concerning the given, or assumed 
for purposes of calculation, is contained in Table 
test borings the pipeline site showed that to a distance from the coast 
of approximately 11,000 ft the foundation consists of loose grey or tan fine 
co with frequent shell fragments and silt lenses. ‘Beyond this distance, the 
i medium comprises very soft grey clay with occasional sand lenses . The av- 
2 properties of these sediments are recorded in —_ 2; additional all 


results are presented later, 
Relative” distributions of the sediment types are en in cross-sectional 
- profiles along the pipeline site in Figs. 3, 4, and 5. These portray more spe- > 
- cifically the variations within the clay medium of wet density, ultimate un- 
confined compressive strength and ultimate shear strength. The full- a 


contours defining wee of these properties reveal a zone of very ' ry weak sedi- 


casing pipe in a position Immediately above OF below it. Final construction 
excavated trench 1tor / below the mualine and allowed a period for natural 
Type 
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OURED PROFILES OF WET DENSITY OF THE SEDIMENTS ALONG THE PIPELINE SITE 
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deeper made in the area, "rather westward of the line of 
borings, indicate that locally stiffer material may be encountered there in the 


aan loss in strength of the clay when disturbed is evident in 
Figs. 3 and 4 in which the profiles (b) present test results on remolded sam- 
ples. It must be supposed that the backfill material that ultimately confines — 
+) the pipeline in its trench will suffer loss of strength comparable to that shown, 
die In the analysis of the foundation problem that follows, attention will be | 
directed first to the vertical stability of the pipeline and then to the stability — 
of the line under thermal effects. 


VERTICAL STABILITY OF THE PIPELINE 


oa ‘Buoyancy y of the Pipeline.—For the dimensions of the pipeline specified in 


Table 1, the displacement volume of the casing and return water line (with 


: 


3, N Ss 
TABLE: ~BUOYAN cles OF THE SUBMERGED PIPELINE 


Buoyancy Negative Buoyancy 


i in lb per 


(6) 


Oper: rational tional — 
-loaded) 


a _ thes have a buoyancy in sea . water, By, 0! of 119 lb per ft (64 Ib per cu tt), 
-_ 7 For the total dead weights shown | in n Table 1, the net weight, ¢ or r negative buoy- — 


«6 8 (46 Ib ser cu ft to 67 lb per cu ft) in the air-filled condition and from 112 a 
Ib per ft to ‘152 Ib per ft (60 lb per cu ‘ft to 82 Ib per cu su ft) in the year a 
When entrenched and backfilled in disturbed and semi-fluid silt or clay | the 
“force of uplift on the pipeline may be considered to be equal to the weight of 
- interstitial (void or pore) water att, then ( the pipes. | if Wg be the dry density 


“4 


-¢ the sand or clay (Ib per cu ft), then (r Wa) is the weight of water in unit 


volume, in which r is the moisture content (percentage). One cu ft of pipeline — 
in the ‘sediment of (r Wa). In terms of the 
Ww w, this would be an 1b per cu ft) 


Clay 
ig [ins | | ase | | 
ae 
3 | 


SULPHUR PIPE LINE 
‘Taking representative mean values, Wy = 100 Ib pe: per cu cu ft _ r = 100% (Table 
2) for the clay, the uplift force in this sediment is thus Bg = 50 Ib per cu ft. 
It follows that the loaded pipeline would have an effective ye weight ' within * 
‘the clay of 138 Ib to 178 Ib per running ft (74 lb per cu ft to 96 Ib per cu ft), _ 
= qari. In the sand medium, with Wy = 110 lb per cu ft and r = 25% 
22 Ib per ft (102 lb per cu ft to 123 Ib 
“per cu 1 ft). figures a are » summarized in Table 3. 
is to be determined whether, in the circumstances that a hurricane 
could so agitate the backfill over and around the pipe as to make the sediment 
semi-fluid, the pipeline would tend to become unstable and rise to the surface. 
The criterion for this happening could be considered to be a state in which | 
‘the negative buoyancy of the pipeline inthe sand or clay would be less than the © 
negative buoyancy of the solids in the disturbed sediment 
If Gg be the specific gravity of the solids in the sediment it can be shown 
buoyancy, Ws, of per cu ft of ‘sediment 
1 
r+- 
in lby per cu ft. For a value of Gg 0 of 65, Eq. 2 ‘reduces 
“Tbr per cu ft. ‘The value of r for a semi-fluid sand or mud could be expected to ° 
aa the maximum shown in Table 2. Thus, _ if r ranged from m 100% to 150%, 
the negative buoyancy of the solids would be from 27.8 lb per cu ft to 20.4 Ib 
per cu ft and the conclusion must be that, so far from there being any el 


Vertical "Equilibrium in the Sediments. .—The approximate pressure of 
_—* pipeline on the sand foundation in the excavated trench near the shore- _ 

end would be from % 78 Ib per sq ft to 106 lb per sq ft that would be well within 
the bearing capacity of the sand (estimated to be about 500 Ib per sq ft). The 

sand foundation is not considered to present any instability problem except 

: insofar as it may show a tendency to develop piping under the influence of the © 


loss from the pipeline to the sediment conditions. Further 


in which 18 the maximum strength of the sediment, and Ba dimen- 


 sionless coefficient. exceeds 5.14, the pipeline? would definitely sink. 
- Adopting B = 2 for safety,8 the permissible maximum shear strength of of the 
that should be lized (for D = .44 ft) is 


in per sq ft Thus, for W Wm ranging | from 138 Ib} per ft to 5178 ger ft 


7 Theoretical Soil Mechanics, by Karl Terzaghi, John Wiley and Sons, Inc., New York, = | 
8 “Some Reemeeinn wit Engineering Considerations in Marine Pipe Line Con- 

struction,” by R. O. Reid, Proceedings, 2nd Conf. on Coastal Engrg., ‘Council on Wave 
Research, Engrg. Foundation, Berkeley, Calif..1952,p.749, 
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Ib per sq ft to 124 lb per sq ft pressure), the minimum values of T,, a -" 
to assure . safe ‘Support of the pipeline range from 46 lb per sq ft to 60 lb per ~ 
- These values apply to undisturbed sediment that is presumed to exist at 
a the bottom of the entrenchment. Fig. 5(a) shows that the contours of Tm for a 
_ 40 Ib per sq ft, 50 lb per sq ft and 60 lb per sq ft are within 7 ft of the —_ 
ment surface for about half the distance over which the sediment is cohesive, 
Between about 17,000 ft and 27,000 ft from shore, however, these contours 
sharply and extrapolation suggests that the greatest depths attained by the 
46 Ib per sqft and 60 Ib per sq ft values are about 14 ft and 18 ft, respectively. 
‘The implication then is that somewhere in this zone the pipeline will tend to- 
_ sink below the level of a prepared 7 ft ‘deep trench, If positive “sinking be 
taken as occurring when 8 = 5 in Eq. 4, then the critical range of — 
Shear strength would be 19 Ib per sq ft to 25 Ib per sq x, which is encountered — 


_ could be expected to sink here and stabilize itself in a curve codegataniiie 
"paralleling the contours of Fig. 5(a) at depths somewhere between 7 ft and 
14 ft for the lightest _— coat and fases ences 9 ft and 18 ft for the heaviest 
On the basis that the lightest weight- coat t (Table 1) would probably be used 

it is possible to investigate the stresses that might ensue from the pipeline _ 
Sagging without adequate support across the pocket of weak sediment, until 
held like a catenary by the tensions generated at the abutments of stiffer — -_ 
_ Support. It is found8, 9 9 that such a situation could not obtain without the pipe- 7 
_ ‘Hine resting on supporting foundation at about 13 ft or 14 ft beneath the mud- © 


line near the oumer of the omy In the event that the pipeline has two un- s 


Thermal the Sediments. Annual variations at 
the mudline and within the sediment below 15 ft of vater have —_ awe : 
in the Atchafalaya Bay area off the coast of Louisiana. These datal0 (Fig. 6) 7 
‘show that the overall mean temperature within the sediments in that area 7 ' 
tends” to be about 72°F. At 10 ft depth below the mudline the amplitude of the 7 
_ annual temperature fluctuation is about 3°F; near the mudline a design range 7 
of temperature over the year from 50°F to 90°F would seem to be indicated 
and should be applicable to the Grand Isie area. 
‘The fluctuations in temperature within the sediment, measurable in Fig, dl 
at the different levels, both short period (afew days) and long period (annual), 
provide a » a means? for estimating the thermal diffusivity (specific conductivity), © 
of the eotiment. The value of k, so found is 0.014 ft? per | hr, , that at agrees 


ee 9 “Coastal Engineering Problems in the Laying of a Sub-marine Liquid Sulphur Pipe-_ 7 
-_ line in the Gulf of Mexico,” by Basil W. Wilson, Tech. Report No. 187-1 (ref 58-24F) 
_ Texas Agric. and Mech. Research Foundation, College Station, Tex., December, 1958. 
_ ont 10 “Transfer Processes Operating at the Ocean Boundaries; Temperature Data; 
- Data Report (Section 1), Texas exas Agric. an and Mech. Research ich Foundation, , College | Station, 
‘Tex., December, ,1952,0 
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“quite e well with a value of 0. 011 11 £t2p per hr found _— Reid! ¢ fora different 2 area 
if Cy be the specific heat of the wet s sediment, then the thermal conductivity 
«K of the wet sediment is one 
K= = Wy C 
For Wy in Ib per cu ft, k in £t2 p per ” mi Cy in BTU per = °F, the units of K- 
will be BTU per °F per ft per hr. The specific heat Cy, can be expressed in ; 
terms of the specific heat of the dry solids Cq and the a heat of of sea — 
water Cg by the relationship. 
Laboratory determinations!2 of Cq for "samples: of the sediments show an 
_ average value of 0. 22 calories per gram per °C. Accordingly, with Cg ~ 1.0, © 
dependence | of Ce on r for a range of values of moisture content from 48% to 2 


TABLE 4.- —COE FFICIENTS OF THERMAL CONDUCTIVITY O OF SEDIMENTS; = 


EXPERIMENTAL DERIVATION 


1Degree | Conductivity K 
per °F per ft per hr) 


ust — 


ine level 


“90% at the pipeline ‘depth yields values for Cw — 0. 45 to 0.56. These to- 
_ gether \ with k = 0.014 ft2 per hr, permit evaluation of the thermal conductivity 
_K in terms of Eq. 6. Anaverage value of K = 0.71 BTU per °F per ft per hr is 
‘thus obtained from typical boring test data for the clay at the pipeline de depth 
of 6 ft to 7 ft below the mudline. 
_a _ The thermal conductivities K of samples 0 of the ‘wet sediment were later 
determined!3 py refined experimental methods and found to vary slightly 
with temperature from 1. .2 to to 1.4 for the fine sand and 0.7 to 1.0 for the soft 


_ 11 “Oceanographic Analysis of Marine Pipe Line Problems: Engineering, by x. 0. 
a Reid, Final Report, Texas Agric. and Mech. Research Foundation, College Staticn, Tex., 
12 “Thermal | ‘Conductivity of Soils, ” by C. M. Shilstone, Report, Shilstone ‘Testing 
Lab., New Orleans, La., June, 1958. 


13 “Experimental Determination of Soil Properties,” by B. J. Fluker, Tech. Report, 
Texas Engrg., Experiment Sta. Sta., College Station, Tex., August, 1958. —_—s oy 
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that the thermal properties of the pipeline components accord with Table 1, 


_ the equivalent coefficient of thermal conductivity, Ke, over the annular layers — 


ATER SEDIMENT INTERFACE TEMPERATURE, 


AS 4 4 


&S 
x 


ONSHORE 


TEMPERATURE,O, 
2 


z—PLANE 


_ TATION; (b) ADAPTATION TO THE COMPLEX Z-PLANE; (c) CON 
REPRESENTATION IN THE W-PLANE 


_ from inside to outside of the pipeline, is estimated? to be about 0.045 to 0.048 © 
It is now assumed that the heat flow from the entrenched pipeline through © 
the sediment reaches a steady state condition in which the temperature, 0s, 


 &§ Heat Flow from Pipeline to Water-Sediment Interface.—On the assumption [im 
7 
— — 
NQUTSIDE PERIPHERAL PIPE-SEDIMENT — 
= 
7.—HEAT FLOW FROM PIPELINE TO MUDLINE; (a) SCHEMATIC REPRESEN- 
FORMAL 


i 
at the mudline remains constant. of the temperature, 62, ‘pre- 
-vailing at the external periphery of the casing pipe weight-coat, may con- 
= | veniently be analyzed with the aid of the theory of complex variables and a 


h use of conformal mapping. . Thus, the pipeline of Fig. ‘7(a) may be schemati- 
-cized in the z-plane of Fig. To) by a cisele of radius R with uniform — 
= perature @9 at height d above the x- axis, along which the temperature is uni- 

formly 6g. This maps conformally into the w-plane as two circles concentric 
with the origin of the u, v-axis.14 The boundary at constant 62 isa circle ¢ of 
‘TABLE 5. _—PERIPHERAL TEMPERATURES AT THE PIPE -SEDIMENT 7 
INTERFACE UNDER PIPELINE OPERATING CONDITIONS» 


re 
Conducti- Pipe 
vity of Mud] Surface 
K, in BTU in °F 


per ftper | 


| | 137 109 


ing 
Condition 


Laying 
Condition 


radius Ro, that at constant 6, acircleof radius Ry. The flow of heat, per unit 


time, Q, through the annular layer (the sediment) in the w-plane is _— <~ 
7 


oe 


14 Dictionary of Conformal Re resentations, by H. Kober, Dover Publications 
New ‘Yor » 1957. 
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in which D and 1D, are, respectively, the diameters oft the casing pipe weight-— 


fae gooey and of the hot water line (inside). Consequently, on equating 


+ (2) 


(1) 
ie Ke = = - 0.048 BTU per oF per ft per mee 6. 7 ft, ‘D: =2 R = 1.44 ft, Dy, 
= 0.59 ft and K values adapted from Table 4, the external peripheral temper- 
_ atures of the pipeline, 62, are assessed in Table 5 for extreme winter and 
summer conditions. For a mudline temperature ‘of 80°F, the distribution of 
_ temperature and lines of heat flow throughout the sediments follow the patterns 7 
Laying and Prestressing the Outer Casing Pipe.—Having regard to the fact 
~ that the pipeline in the as-laid condition (filled with air) will exert on the 
foundation an average pressure of 59 to 87 lb per sq ft, the question arises 
whether the prestressing ‘of the casing pipe can be accomplished as desired, q 
The conditions at the shore-end where the foundation is mainly fine sand are > 
In Fig. 9(a) the pull P, mn the casing pipe is shown resisted by a distributed 
on frictional force q per unit length, which is assumed independent of the elon- 
gation of the pipeline under the pull, Taking the origin at the shore of an x- 
axis: along th the pipeline, the tensile stress, Ox, ‘developing in the casing at any 


‘distance x, will be 


) 


* (13b) 

in which Uc is the coefficient of friction of the casing on the sand, We the 
negative buoyancy per unit length of the pipeline in sea water, and $ the angle © 
of internal friction for submerged sand. For We = 85 Ib per ft to 125 lb per. 


=— 


— 
e heat, flowing trom the hot water to the outer periphery of 
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ft, d= = 20° | is found to be from 31 Ib pe per ft to 45 lb per ft and the oes ] 
a (Fig. 9(a)) in which the non from the “m0 Po= = 200, 000 1 Ib vanishes La 


therefore from 6500 ft to 4500 ft. 


Surface 


Surface 


=? 


‘ ‘resistance q 1 will pe on the ochasten of the clay and will not in general be 
independent of the elongation of the pipeline under the pull. The situation is : 


| 
* 
4 
4 — 
7 
— 
FIG. 9.—_PRESTRESSING OF THE CASING PIPE BY APPLICATION OF END PULLS» 
— SHORE-END (b) OFF-SHORE-END OF THE PIPELINE 
— 
— 


"represented in Fig. 9(b) in which the origin is again taken at the end of the 
pipe. Over a distance a from the origin, at which the elongation will be large, 7 ial 
the resistance q per unit length is assumed to be directly — to the — 


maximum shear — 


assu ming th the pipeline t to be amend by the soft clay at the mudline . Beyond 
= a, the resistance q wit be some function of the elongation, u, the — 


clay ‘(moisture content 0%) from the pipeline site, a test was: 
"devised, which is self descriptive in Fig. 10(a). Extraction of the model pipe 
yielded the test results of Fig. 10(b) from which a reasonable value of ¥ would. 
be about 4300 Ib per sq ft. The ultimate shear stress developed was of the — 
a order of 0.14 psi to 0.15 psi which is in fair agreement with the cohesion 7 
fee Mohr circle stress analyses of typical data from the borings 
fae 11). The test was decisive in indicating that the influence of Eq. 15 would 
over buta distance b (Fig. 9(b)) which, relative to a, could” 


5 


The tensile. stress from the end pull P, may thus be expressed as 


yt 
In this, Tm ,isa a function of x, as indicated in Fig. 12(a), which is derived from 
«Fig. 5. The integral of T,, with respect to x is represented then by the area 
~—_ the full-line curve of Fig. 12(a) for undisturbed sediment or the dotted- 
7 line curve for disturbed clay. The distance a from the mine- -end of the pipe- 7 
line over which the | prestress 0x declines to zero is readily found from Eq. 
16 by trial and error and graphical integration of the areas. For undisturbed 


z sediment ‘it is S only about 1400 ft and, for hah that has been aatiiaie — dis- 


may be evaluated from the basic stress- -strain relationship 


in which E is the modulus of elasticity in tension of the steel. In the case of 


the shore- end of the casing pipe, » Eqs. 12 and 17 yield 


which, for. E = = 30 108 Ib Ib per ‘sq 


TE A.’ 
indicates an Up of from 1.3 ft 0. tt ft to the 
estimates of a of 6500 ft and 4500 ft, respectively, cited » the 
_ case of we mine-end of the casing pipe, Eqs. 16 and 17 give — 


Eq. (16 gives, for x x= 
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which, , if 7; Tm conforms to the strength a disturbed sediment (a = = 9300 ft), ; 
‘results in an elongation of 1.7 ft, or, if the sediment is essentially undisturbed 
(@=1400ft), only0.24f. 

_ These computations suggest then that any effective prestressing of . 
casing be unattainable because of the high cumulative resistance 
of the sediments along the pipeline. (During laying of the pipeline in June, — 
1959, end pulls from 260 ,000 Ib to 300,000 Ib were applied to the casing and / 

. two external water lines (roughly the same stress condition for the casing — g 
and return water line treated in the text. ). Measured expansion at the shore- _ 

end was only 7 in. (prediction 10 in. to 15 in.), at the mine-end | only 2.5 ft 
(prediction 1.7 ft), instead of the hoped-for 12 ft.) End- pulling the casing pipe 
would leave the central portion entirely unstressed and, after anchoring of 
the ends, slow yield of the sediments would ultimately allow the non-uniform _ 
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13. 13.— FIXING OF THE INTERNAL LINES IN THE CASING PIPE; SCHEMATIC 


- prestress to ‘equilize at a very low value of stress along the entire ‘Pipeline, 

tending, essentially, to defeat the purpose of the prestress. 
7, Fixing of the Sulphur and Hot Water Lines.- —The fixing of the inner lines is 

portrayed schematically in Fig. 13(a) with reference to an origin taken at a : , 
offshore end of the pipeline. _ Expansion: of the inner lines generates a small 
amount of friction in the rollers which carry them within the casing pipe; this: q 
is assumed constant and designated q per unit length. The cumulative resist- 

ance forces push on the offshore anchorage, in consequence o of which a com- 


_—_ stress develops in the inner lines that at any point x will be 


4 
— Wate? & Sulphur Lines 
ag! 
| 
if 


= which Li is the total length c of the casing pipe and A, and _— are ‘the sete | 
sectional areas of the sulphur and hot-water lines, respectively. (The negative 
sign is to indicate compressive stress. ) Taking q= Uy (We + Wh) 


cient friction of of roller bearings 
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FIG. 
We and Wy) are the dead weights (in air) of hot-water-loaded hot- 
water lines respectively, we find, for the particulars given in Table 1, with» 


= 0.03, a linear distribution of stress from 5400 lbs per | sq in, atx = 0 to 
zero at x= shown in n Fig. 


> 
> 4 
— 


The elongation of the inner lines under _ temperature e rise, 6 9, from 80°F i. 
to 190°F may be shown to 
in n which oj is the coefficient of the steel of the inner 
(see Table 1) and is computed to be 26.7 ft. (It is understood that actual ex- 
pansion was of the order of 27ft.) 
Stresses in the Pipelines under Normal Operation.—Under 
ating conditions of the pipeline with the inner lines fastened to the casing, the a 
ends anchored and the pipeline buried, it is assumed that the anchorages do 
not permit any sensible breathing or elongation of the pipeline. “Accordingly, — 
7 the tendencies towards expansion of the inner lines under the rise, 0; from 
the fixing temperature of 190° F tothe operating temperatures of 300°F : 20° 
= are thwarted by the end-forces P of the anchorages (Fig. 14(a)) and result no 
thermal compressive stress throughout the line distributed as shown in Fig. 
(14(b) with a variation from 35,600 psi at the mine-end to 20,900 psi at the 
-shore— end (inclusive | of the fixing prestress). The stress Oj from these end = 


constraints is — from 

The | casing p pipe e (and return water line), however, will also be subject toa 
temperature Tise, Oc, from the fixing temperature of 80°F to the values" 
shown in Table 5, thereby acquiring stresses, do, , given by Oc=- 
which | OM is the coefficient of linear expansion of the « casing steel (Table 1). a 
Under extreme summer conditions, (Table 5), thermal compressive stress in 
the ¢ —— — vary from about 9900 psi at the mine-end to — psi at the 


casing pipe be quite small (Fig. 
ee. In the event of an unavoidable shut-down of operations in water, resulting _ 
ina drop of all line-temperatures to 65°F (Fig. 6), the induced thermal tensile 
stress in the inner lines would vary from _ about 23,700 psi at the mine-end 
to 29,100 psi at the shore end (Fig. 14(c))} 
Longitudinal Stability of the Pipeline under Thermal Forces.—With end 
anchorages inhibiting expansion of the pipeline the latter would be prone to 
- buckle unless prevented by the lateral restraints of soil forces. The condition | 
is posed in Fig. 15(a) in the rather hypothetical case of half-wave ‘buckling of 
the pipeline over its entire length L. The deflections y y ‘under the compression 
of the forces Pj and P,, pertaining to the inner lines and casing pipe respec- 
- tively are, , of course, not necessarily vertical. Buckling is resisted by normal 
- soil forces p per unit length, considered proportional to the e deflection. The 
modulus E, appearing in this 


some departure from constancy i in sediment as will be 
_ cussed later. The end conditions of the line are assumed to | approximate 
hinges because any er end moments will be of small eaciaianaid due to the length 


ofthe lines | 
lines 


- 
= 
hes, 
— 
_ 
af 
» 
| 4 
i stress there will be actual elongation of the * f 
pth L. Thermal stress willtend tobe relieved 
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by the buckling so that ‘the forces Po _ cé will | diminish | as the deflection 


y increases. If 6 be the increase in length of the x pipeline resulting from the e 
buckling, | then the relief of force will amount to E A, 6/L for the casing and 
EAj 6/L for the inner lines” (in which Ag and Aj are respective cross- 
sectional areas), so that the work done by the induced forces w will 
(Pj -2 E A; 6/L)é on the inner lines and (Po -2 §/L)é on the casing 
= . Following Der Den 1 Hartog, 15 ‘it may be sh shown 7 the ae 6: is given by 


dx 8 L 


in: which I, 2 and I, are the section ‘of the casing pipe and the internal 


ey given by Den Hartog!5 for the case of a hinged member under 
‘compression from constant end forces P is the same as Eq. 29, but without © 
- the second term on the left hand side. _ critical buckling load Per: as given | 
y Eq. 29, is | 
/L E A 
7 
“The third term on the aie hand side of Eq. 30 arises ses only in the case of 
compression from thermally induced forces and represents additional com- - 
pression that the pipeline could withstand because of the tendency towards © 
mediate stress relief on buckling. If it can be shown that the forces P (Eq. 
St = computed from just the first two two terms of Eq. 30, then — 1 
Pood: pipeline will be quite definitely _ stable against any condition of buckling. | 
_ As shown by Den Hartog, the liklihood of a very long pipeline buckling into — 
a single h: half-wave is quite remote. The tendency is for the line to buckle into 


a series of alternating half-waves, Fig. 15(b), of length 1. The condition in fact | a 


15 Advanced of Materials, by J. P. Dea MoGraw-Hill Book ¢ Co., Inc. Inc., 
York York, 
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‘the energy of *ibending;, that are, intnat order, 
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IG. 1 5.—STABILITY OF THE PIPELINE AGAINST BUCKLING UNDER 

_ THERMAL STRESS; (a) SCHEMATIC REPRESENTATION; (b) | 

MULTI-HALF-WAVE BUCKLING; (c) DIMENSIONLESS | 
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SULPHUR PIPELINE 
- that a very long member buckles into a half-wave of specific length 1 is de- 


7 fined by the minimum value of P,., which will just make that length unstable. * 
This condition, , represented | by the point M in 
VE, EI E 
— 
which r ‘result is | given also by xy Hetenyilé by use use of another r method ee 4 
_ Evaluation of the Soil Modulus, E .—The use of a soil modulus in the re-— 
lationship of Eq. 24 is now firmly established in foundation engineering. 17,18. 1 
_ However, it is generally conceded that Eg usually increases with depth and 
at a given depth it decreases as the deflection y increases,17,18 
value of Eg to be adoptedin Eq. 31 must therefore take cognizance of the oer - 
of the pipeline and of the remolded nature of the material Bea 
_ As remarked previously, the cohesive sediments covering a — wide 
area of the coastal region west of the Mississippi delta are very similar. 5 
comprise mainly -montmorillonitic clay with inclusions of Illite 
‘Kaolinite. 20 For such material F. D. Gaul,19 A.M. ASCE, has found a value 
Eg = 6900 ib per sq ft to be applicable in laboratory experiments of . 
laterally loaded pile; this is thought to be rather high. Fortunately, test pa 
exist from field measurements of a laterally loaded pile in the neighborhood - 
of Timbalier Island about 4 miles off the Louisiana coast, some 25 miles — 
west of Grand Isle. These data, comprising measurements of p and y (Eq. 
24), have been correlated with standard laboratory tests of the sediments ob- 
tained from borings made on the site by B. McClelland, M. ASCE, and J. a 
Focht, 18 M. ASCE. These authors show that the relationships of p to y for dif- 7 
ferent depths follow trends that are similar to the relationships of deviator 
stress, 0, to strain, ¢, in the quick triaxial shear-strength tests. According — 


to their correlation, w which they hold to be generally - applicable to cohesive 
materials, the soil modulus can be from 
| 
(82) 


oil 16 Beams on Elastic Foundations, by M. Hetenyi, Univ. .of Michigan Press, Ann Arbor, | 
, 19%, p. 15 
’«Non-Dimensional Solutions for Laterally y Loaded Piles with Soil Modulus Assumed 
; Proportional to Depth,” by L. C. Reese and H. Matlock, Proceedings, 8th Texas Conf. 
on Soil Mechanics and Foundation Engrg., Special Publication No. 29, Bur. of Engrg. 
Research, Univ. of Texas, Austin, Tex.,1957. 
_ 18 “Soil Modulus for Laterally Loaded Piles,” by B. McClelland and J. A. Focht, Jr. “” 
Proceedi s, ASCE, Vol. 82, No. SM 4, October, 1956. eee ee eee 
a "Movel Study of a Laterally Loaded Pile,” by R. D. Gaul, Proceedings, ASCE, 
84,No.SM1, February, 1958. 
20 “The Clay Minerals of Recent Marine Sediments to the West of the Mississippi 
al Delta,” by R. F. McAllister, thesis presented to the Agric. and Mech. College of Texas, 
in College Station, Tex., in 1958, in partial fulfillment of the requirements for for the de- 
' 21 “An Investigation of Lateral Loads on a Test Pile,” by A L. Parrack, Tech. Re-- 
port Proj. 31, Texas Mech. Research Foundation, Station, Tex., 1952, 
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Here, » dentate. because g/e is not a linear function and varies with depth, 
Eg must generally be a function of both depth and strain. From the nature of 


‘the g -¢ relationships obtained from the boring tests along the pipeline site, 


would that can be curve-fitted by an equation of of 
form 


in which C is a constant and n a numerical exponent. Evidence of this is pre-— 
sented in Fig. 16 for typical stress-strain curves taken at random, The values 
of C and n necessary to obtain the best possible curve-fits are plotted in Fig. 
17 against distance along the pipeline. The conditions of immediate interest 
concern the values of C and n for the remolded clay at the pipeline depth, ’ 
approximated by the test data at 5. 5 ft below the mudline in Fig. 17(c). The 
constant C is found to vary from 15 to 37; that of n from about 1.9 to 2.3. It a 
» reasonable to adopt mean values of C = 25, n= 2. 0. 


for the values of C and n quoted 


The a, of Eq. 35 is shown in Fig. 18 and reflects th the installa of E 3 on 
the: deflection. ‘Because the ultimate 2 strength of the cohesive material 
reached in many cases when the strain exceeds about 6%, a mean value of Eg 
up to this limit of strain would seem to be a satisfactory criterion for adop- 


tion. We take therefore, accordance with» Fig. 18, Eg psi or 2450 Ib 


This value may be compared with the result of the extraction test (Fig. 10) 


for which 7 was found to be 4300 Ib per sq ft in the relationship of Eq. 15. _ 7 
Restating Eqs. 15 and 24 in terms: shear stress, T, and normal stress, 


-(36b) 


For u: u= =y , and ; a a shear str strength ¢ taker 
= 27), the ratio p/q yields 


* 
this one might thus expect be (Q/rx 4300) or about 2740 be per 


sq ft, which is in fair agreement with | the value of 2450 lb per sq! ft given pre- an 
_ The Critical Buckling Load.—The pipeline has the choice of buckling in any > 


- direction. If the return water line is strapped below the casing pipe, the second 7 
- moment of area of a cross-section of the system about a vertical axis will be 
less than that about a horizontal axis and will therefore favor lateral buckling © 

q of the pipeline. Lateral buckling would seem | to be favored also by two. other 7 
_ circumstances, namely, that vertical buckling would require penetration into — 

- Y atiffer material at greater depth in the case of a downward half-wave at the 

= time that an adjacent upward half-wave would necessitate the lifting of © 


the Pipeline ‘against its own submerged weight as w well as the overburden soil 
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ry Assuming lateral buckling then, pang second moment of area @ of the pipe- pipe- 


sq ft and E = 30 x 108 psi, the critical load | Por is ‘found to be 956 ,000 lb, and 

_ the critical half-wave length 1 ; = 44 feet. (For such a hinged length, without — 
~ lateral ‘support, the Euler critical load would be 72 E 1/1? or 476,000 lb.) Be- a 
cause actual ees P from the induced thermal stresses at the mine-end of the 
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is (4 Ay + we find ind (for 04 = 35,600 psi, Aj = 12, 07 sq in. 
= 9900 psi and A¢ = 19.73 - ‘in.) P = 625,000 Ib. It seems, therefore, that ‘the 


Effect of the Heat on the Sediments in Contact with the Pipeline. —All pre- 
- vious considerations of this paper up to this point have been based on the sup- 
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- position that the foundation behaves in accordance with the properties deter-— 
~ mined from the test borings. The remaining important consideration concerns 
the influence of the heat flow from the pipeline on the properties and behavior 
the possibility can be envisaged that in the fine sand sediments at the 
shore-end of the pipeline convection currents may be set \ up in the pore water 
which, through reduction of intergranular pressures, could lead to semi- or 
quicksand conditions and | promote limited settlement of the line. The settle- 
_ ment could presumably continue until the convection currents became diffused 
to ineffective ) magnitudes at the greater Gopths. The amount of such settlement _ 
would seem t to be unpredictable in the present state of our knowledge. 7 
In the soft grey clay sediments prevailing along two-thirds of the length of _ 
_ the line, the chief concern might be that thermal osmosis will result in a mi-— 
7 _ gration of interstitial water from the heated zone, leading to a reduced water 
content in the clay in contact with the pipeline. This, it is thought, could pro- 0 
- some degree of desiccation and shrinkage of the medium that might in- . 
_ duce tensile stresses in the casing pipe and lead to subsidence of the pipeline. 
Regarding the imposition of tensile stresses in the casing pipe from 
shrinkage of the clay medium it would seem that the process of expulsion of 
void water from the clay in the neighborhood of the pipeline along the heat — 
flow lines would be a fairly gradual process. The developing shrinkage stresses 
should therefore have time to disperse through the same process referred _ : 
~ in Section 8 in respect of the prestress, namely, that of a continuing slippage 
- or creep taking place at the ultimate shear strength of the heated clay, a 
= reduced from the liquid to the plastic limit. It might also be thought 
‘aa on actual tensile stress imposed on the casing pipe as result of clay _ 
ay might even be beneficial in that it would help to provide | a condition © 
tension to offset the compressive stress. in the line and thereby achieve 
naturally the prestress condition that is sought from preliminary endpulling 
= the casing pipe. It seems that the poeniom of the effect of clay shrinkage on 


; = of void water by the thermal osmosis would lead to a 1 consolidation co z 
the clay in reverse along the heat flow lines, Because upwelling of material 
from below the pipeline cannot be expected, the line must ‘necessarily sink, 
_ The degree of settlement likely to ensue, however, must obviously be limited - 
y by the increasing strength of the medium encountered beneath the pipe and the - 
— at some stage for the clay thereto reach stability : from consolidation 
brought about by progressive bonding. Again, it does not seem possible to 


amount settlement ‘that t might be involved. 


CONCLUSIONS 


P 2 Laying of the 7-mile sulphur pipeline at a depth of 6 ft to 7ft below the 


-mudiine off the eastern end of Grand Isle, Louisiana, its burial in 
very fine sand and silt over a distance of 11,000 ft from the shore and in ex- 


tremely soft grey clay over the remaining distance out to the mine structures. 7 
= In these media, the lightest density of weight-coat for pipeline protection, 

(135 Ib per cu ft in the range of densities treated in this paper), , would seem 
to be be best from the point of view of relieving f 
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of "weakness without danger of flotation of the line in the | event of agitation of | 
the backfill by hurricane waves and flood. = 
_ The sand sediments are expected to provide adequate support for the line 
a settlement is possible as result of convection currents induced 
in the pore water by heat flow from the pipeline through the soil. In the co- | 
clay further offshore, the pipeline shouldbe generally stable at 
bottom level except over a distance between 17,500 ft and 23,000 ft from the 
_ coast, where it will almost certainly sink toa depth somewhere between 7 ft b 
and 14 ft below the mudline in a pocket of extremely weak sediment. The 
gradients involved in this sag would be very slight and combined tension and © 
flexural stress “resulting from the depression should not exceed about 4300 
psi. Heat losses from the pipeline may induce other limited subsidence in the © 
_ Clay as a whole, asa consequence of thermal osmosis and partial desiccation 
7 of the medium in contact with the pipe. Progressive bonding and consolidation 
round the pipe, however, should ultimately arrest this movement. 
_ 4 Analysis of the heat flow from the pipeline through the soil suggests that 
casing pipe temperatures can range from 74°F to 121°F (Table 5). These 
: temperatures affect the stress condition in the casing pipe, but not in the — 
lines, if ends are anchored. The highest thermal stress anywhere is 35, 600 
psi compression at the _mine- of the inner lines under operational 
tions; in the casing, maximum thermal stress should not exceed 10,000 psi a 
compression. It is shown that even the weakened soil forces to be expected of 
: disturbed backfill should be adequate to to me pipeline 


_ It is a final conclusion of this study that uniform onenen of the c casing 
pipe would inevitably be frustrated iby the cumulative soil resistance mobilized 


any ‘anchorages unnecessary, because the friction is sufficiently great cumu- 
latively to be able to ‘restrict end- expansions to small final 


thermal stress at the ae that it would afford. The _— nature of this, 


— 
In certain areas of this study ‘the valuable advice and assistance of R. oO. 
Reid, R. E. Schiller, M. ASCE, and S. J. Buchanan, F, ASCE, were of great 


express his appreciation to the Freeport Sulphur Company, for all 


study was undertaken, for permission to publish this material. peer 
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Ag = = cross-sectional area a f the steel a the e sulphur line; 


= buoyaney of ¢ pipeline entrenched in 1 the si sediments; 


dimension of length (ie. 
= constant of proportionality, Eq. 33; 


ue specific heat of dry solids in sediments; 

heat. of sea water; 

specific heat of wet sediments; 

= outside of casing pipe, inclusive of welght-coat;, 


diameter of hot. wer line; 
depth below mudline, 7; 


modulus of elasticity 0 of steel; 
= modulus of f soil 1 foundation, Eq. 2 
= specific gravity of sediment eatin: 
= film coefficient of heat conduction, (Table 2); 
second moment of area of steel cross-sections for all components of © 


a pipeline, about vertical axis; 
ey value of I for casing oaien and return water line; Sor 
= value of I for concentric : inner lines; 
K = = coefficient of thermal of wet 


= thermal diffusivity (specific ‘conductivity of heat) of wet 


= length of pipeline; 
= critical half-wave length of buckling; 


= total restraining force acting on end of pipeline; 
value of P applicable to and return 


Por bucking load for | 


value of P applicable. to inner lines; 
lateral soil force on pipeline per ‘unit length; 


a quantity of of heat flow per unit distance per unit time; 
per unit length; 
= radius of casing pipe to exterior of weight-coat; 
= radiu e toe 
R xterior of weight-coat; 


A. cross-sectional area | r pipelines (= An + As); 
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_ SULPHUR PIPE LINE 
= moisture content of sediment; — 


= = degree of saturation of sediment, (Table 4) )s 
longitudinal elongation of pipeline; 


= value ‘of u at origin 2 


= = value. of u at any point x ‘from origin; aa 
= negative b oyaney of pipe li ine ins sea wa it 
negat uoyancy f pipeline in sea water er per u uni t length; 


= density of of (ie perc cu 


= = density of solids per unit volume of sediment (Ib per cu u ft); 


= variable lateral on 


a; = coefficient of linear steel of inner lines; 


= dimensionless strain, ‘Eq. 335, 


difference, rise or fall in temperature; 


= rise or fall of temperature of casing | pipe relative to its to its fixing tem-— 


peratures 


« a _ = rise or fall of temperature of inner inner r lines nes relative to to their fixing tem- 


= temperature of water-sediment interface; 


of hot v t 
temperature of te) wa 
0 om temperature at exterior surface of casing pipe weight-coat; 


Uo = coefficient of friction of casing pipe on sand sediment; 


effective coefficient of of roller bearings inner lines 


within casing pipe; pO 
: 
universal constant, 3. 1416; 


"aa 4 = = normal ineeie in sediment, Eq. 32; 
= _ = longitudinal thermal stress in casing | pipe and return water line; 


ong thermal stress in inner lines; 


oe "= longitudinal sts otress in pipeline | at any point distant x x from origin; 


stress in cohesive sediments; 


ts = ultimate shear strength of clay sediment; and 


7. : 
uy 
ensity of sediment (lb per cu ft); 
x = ‘ 
— return water 
constant of proportionality, Eq. 15; — 
— 
— 
| 
| 
angle of internal friction of sediments in submerged state. | 
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A a was made to determine the effect of non- -hydrostatic consolidation 7 


é the cohesion and friction in two saturated clays using the CFS-test proce 
dure. The results showed that cohesion is unaffected by non-hydrostatic con ; 
ae , The value of friction at zero CFS-test strain is increased by non- 7 

- hydrostatic consolidation but the variation of friction with strain is the same 
in ahydrostatically consolidated sample. In addition a time transfer of co- 


hesion to friction was observed. = ae 


INTRODUCTION | 

Notation. —The better evel symbols adopted for use inthis paper. are » defined where 
they first appear, in the illustration or in the text, and are ae alphabeti-— 
cally, for convenience or reference, in the Appendix. 
Review of CFS-test.—A previous paper3 reported the first published 


a of the basic research on the components of soil strength being ——- in 


Soil } Mechanics Research Laboratory of the ‘University « of Florid 
erence reports the development of atriaxial compression test, meieet the CFS 


an abbreviation of cohesion-friction-strain) -test, designedto permit the com- 


Note.—Discussion open until January 1, 1962. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 

ot the copyrighted Journal of the Soil Mechanics and Foundations Division, Proceedings _ 
of the American Society of Civil Engineers, Vol. 87, No. SM 4, aan, 1961 
2 Asst. Prof. of Civ. Engrg., Univ. of Florida, Gainesville, Fla. . 


2 Graduate Asst., Civ. Engrg. Dept., Univ. of Florida, Gainesville, Fila. 
3 “An Experimental Study of the Development of Cohesion and Friction with . ala 
s in Saturated Cohesive Soils,” by J. H. Schmertmann and J. O. Osterberg, June 
1960 Research Conf. on Shear Strength of Cohesive Soils, ASCE, duly, — ‘ul Sad 
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putation of cohesion a and at strain — during compression and 
_ thereby establish curves of the variation of cohesion and friction with strain. — 
Since an understanding of the definitions used for cohesion and frictionand the _ 
techniques of the CFS-testare necessary in order to follow an work reported - 
herein, they are briefly reviewed in this introduction: ir ve Pa 
The angle of internal friction ¢¢ atany strain is the angle for which the 
_ gent is the ratio of the change in shear stress to the change in normal effec-_ a 
tive stress occurring» on the plane of Mohr envelope tangency at that strain, 
during a stress change occurring without significant change in soil structure. — : 
‘The cohesion ce of a soil, at any strain €, is the shear stress developed on — , 
the plane of Mohr envelope tangency at that strain, if the effective stress on 
that plane could be reduced to zero with proportionally the same change in soil , 
structure as occurs in the determination of de. ‘a ear 
_ For saturated soils, the effective stress ata point is defined as the total 
stress less the pore water pressure measured by a piezometer at that point. _ 
The plane of Mohr r envelope tangency is the plane of maximum stress obliquity 
for the frictional component of strength. The expression ° “without significant 
change in soil structure” is, as of 1961, | interpreted to be a void ratio change 
of less than 1% in the same soil, caused solely by a change in pore pressure. 
From: the foregoing, one may see that the cohesion and friction definitions © 
used herein have a physical significance which one can “feel” and work with | 
inthe laboratory. The expressions “cohesion” and “friction” are retained here- 
infor continuity withthe previous work. 3 However, other expressions for these 
components, such as “pore pressure insensitive” and “pore pressure sensi- 
tive,” may eventually be preferred. There is no claim that these are basic 
mechanical properties of a soil. The definitions were chosen with the hope 
they will eventually prove useful in engineering application and thereby | 
justify their existence. The writers are optimistic that this is the case. This 7 
work presents only part o of the results from a research program active at the 
University of Florida, under the direction of the senior writer, and thoveiese 
z= only a part of the experimental evidence that will eventually be pre-- 
_ sented forthe evaluation of the engineering significance of the | CFS- -test defini- 
_ Fig. 1 helps toexplain the CFS- -test technique. An axial deviator stress, 0, 
- imposed on a soil ‘Specimen inthe triaxial chamber by subjecting the speci- 
mentoaconstant rate . of compressive strain. During this strain the pore water 
pressure is externally controlled to maintain a constant value of major princi- 
pal effective stress, 01, on any horizontal plane : suchas 1-1 in Fig. 1 (a). The 


- control is by periodic manual adjustments of p pore e water | pressure, u, to: main- 

in y which ! is the constant cell pressure and 3 is held constant at 
selected value. Thus 


One of the pre-selected 01 values is the  r)nigh | indicated in Fig. 1(b).. Two 


‘Successive deviator stress-strain points are obtained at (01) high, 
as represented by points a and b. ‘The pore pressure is then increased to re- 
duced he to the pre-selected value of (91)1ow- . In the work reported herein the 
~— (6)nie high Value was 95% (0), imposed during consolidation, and the (61) jow was 


x Fie The e increased pore pressure reduces the ‘frictional component of 


Pon 


2: 

> 

: 
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Strength and thus reduces the deviator stress the speciman can support at that 
strain. Because the test is strain controlled, the ‘result ij is a measured reduc- — 

tion” in deviator stress instead of an increase in strain. After sufficient time 

™ for % equilibrium at the lower value, two additional stress-strain data points 


Hydrostatic 


(a TRIAYIAL STRESSES 


- actual stress-strain path 


are obtained, as represented by points c c and d. ‘The po. pore pressure is then re- 
duced to raise Gy tothe (03) high value which results in increased friction and 
measured , deviator stress: Points e and | f are obtained to define two stress- 


strain curves, one for each 0] as shown by the dashed lines in | Fig. 1 (b). 


sM4 

= 
ir 

— 


At any value of strain, such as € in Fig. 1 (b), a Mohr circle c: can be drawn 
representing the stress conditions for the ‘Specimen at each 0, value. Thetwo 
circles” obtained are shown in Fig. 1 (c). The angle made by the tangent line 
1-2 withthe 0 axis is ¢¢, which is then usedin Eq. 3 to obtain c,. Eq. 3 is de- 


rived from the shown in 1(c)s 


The values of de ond Ce — are 
of cohesion- ~friction- strain behavior for the soil specimen tested. _ ai 

2 Purpose. .-In any field situation the soil is actually consolidated under non- 
‘hydrostatic stress conditions rather thanthe convenient hydrostatic consolida- 
‘tion usually employed inthe laboratory triaxial test. It is therefore important 
_toknow if non- 1-hydrostatic consolidation results in in materially differ different st subse- 


TABLE 1, 1, _—SUMMARY OF RESEARCH PROGRAM co 


7/1.0, 1.2, 1.4, 1.6 Standard (see pro- 
cedures) 
1. 0,1. 2, 1.4, 1.6 mented 
13.65 0, tandar 
K 


"principal effective stress, 01, seemed much more important than 
_ structure in determining cohesion behavior. _ However, the possible importance 
_ of anisotropic structure was not studied. - Such a | study was one one of the purposes > 
a cohesion is anisotropic due to anisotropy of soil structure, then t the work 


» AL \. Casagrande, F. . ASCE, and N. Carrill indicates? that fora constant value 


the cohesion developed on the potential failure planes should vary measurably 


_ with significant variations inthe degree of structuralanisotropy. In this study, 


4 “Shear Failure of Anisotropic Materials,” by A. Casagrande and N. Carrillo, Con— 
tributions to Soil Mechanjcs 1941-1953, Boston Soc. of Civ. Engrs., p. 122. 
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quent strength behavior when compared to hyd 
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____ The writers were especially interested inf |" 
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atest s series s by keeping initial soil structure, e, 04, the plane of 0} SS 
- andvoid ratio constant within the series. The degree of anisotropy of the soil 7 
structure was varied by changing the 01/03 ratio to which the specimens joni 
consolidated prior to CFS-testing. This ratio was varied from 1.0 to 1.6 and 
it is assumed that significant variation in structural anisotropy was produced. ; 
Comparison of the results of such tests were then used to judge the possible = = 
anisotropic behavior of the cohesion component of soilstrength, wea 
Scope.—Table 1 summarizes the scope of this investigation this investi- 
gation two machine-remolded, near-saturated clays were used. Kolinite 
chosen as the main clay | because previous showed it had strength 
properties similar to a variety of clays natural ‘composition and had advan- 
tages of (a) additional structural sensitivity, in spite of machine remolding, 
Cb) known minerology, and (c) relatively high permeability which aids the con- 
venience of CFS- -testing. be Four of the five series used duplicate specimens of 
‘koalinite sothat these series may be compared directly. To evaluate possible 
effects of mineralogy and grain size distribution, Series II was performed in 
. an identical manner to Series I except that Boston blue clay was used. This 
= is an iliite-type clay with considerable silt. It is important to note that both 
_ kaolinite and illite are platey clay minerals and the results presented herein 
may be limited to clays with similar r type clay minerals. 
Series I and Ilare duplicates, except that the stresses in area approximately 
_ double those of I. The stress was doubled to determine whether the magnitude ~ 
_ofstress was a significant variable. The values" of 1. -80 kg per sq cm and 3.65 
kg per ‘sq cm were chosen for convenience. The ratio (04 /03), was limited 


to i. 6 of excessive compressive strain at higher values. 


- Table 2 presents additional detail about | test and and specimen cc 2n conditions. 


—— 


Equipment Used.—The triaxial and anisotropic loading equipment was 
dev eloped by the Norwegian Geotechnical Institute. A complete description of 
the equipment shown in Fig. 2 is available. 9 The NGI cell piston design was 

modified to accommodate a 1/4- in. . piston running in ball-bearing sleeves.6 
> An SR-4 loadcell was developed to replace the use of a proving ring. - ‘This 
cell hasa compression of 0.006 mm perkg versus 0.090 mm per kg fora prov- 
ing ring of the same load capacity (30 kg). This increase in rigidity enabled 
the tests to be run at a more constant rate of strain. The design of the load 


5 “Triaxial Equipment Developed at the Norwegian Geotechnical Institute,” by A. 
Andresen, J. Bjerrum, E. DiBiago, and B. Kjaernsli, Norwegian Geotech. Inst. Publica- 
tions No. 21,08lo, 1957, 
6 “Prestress induced in Consolidated-Quick Triaxial Tests,” by A, Casagrande and 
ey Wilson, Proceedings, 3d Internatl. Conf. on Soil Mechanics ‘and Foundation Engrg., 


7 “An Experimental Study of the Effect of Anisotropical Consolidation on on the Cohesion 
and Friction in Saturated Clay,” by J. R. Hall, Jr., thesis presented to the Univ. of - 
Florida, at Gainesville, Fla., in June, tates partial fultillment of the requirements 


for the degree of of Master of Science. 


— 
— 
| 
— 
» — 
ee — 
— 


ER ANISOTROPIC CONSOLIDATION 


-TEST AFT 


G CFS 


oF 


2. SPECIMEN DURIN 


|| 


“August, 19 1961 


= preparation of duplicate remolded specimens with ahigh degree of sat-_ 
uration and structural duplication, a“Vac-Aire” extruder was used. 8 A helical 
_ —particle- -alinement structure results from the use of this extruder, and this, _ 
unfortunately, complicates the mental picture of the structural effect of sub- — 
sequent non-hydrostatic consolidation. However, non-hydrostatic consolida- 
tion should develop s« some ome degree of anisotropic structure when compared to the 


etartingcondition. 
Specimen men Preparation.- —The initial dimensions of all specimens were 8. + 
‘em high and 3. 38 cm in diameter. Two clays were used. Table 3 lists the 
Atterberg limits, specific gravity, and grain size analyses of both. i 
_ Kolinite. —The kaolinite specimens ° were prepared from “as received” ka 7 
_ linite powder. This was a particularly pure commercial kaolinite. The pow- 
- der was mixed with distilled water to give a water content between erage 
41%. A detailed description « of this clay and a study of its geologic origin» is. 
available. 9 These specimens were designated DWEPK. = = 
Boston Blue Clay.- —The Boston blue clay was obtained from a pit | in Cam- _ 
- pridge, Mass. In remolding, the clay had tobe allowed to partially dry in order 
to gain sufficient for extrusion into round bars. It was finally extraded 


rated BBC. 


TABLE 3.—CLAY PROPERTIES 


38 | 19 | 2810 | 98 


7 oo wool drains, together with filter paper strips along the sides ai and 
filter paper caps onthe top and bottom, were used intriaxial tests. The : speci-_ 
- mens were Sealed from the cell pressure media (water) with two latex mem- — 
branes (each 0. 056 mm m thick) separated by oe of grease. aeaad details — 
Non-hydrostatic Consolidation.—In order to reduce the void ratio variation | 
within atest series, the tests were run so that the only variable was the minor oe 
principal stress applied during consolidation. _ The major principal stress 
applied during consolidation was held nearly constant within each series, giv- 
ing essentially the same void ratio after consolidation in spite of the greatly | 
different (61/03) . values of 1.0, 1.2, 1.4, and 1.6 (Table 2). With e and 0; con- 
stant after consolidation, the effect of the non- hydrostatic consolidation could 
be better evaluated. 
The method used for consolidating the specimens in Series I, , I, Mand IV, 
_ which is referred to as standard in Table 1, was as follows: Using the aniso- on 


ee -%, “De-Aired, Extruded Soil Specimens for Research and for Evaluation of Test Pro- 
‘ceedures,” Matlock, Jr., C. W. Fenske, and R. F. A.S. T. M. Bulletin: No. 
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tropic loading device, the weights are ona hanger which transmits 
the load to the piston of the triaxial cell through a system of knife edge sup- 
ports with alever-arm ratio of 5 (Fig. 2). The load tothe piston is five times 
_ the load applied to the hanger. After noting the characteristic shape of the 
semi-log plot of the time-consolidation curves for kaolinite under a hydrostatic 
pressure of 3.65 kg per sq cm, it was observed that at approximately 6 ml de- Eq 
crease of volume, as meas ired by the amount of water drained from the speci- 
men into a burette, the consolidation curve starts the turn into secondary con- 

- solidation. . The total volume change was slightly more than 7 ml. It was de- 
cided to add the weights at equal intervals of volume change and to have all of — 
7 the peg asin added by the end ae ml or ip of the total consolidation. The 


weights available. This load increment was usually 1 kg on the specimen, and 
from four to eleven wey ery were applied during the approximately 25 min 
required to consolidate 6 ml. The cell pressure was applied in one step at the 
- beginning of consolidation aed remained constant throughout the piston loading 
cycle. A similar procedure was used forthe BBC in Series III andthe kaolin- 
ite at = 
7 i, All of the > anisotropic load was not added at the beginning of consolidation 
because the specimen was not able to support the total load without excessive a 
strain or possible failure. By adding the load after equal increments of vol- 7 
ume change, ‘the specimen was given time to consolidate and thereby increase 


strength to support the following load increment without obviously excessive 
Two other methods of consolidation were used in nee Vv. For t test H-27 
the specimen was first consolidated hydrostatically under a pressure of 2.28 
kg per sq cm. Afterthe condolidation had proceeded 1,370 min in the second- a 
ary stage, following an approximate 50 min primary, the axial load was then 
applied = increments until the maximum stress, o, was equal to 3.71 kg per 4 
sq cm. In test H-28 the axial load and cell pressure were both increased q 
six increments such that the stress ratio O75 /93 was equal to 1.6 at all times. 

The specimen was allowed to consolidate for 8 hr between each increment. The 
final stresses on the specimen | after consolidation were Oy = 3. 64 kg per sq 
cm and 03 = 2.28 kg | per sq cm. 

In all tests for Series I, I, and III the time allowed for secondary consoli-— 

oe after the final consolidation load increment varied between 17.4 hr and : 

Strength Testing. —As mentioned in the Introduction, the CFS- test was used 
for determining cohesion andfriction as a function of strain. ‘The compression 
rate of 0.0053 mm per min (0.0002 in. per min) was chosen the same for all 
tests. This rate permitted obtaining the required data for computing cohesion 
and friction over the range of approximately > 5%axial str strain com- 
of approximately 12 hr. 


non- ¢ consolidation ‘the specimen is supporting a a deviator 
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bal axial strain during non-hydrostatic consolidation 


axial strain due to = 
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FIG. 3.—DEFINITION OF SPECIMEN EN STRAINS 
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and has undergone axial strain as shown in Fig. 3. mm order to support this 
_deviator stress, the values of cohesion and friction zero. If 


1.82 ke/qm= 


‘CFS-test Axial Compressive Strain, in 


FIG. 8.- —VARIATION OF COHESION AND FRICTION WITH 1TH CFS-7 TEST EST STRAIN F¢ FOR 


on, 


Cohesi 


ee 


Friction 


otal Axial a1 Compressive Strein, Cons 9 


™ 9. 9. VARIATION OF COHESION AND FRICTION WITH TOTAL STRAIN FO 
specimen is to a hydrostatic state of stress, as in the case of hydro- 
- consolidation, then at zero CFS-test strain both the cohesion and friction 


were assumed to be zero. 
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 Forthe the non- pareinatic case, e, the ins method was used to compute the 
‘values of cohesion and friction at zero CFS- test strain: The values  ofcohesion 


int 


a" 


Cohesion, | 


(CFS-test Axial Compressive Strain, ing 
FIG. 10.—VARIATION OF COHESION AND FRICTION WITH CFS- TEST STRAIN | FOR 
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6? 


1.51 
“crs 
FIG. _VARIATION OF COHESION ‘AND FRICTION WITH T TOTAL STRAIN FOR 
SERIES IE 
“fully defined and Ce and are plottedagainst that strain. This may be seen 


Figs. eS 5, 6, and 7, that present the detailed results for Series II. Because it 
is impossible, from the data obtained, to draw the two curves for different ef- 


on 


Fricti 


friction are Computedonly at strains w wo required 0; Curves 
— 
aa 
— 
— 


fective stresses at small values of CFS-test strain (less than 0. 25%), an ex- 
‘oe procedure wasused. The extensive testing experience with DWE PK 7 
and BBC specimens hydrostatically consolidated for approximately 23 hr prior — 


‘OFS-test Axial Compressive Strain, +... in int 


12 VARIATION OF COHESION: AND | FRICTION WITH CFS- TEST STRAIN FOR 
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Cohesion, 


a a3 5 9 20 
‘Total Axial Compressive Strain, € 
 *FIG. 13.- _VARL ATION OF COHE SION AND FRIC TION \ WITH TOTAL STRAIN FOR 


o shear testing indicated that at CFS-test strains of less than 1% the frictional 
component of strength varies much less rapidly with strainthan the cohesion. 3 

Thus, the friction was extrapolated to zero strain instead of the cohesion. This” 
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S to compute the value of cohesion at zero strain by Eq. 3. ‘The initial 
as of cohesion and friction as computed by this method are shown on the 
_ curves plotted for the individual test (Figs. . 4 to 7) as wellas the curves show- | 
ing each complete series (Figs. 8 to 13). 
Series I, II, and IIT.—Series I, and lI were performed inthe same manner 
andare therefore presented together. , This analysis began with a review of the 3 
strain history ofeach specimen, as shown in Fig. 3. Afterthe initial cmreston 
_ the specimen was consolidated with (5/63) ratios up to 1.6. The applica- . 
tion of a deviator stress during consolidation resulted in axial compression 
- greater than the height reduction measured during hydrostatic consolidation to 
the same value of 01. This additional strain is designated €,. Then astrain 
<CFS_ was imposed during the CFS-test. Although the foregoing two strains 
were computed witha different reference height than the total strain €, , the per- 
centage difference in the heights ss , and ho)p is small and it is ‘considered 
adequate to estimate the total strain due to deviator stress val wy. 


The cohesion- friction- strain curves for these series are presented in. Figs. 
i to 13. _ Two strain coordinates are used. Figs. 8, 10, and 12 are based <la 
the CFS- test strain, €crg; Figs. 9, 11, and 13 on the estimated total deviator aa 
strain. Ce, . From both these sets of figures it may be seen that the cohesion be- 
" havior is a function of (01), and is unaffected by (33),. Thus, these results| 
the concept that cohesion is an isotropic property 
Nz Considering only Figs. 8, 10, and 12, it appears ‘that the friction- ‘strain be- 


ing values of de 
strain that occurre during the 1 non- -hydrostatic. consolidation. When 
on the basis of the estimated total strain, as in Figs. 9, 11, and 13, there is 
- some displacement toward lo lower friction values with increasing (51/93 er 
that the general friction- strain behavior is similar to that of the hydrostatically 
consolidated specimen. _ The values used for are critical in the foregoing 


apply ‘the deviator stress during ‘consolidation the estimate of. 
The estimate used herein, as shown in Fig. 3, is probably not sufficiently re-— 
 finedto permit a conclusion regarding any displacement effect of non-hydrostatic | 


7 > A conclusion is also not possible with regard to another important aspect 
friction behavior, pore pressure development. Pore pressure is controlled 
inthe CFS-test and thus its development in response to other factors possibly | 
associated with non-hydrostatic consolidation was not t considered in this re- - 

Series IV. —Three tests were run nto determ ine the effect of the length of _ 

time permitted for non-hydrostatic consolidation on the kaolinite. . The three 
tests showed no important d differences intheir results. . When the value of fric- 

- tion was extrapolated tozeroCFS-test strain to compute ‘cohesion at this strain, 
the Caen value of cohesion ranged from -0.08 kg per sq cm to 0.00 kg per 7 
sq cm. The © value (of friction was then mathematically adjusted to suit a co- 

The three Series IV tests were allowed to stand for 66 hr, 131 hr, and 299 
hr, respectively, after after the of consolidation. . From. the 
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COHESION 


sults ot the Series IV andprevious tests it was concluded that as the non- -hydro- 
statically loaded specimen continues to consolidate inthe secondary stagethe | 
: deviator stress is initially carried by both cohesion and friction but is eventu- 4 
ally carried totally by friction as the cohesion drops to zero. However, a zero 


— 


cohesion is only possible for — able todevelopthe required friction to carry 

_— The conclusion of the cohesion to friction transfer is strongly supported ww 

the unpublished experimental work of R. G. Bea,10 A. M. ASCE, and his im- 


on an 
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FIG. 14.—FRICTION ANGLE AT ZERO CFS-TEST STRAIN VERSU! 
CONSOLIDATION STRESS RATIO 


contribution to its is recognition is acknowledged. The present w ork on. 


 ‘Extrapolated fricti 


his subject is ; preliminary. Extensive additional research is in progres. 

_ The tests in Series I, I, and III were conducted after secondary consolida- 
tion times between 17. 4 hr and 19.7 hr. ‘This was insufficient time for the co-_ 
hesion to drop to zero, as seen in ‘Figs. 9, 11, and 13. The value of friction 
angle at zeroCFS-test strain, versus consolidation | stress is plottedin Fig. 
14 for test Series I, II, and I, together with the curve representing the value 


_ 10 “An Experimental Study of Cohesion and Friction During Creep. in! Saturated Clay,” > 
by R. G. Bea, thesis presented to Univ. of Florida, Gainesville, Fla., in June, 1960, in 
_ partial fulfillment of the requirements for the cctitiainnianss of Master of Science. a a 
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of friction that would be obtained had the | specimen consolidated over an ade- - - ‘ 
: quate length of time for the cohesion to drop to zero. wif cohesion is equal to 7 
then it from | Eq. 3 that 


The required time for the cohesion to drop to zero seemed to depend on 0; as 
: = as the type of clay. The BBC transferred cohesion to friction more ra- 7 


the same e clay at = 1.80 kg per sq cm. 
Series V.—Tests and H- 28 were consolidated | as described under 
__ heading “Testing Procedures” to determine the effect of the loading method — 
used to arrive at the final principal stresses applied during consolidation. 
There was no significant difference between the results computed from these 
ES and the comparable tests of Series Il which used the standard method. 
Considering the e differences inthe loading method used to consolidate the speci- 
mens it was concluded that the methodof consolidation had no affect on the co- 
hesion and friction of the | kaolinite tested. _ However, it must be noted that the 
time allowed for secondary consolidation after the application « of the last con- : 
Solidation 1c load about the same for all tests compared. 


7 

The followi ing is a. summary | : of the conclusions in regardto non-hydrostatic © 
consolidation and its affect on the cohesion a and friction for the _—" and ex 


consolidation did not affect subsequent cohesion-strain behavior. . This behavior 
also seems independent the consolidation procedure equal secondary 
- consolidationtimes are used. The results suggest that cohesion is anisotropic 
property in the two platey clay mineral soils tested. 
The value of friction at zero CFS-test strain, after a approximately equal 
PR Pho conditions, is increased by an increase of consolidation stress _ 
ratio, ( 31/83) ¢ - However, the variation of friction with strain is similar ical 
that of a hydrostatically c consolidated specimen. 
_ 3. During non-hydrostatic consolidation the deviator stres stress is carried at 
first by cohesion and friction but cohesion reduces and ee to zero” 


the soil has adequate friction capability. 


= . The anisotropic soil structure assumed to result | from n non- Si 


=n" 
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= cohesion at strain €; 


= void ratio; 


= specific gravity of soil solids; 
= initial extruded height of specimen (Fig. 8) 


tee height during CFS- test (Fig. 3); 
(oa = = = specimen height after non- hydrostatic consolidation; start of CFS- 


(Moh = = specimen height after hydrostatic consolidation (Fig. 3); 
7 excess hydrostatic pore water p pressure; 


= axial compressive strain due to CFS- test (Fig. 3); ¢ 


CF 


= micron, 0.001 mm; 


= effective normal stress; 


major principal stress (horizontal plane a 
effective n stress; 


3 during CFS-test (Fig. 1) 
Meter curve d ul ing (F ig. 
. iow = lower valued 0, curve during CFS- -test (Fig. 1 
minor principal st stress (vertical 


= effective. minor . principal : stress; 
effective minor principal st stress developed during consolidation; 


= = - - applied through vertical pis 1 piston; 
_ = hydrostatic stress in ‘triaxial cell; 


fo be symbols have been adopted for use in the paper. 
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shear stress on plane of of maximum stress for frictional 


angle of internal friction developed at start of CFS-test st = 


SM 4 = 
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SOIL MECHANICS AND FOUNDA’ TIONS DIVISION 


Proceedings of the American Society of Civil Engineers 


KY. 


‘Results of long term c tests on five. remoldedand two o undisturbed 
clays are reported. Most of the experiments were run until final settlements 
were attained. Basing on the present and past experimental evidences, the 
" settlement- time curves of secondary compression can be classified into three — 

7 types according to their characteristic behaviors. Briefly, for Type I curve, : 
_ the rate of secondary compression gradually decreases with time; in Type I 
curve, the rate is proportional tothe logarithm of time for a considerable range a 
of time, then rapidly decreases, , and for Type III curve, the rate increases with 
then vanishes. All experimental curves fall into one of these 


* was found that thistheory is adequate for clays behaving according to curves - 
Type Iand II. For Type Ill curve, however, significant deviations of the pre-_ 
dicted curve by this theory from the observed curve was noted. An extension — 
of the theory was made and c mparison of the theoretical and experimental — 
curves shows that the modified theory is satisfactory forthis special case. 
result, all the observed phenomena can be accounted. 
_ Observations on the entirely different and distinct behaviors between sam- 
ples inthe undisturbedand remolded state are described. Data for the relevant 
; > soil parameters, in particular the ratio of secondary to primary compressi- — 
bility and the rate factor, together with t their variations 1 water content 


Note.—Discussion open until 1962. To the closing date one month, 


of the copyrighted Journal of the Soil Mechanics and Foundations Division, ae J a 
of the American Society of Civil Engineers, Vol. 87, No. SM 4, August, 1961. — 
1 Civil Engr., Norwegian Geotechnical Inst., Blindern - Oslo, 
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the subjects v which much controversy in 1 the field of SoilMe- 
chanics is secondary consolidation. (Theterm secondary consolidation isopen _ 
to objections. As the compression proceeds under non-measureable 

¢ pore pressure, the process can hardly be designated as consolidation. Other 

y. terms like secondary time effects and secular effects have been used. It ap- 


~ in ‘the literature and will therefore be adopted indiscriminately in this paper 
as synonymous to secondary compression.) While for some clays oedometer 
tests have shown that secondary compression follows approximately a linear 
——: when plottedagainst the logarithm of time,2 deviations from Buis- 
‘man’s Law have been noted. In Mexico City Clay, for example, some experi- 
cL settlement-time curves exhibit a tendency for the compression to in- 
"crease with the logarithm | of time within the duration of the experiments. 3 For 
. two remoulded clays, a Residualanda Glacial Silty Clay, the settlement curves 
reported by G. A. Leonards, M. ASCE, and B. K. Ramiah, M. ASCE, show4 a 
distinct nonlinearity and approach horizontal asymptotes. The apparent con- 
-_ fusion arises mainly from the lack of a simple and rational theory on which — 
delineation of observations may be based, the extreme susceptibility of second- _ 
consolidation to environmental changes, (particularly temperature). 


The literature on secondary consolidation is extensive. review has 
i. been included in an earlier publication. 5 Of the theoretical treatments, the q 
works of T. K. Tan® and M. A. Biot? appear to be most significant. — 
; theories are, however, mathematically so complicated that the soilparameters _ 
involved cannot be conveniently determined, thereby both their experimental 
verification and practical application | are hindered. A theory based on the as- 
: sumption that the soil structure behaves as a rheological model has been de- 
_ Waleped. 5 The model usedwas originally presented by Tan and the theory may 
be shown tobe a special case of the generaltheory of Biot. The advantages and 
ed disadvantages of the present theory will be evident as the analysis proceeds. 


2 “Results of Long Duration Settlement Tests,” by A. S. K. Buisman, Proceedings, e 
First Internat]. Conf. on Soil Mechanics and Foundation Engrg., Vol. 1, 1936. 
7 3 “Consolidation of Mexico City Volcanic Clay,” by L. Zeevaert, Conf. on Soils for © 
Purposes, ASTM, STP, No. 232,1958. 

4 “Time Effects in the Consolidation of Clays,” by G. A. Leonards and B. K. Ramiah, 


$8 “A Theory of Consetintien of Soils Exhibiting Secondary Compression,” by R. E. 
Gibson and K. Y. Lo, Acta Polytechnica Scandinavica, 296/1961,Ci10. oe 
a 6 “Secondary Time Effects on Consolidation of Clays,” by T. K. Tan, Publ. Acad. 
* Theory of Deformation of a Porous Viscoelastic Solid,” byM. A. Biot, 
_ Journal of Applied Physics, Vol. 27, No. 5, 1956, p. 459. om (oe 
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SM 4 LAYS 
The results of long term tests on remolded clays and two undlstushed 
lays are reported herein. These data are analyzed inthe light | of the proposed | 
theory which was found applicable to most of the clays tested. An extension of 
the theory is also made to include the rather unusual behavior of some clays 
revealed inthe tests. . Consequently it was found possible to describe the time- 


= 


ia Notations. — —The letter symbols used in this ‘paper are defined where 
=§ first a appear, in the illustration or in the text, and ¢ aed aeranges ee 


“RECA 


an For the sake of completeness and dni of the method of analysis 


of test results, the basic ideas of the theory previously presented® will be ‘al 


7 


(b) MODEL OF MODIFIED. 


THEORY 


FIG. 1 _—SCHEMATIC. REPRESENTATION OF 


This: theory is stric’\. _oae- -dimensional, that is, motion of fluid and soil 
skeleton occur in one dir ‘ion only. The assumptions involved in Terzaghi’s | 
theory of consolidation are retained with the exception that the soil skeleton 
_ assumed to behave as a rheological model shown in Fig. 1(a). a 
When atime dependent effective stress o'(t) acts onthe model, the ll 
spring a compresses instantaneously, but the settlement of the Kelvin element 
; b-A is retarded owing tothe Newtonian dashpot ». Due to the low permeability 
of clays, the effective stress o' (t) increases gradually withtime from zero to 
a. the full value of the applied stress. Hence the compression of spring a is also 
gradual and is fully wheats applied stress has become 


effective. 
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_ Under the gradually increasing effective stress, ‘the Kelvin body commences 
to compress. Initially the full load is taken by the Newtonian dashpot A and is 
. —_ progressively transferred to the Hookean spring b as compression proceeds. . 
This phenomenon oftransference corresponds tothe process of secondary con- © 
-solidation occurring under sustained effective stress. After a long time has 
: é elapsed, the full effective ‘stress will be taken by a and b, the dashpot A sus- ; aah 
taining noload. It is also clear, therefore, that secondary compression is pre-— 7 


sent during the range of primary consolidation. 

a Through the model, an intrinsic time dependence of the effective stress- 

strain relation has been introduced. The consolidation phenomenon thus co! an 
4 sists of two time-dependent porcesses such as the dissipation of pore water — 
_ pressure and the creep of the soil skeleton under sustained effective stress. 


|The effective stress-strain can be in the form 


in which e is the compressive strain, a denotes the primary compressibility, 
refers tothe secondary compressibility, is the viscosity of the soil struc- 
ture, t becomes the time after application of load, and ¥ represents a a dummy 
“tan the continuity of flow of water through an element of soil is considered, — 


in n which z z is the : space e variable. - 
_ Eqs. 1 and 2, together with the appropriate initial and boundary conditions, 
define mathematically and one-dimensional problem completely. eS ee 
= The classical problem considered by K. Terzaghi, Hon. M. ASCE, isS a fin- 
ite layer of clay subjected to a surface loading of infinite extent. The pressure 
_ is instantaneously appliedand maintained constant thereafter. A complete so- 
7 lution through the use of integral transform method has been obtained for this 
considering ‘secondary ¢ consolidation. The initialand boundary condi- 
it for this problem are identical with those in the oedometer. This ager 


e is only concerned with the verification of the basic assumption made in the 


sults so that the following tic 


ers’ 


¥ 


in which Ao is the pressure increment _ _ symbols have been previously 
defined Eq. 3 is only valid afteran elapsed time ta, at which the applied stress 

has become fully effective. Its" usefulness lies | ‘in its sake for the purposes — 


4 _ 8 “Die Berechnung der Durchlassigkeitsziffer des Tones aus dem Verlauf der hydro- | 
‘dynamischen S Spannungserscheinungen, ” by K. Terzaghi, Akad, Wiss. Wien Abt. Ila, Vol. 
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large values oftimeis 
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4 stated. To process, 3s, the complete 
solution given in another has to be used. 


-Onthe basis of Eq. 3, two methods of determining the parameters a, b, and 
have in the analysis of the test re- 


tested |L.L.| P.L. ] 
(brown) 75 0 5} 52] 10 | 30-70 terme 
‘Clay | 40 - -2 .0 t per sq all 
‘Sodium 
Bentonite 
Fornebu Clay, 


“undisturbed | 44 | 21 sol 15 - 


Clay, 


Lilla Edet 


Clay, re- 


disturbed F 25 - 4, 0 kgpersqem 


quires t the value of final settlement _ to be obtained or r estimated. The second 
method uses a a plot of logio (52 - - 6 1) against ty, in which 59 and 64 are the 
- "settlements at times ty and tj, respectively, the interval At = tg - ty being kept | = 
_ Pay constant at any convenient value. This method can be used without resorting 
to final settlement, but involves the assumption that A/b being. constant during 7 
‘the working range, 

3 is purely a theoretical result, a consequence of basic 
of the validity y of the model. For ra clay to behave. - according to the theory, | 
therefore, it requires that Eq. 3 should be satisfied within its range of validity 
as indicated theoretically. | This, in turn a that each of t the | two plots i in 
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‘a question Uld pive 4 pine ranpe econ 
7 whichthe excess pore water pressure has practically dissipated. The theory 
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therefore conducive to direct ceca verification, and its validity and 


applicability easily assessed within its framework of assumptions. 


_ Longterm oedometer tests were conducted on two” British clays; r a 
London Clay and Grangemouth Clay, Commercial Sodium Bentonite and on two 
- Scandinavian Clays, Fornebu Clay and Lilla Edet Clay, both undisturbed and 
remolded. The index properties of the clays tested are shown in Table 1. 
- ati The geotechnical properties of London Clay have been investigated exten-— 
sively. 9 The geotechnical properties of Lilla Edet Clay have been considered 
by L. Bjerrum and T. H. Wu,10 M. ASCE, and the samples tested are nearly 
normally consolidated according to oedometer tests. The Fornebu Clay is 
_ normally consolidated and some relevant properties have been reported by | 


Bjerrum, N. Simons and I. Torblaa. _ 


_ The experiments on the first three clays shown in the table were conducted 


3/4: in. thick in a brass" ring. The last two clays 1 were tested at the 1e Norwegian 

_ Geotechnical Institute in a constant temperature room. The dimensions of the 

_ Sample were 20 cm2 in area and 2 cm high. ‘The rings w were lubricated v with | 

silicon greaseto reduce side friction, and filter papers were placed on top and 

bottom of the sample to prevent the clay from being forced into the porous stone. 

In all cases drainage was permitted at both ends of the samples. = as a 
The majority of the tests was performed | until ultimate a were 


reached and the temperatures recorded during the tests. 


FACTORS AFFECTING THE VALIDITY OF EXPERIMENTAL RESULTS 


Before proceeding to a description of the test. results, it is necessary 


4 


3 


_ briefly examine the factors which exert considerable influence on the accuracy — 


of the experimental data. 
“Side Friction.—The effects of side friction onthe consolidation behaviors of 


clays were investigated by D. W. Taylor who reported!2 that the total friction 


at the end of primary consolidation ‘measured was less than 15% of the total 


for wi Blue 


“same clay” in remolded state for a range of pressure ‘up to 8 ke per sq cm. 

Basing on these data, T. W. van Zelst, M. ASCE, estimated13 that the side fric- 
tion is approximately 10%for a ring of 1.5-in. height and 100 sq cm2 area. The 
7 results of their investigationled to the recommendation of the use of a ring of - 


“Tests on ‘London Clay from. Deep Borings at Paddington, Victoria and the South | 
Bank,” by A. -W. Skempton and D. J. Henkel, Proceedings, Internatl. Conf. on Soil Me- 
chanics and Foundation Engrg., London, Vol. 1, pp. 100-106. 
_ 10 “Fundamental Shear-Strength Properties of Lilla | Edet t Clay,” ihc L. ‘Bjerrum and 7 
Tr. H. Wu, Geotechnique, Vol. 10, No. 3, pp. 101-109. 
_ 11 “The Effect of Time on the Shear Strength of a Soft Marine ° Clay,” by L.Bjerrum, 
N. Simons, and I. Torblaa, Proceedings, Brussels Conf. on Earth Pressure Problems, _ 
(12 “Research on Consolidation of Clays,” by D. W. Taylor, Publication Dept. of Civ. 
and San. Engrg., Serial, 82,1942, 


13 “An Investigation of the Factors Affecting Laboratory Consolidation of Clay,” by 
. W. van Zelst, Proceedings, Internatl. Conf, on Soil Mechanics and Foundation Engrg. — 
Vol. 7, pp. 52- 
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3 suitable diameter to thickness ratio (perhaps 4) ‘in order to reduce the side - 
friction between the sample and the wall of the ring, 
Measurements of side friction have been further conducted byS. Hansbol4 
and G. A. Leonards, M. ASCE, and P.Girault,15A. M. ASCE. For steel rings 
of 4-7/16 in. diameter and 1-1/2 in. untreated with lubricants, Leonards and 
Girault demonstrated that the friction per inch height is of the order of 17% of 
applied loadfor undisturbed Mexico City ‘Clay. In all tests reportedherein, 
- brass rings were lubricated with silicon grease. It is likely that the side 
friction did not exceed 10% of the applied load in the majority of the tests. Be- 
cause in general, the side friction increases from the start of the experiment ; 


_ with the effective stress toa value which remains riley constant in 


“that the side friction does not vary to any significant extent in the range of sec-_ 
ondary compression, it will not interfere with the analyses which follow in the a 
However, it does alter the magnitude of the soil para- 
_ meters computed from the test data as the side friction reduces the magnitude , 
the effective stress acting onthe sample. 
a Vibration.—The effect of a sudden shock on n the process of secondary com- 
pression is clearly indicated inthe lowest curve on Fig. 3. (Different ordinate 
has been used for different settlement curves. The appropriate scale may be | 
ae d by referring tothe same symbol indicated in the curve and th the scale. ) 
_ The curve is displaced downwards abruptly, but retains the right shape. lt was 
found that intermittent vibrations would disrupt the continuity of the curve and 7 
render the results useless. _ Any activating force is therefore undesirable dur- 


_ Temperature.—The most important factor which affects the experimental 
results on secondary consolidation is temperature variation during the tests. r 
“This influence is clearly depicted in Fig. 2. . The. change of temperature 2 with : 
time is shown in the upper part of the diagram and the close response of sec- 
ondary compression to temperature fluctuations may be noted. An increase | 


” in 1 temperature of approximately 3°C completely alters the shape and trend of 


-—s-«‘It is evident therefore, that the importance of temperature control in all | 

tests on secondary consolidation cannot be overemphasized. 
susceptibility of secondary compression to temperature variation has 
an important bearing on the physical mechanism of f secondary consolidation, 


GENERAL ( CHARACTERISTICS OF SETTLEMENT- TIME CURVES 


-- number of tests performed on each clay, and the ranges of water con- 


tent and applied pressure c covered by the experiments are shown in Table 


, The presentation of all the results in graphical form would occupy too much | 
space, and therefore only typical curves for each bend have been included. The © 
7 4 settlement-time curves on a semi-logarithmic plot for the remolded = 


14 “Consolidation of Clays, with Special Reference to Influence of Vertical 
—, by S. Hansbo, Proceedings, Swedish Geotechnical Inst., No. 18, 1960. Nake 
15 “A Study of the One-Dimensional Consolidation Test,” by G. A. Leonards and nd P. 
Giraull, Internatl. Conf. on Mechanics and Foun Foundation Engrg., 
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are s shown in Figs. 3 to 7, and those for the undisturbed clays in Figs. .8 and 9. . 
For the settlement scale in Figs. 4, 5, 6, and 8, different ordinate has been — 
used for different settlement curves. The appropriate scale may be. obtained — 
by referring to the same symbol indicated in the curve and the scale. Ih Fig. 
8 the preconsolidation load pe is approximately 0.6 kg per sq cm from e- log 
. curve. In Fig.9 the secondary compression is too small to give any plot for 
= first two increments. For convenience, only the latter parts of > ortoacy nl 
_ dation curves are presented, needs eugtaning with the reading taken at the 
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FIG. 6 _—SETTLEMENT- TIME CURVES FOR FORNEBU 


CLAY (TEST SERIES I, SAMPLE NO. 
ss of the first a The s scale for settlement has been exaggerated for t the . 
sake of clarity of presentation, and the shapes of the curves are therefore much | 
distorted in the vertical direction. The numbers at the ordinate refer to the 
= occurring under individual increments eaineianees zero for each 
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moulded clays, irvenpective of the type of clay, ‘the settlement-time curves 


| during the secondary range have a similar shape throughout all the loading in- 
few exceptions, , the u ultimate settlement may be reached with- 


Ww, = % 


FIG. 1. _SETTLEMENT- TIME CURVE FOR REMOLDED LILLA 


in 3 weeks, the time ‘required increasing slightly with the loading increments. : 
The undisturbed samples, however, exhibit pronounced changes in the shape of 
the secondary compression curves for different load range For in- 
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‘FIG. 8,—SETTLEMENT-TIME CURVE FOR UNDISTURBED 
-FORNEBUCLAY 
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below the preconsolidation load, the ‘secondary time effect is is 
o. the final settlement is reached within a period of less than. a week. . At 
ot the co near the preconsolidation load, there is an abrupt change in the shape 
at the compression curve for the second and third increments in Figs. 8 and 
and secondary compression becomes. important. The shapes of the subse- 


quent curves are si similar to each other; some resemble for 


TYPES OF SECONDARY 
exar examination of the secondary c curves ina versus 


the a of time plot in all the tests conducted, together with a study of 


Tey 


\ II (b) 
—Type III (a) 


FIG. 10. PES OF COMPRESSION CURVES 


the available literature on the subject (Buisman,2 Koppejan, 16 Zeevaert, San and 
‘Leonards and Ramiah4) , indicates that there are three types of secondary com- 

_ pression curves which can be classified according to their characteristics as 

aa : shown in Fig. 10, in which curves for the whole consolidation process are il- 
lustrated. Curve Type I has a gentle curvature concaving upwards, the rate 

"secondary consolidation decreasing with time and the curve becomes horizontal 
-_ when the ultimate settlement is reached. The time required for this type ot 

curve to reach final settlement is comparatively short, being of the order of 
7 or so. Most of the remolded samples of different clays and 


16 “A Formula Combining the Terzaghi Load-Compression Relationship and the 
- Buisman Secular Time Effect,” by A. W. Koppejan, Proceedings, 2nd Internatl. Conf. : 


on Soil Mechanics and Foundation Engrg., Vol. 3, » Pp. 32-37, | 
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Fornebu Clay tested fall within this group. | — of Type I are characterized 
by the straight portion of the curve following approximately the logarithmic 7 
law for an appreciable range of time. The rate of compression decreases ra- 
pidly near the ultimate ‘stage and becomes zero as the final settlement is at 
_ tained. (It has been reported by R. Haefeli, M. ASCE, and W. Schaad,17 that 
z compression curve continued to be straight after a testing period of 5 yr. 
The clay used was Kaolinite. It appears, however, that the settlement must 
_ terminate finally because the sample cannot be compressed to nothingness in 
7 _ ‘the oedometer. In fact, the validity of the logarithmic law must be bounded, 
: - otherwise it would lead to the absurdity that the sample could vanish in a 
- laterally- confined condition. Letting the settlement S equal the initial height | 


2 
of sample h h in Buisman’ 


which is clearly finite. ) 
Type II curve is concave the rate of secondary compression 
creases with the logarithm of time, then slows down gradually and finally van- — 
ishes. The acceleration of compression canbe either gradualas in Type I(b) 
or abruptas in Type II(a). Bothtypes have been observed in undisturbed Lilla 
Edet Clay andexamples may be found in Figs. 9 and 18(b). The ultimate settle-_ 
ment for this type of compression takes a long period to attain and tests have - 
been running up toi a 100 days or more to reach equilibrium. 
Type ITand Type Il curves are those generally encountered, while Type i 
curve needs some special consideration, = 
environmental factors which may affect secondary compression have 
been presented. - During the tests on Lilla Edet Clay, no vibration from any 
_ accidental cause was noted. The temperatures varied but little throughout the 
duration of each increment, and this variance of 1 the order of + 1/2° C is thought _ 
tobe insufficient to cause a sudden break. A comparison of these curves with | 
Figs. 2 and 3 in which the effects of temperature and vibration are depicted 
gives a convincing evidence that this phenomenon must be a true behavior of — 
the clay. _ Test results on remoldedsamples of the same clay did not show this 
peculiar behavior and therefore further confirms this point | 
act It is that the peculiar behavior of secondary compression as 
ture. . The same phenomenon has not been observed in undisturbed Fornebu — 
Clay, which although not quick, has a sensitivity varying from 4to 8. On the © 
- other hand, the sensitivity of Lilla Edet Clay varies between 20 and 250. — 
all remoldedsamples of Lilla Edet Clay and the other claystested didnot show - 
this phenomenon, it must be a structural effect, which was developed during a 
_ time of geological scale after the deposition of the clay. This structural effect 
may thenbe destroyed in the course of remolding resulting inthe absence of _—_. 
oe phenomenon in remolded samples. It is quite likely that in the loosely de- _ 


17 “Time Effect in Connection with Consolidation Tests,” by R. Haefeli and W 
Schaad, Proceedings, aad nee. ome. on Soil Soil Mechanics and Foundation on Engrg., Vol. a 
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“posited — . clay y particles are are kept in 1 position by sc some form of cemen- 
tation which develops into a bond by chemical processes, aided by the effect of 
time and overburden pressure. This loose structure is susceptibleto a break- 
down when a critical stress (or more likely strain) is exceeded leading to an 
pet compression. (The question which inevitably arises is: Which 
5 observed quantity governs the phenomenon of structural breakdown, the stress 
or strain? Experimental evidence obtained so far does not allow a definite 
_ conclusion to be drawn, but it is hoped that additional tests will clarify this | 
point. Hereafter, therefore, we shall speak of a critical pressure loosely to 
mean a stress which has to be attained before this breaking of bonds can | 
occur.) The recognition of this phenomenon has been noted by Bjerrum and a 
‘Wu, 10 Hs reported the abrupt change in true cohesion (Hvorslev parameter) 
with water content at pressures near the preconsolidation load. The existence 
of the phenomenon of structural breakdown in other clays has been indicated 
sg Bjerrum18 as an implication from the low friction angle found in in-situ 
shear tests. It is suggested that a gradual structural breakdown ve lead to. 
a continuous curve concave downwards as shown in curve Type IlI(b) while a 


sudden br breakdown would lead to curve Type 


' = COMPARISON OF TE TEST | RESULTS | WITH THEORY — 


Figs. 111 to 17showthe analysis of test London Clay, Grangemouth 
Clay, Bentonite, Fornebu Clay and Lilla Edet Clay. (In Fig. 14 adifferent ver- _ 
tical scale other than that shown has been used for increment No. 3. In order — 
not to confuse the figure, this is notshown.) These plots give the relevant soil 
a parameters; the primary compressibility a, the secondary compressibility b, © 
and the viscosity of the soil structure 1~A. _As predicted by the theory, the 
(an closely approximate to straight lines at large times, while the devia- 
tions from linearity represents the hydrodynamic period of the consolidation : 


Fornebu Clay (both remolded and | and Lilla Edet Clay, 
- all behave in satisfactory agreement with this theory. For each clay at some 
typical applied load, a theoretical curve has been worked out and compared 
with the observed curve (Fig. 18). . (In Fig. 18(a) lines how observed curves _ 
and prints indicate corresponding theoretical computations.) The agreement — 
for allcases are again not unreasonable, especially in view of the e exaggerated _ 
__-Undisturbed Lilla Edet Clay requires special consideration. Type Ill curves 
"were observedin some loading increments for this clay, and the simple model > 
of Fig. 1(a) is inadequate to describe the phenomenon. In the following section | 


a an extension of this model will be made wh which will ich will permit analysis of the Type 


The model for the behavior of clays showing Type 
"curve characteristics is shown in Fig. 1(b). This model is a generalization of 
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‘the model of the theory previously presented. 5 To be satisfactory the model 


‘must be descriptive of the following experimental findings in undistrubed clays: 


(a) No appreciable secondary consolidation occurs at pressures below the 
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? (b) ‘When the applied pressure reaches a critical value, large amount of 
secondary compression occurs. The first stage of the secondary consolidation 
behaves in with and | Lo’ Lo’s theory. . Subsequently ‘signifi- 
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é CLAYS 


cant ot deviations ns from the the previous mentioned theory occur, the ‘settlement being 
underestimated by thistheory, 


= (ce) The sample behaves in subsequent loadings according to Gibson 


The model in Fig. 1(a) has been previously described. The spring elements — 

a and b behave according to Hooke’s Law, while dashpot A follows Newton’s Law 
of ideal fluids, that is, the rate of strain is proportional to the acting stress. 
The 1 spring a accounts for primary compression an and the he transference of s stress: 
from the dashpot to spring b constitutes the secondary effect. This model is 
-coupledin seriesto another similar Kelvin body through an element S as shown | 
in Fig. 1(b). The element S has the following mechanical properties: Below 
a certain critical value ofan observed quantity such as stress (or strain), the 
element S is rigid and therefore sustains the stress’ 7 '(t) fromthe b-A ele- 
ment, transmitting it to tothe by - Aj ele element. | However, 


1/A= 2.6 x 10° oom min Z *7.77x10°min’ 


|A/b- 3.8 x 10min" 


= 

= ‘FIG. 17. —DETERMINATION OF SOIL PARAMETERS OF UNDISTURBED LILLA 

_ EDET CLAY (TEST SERIES Il, SAMPLE NO, 2, 3RD INCREMENT; 

METHOD 1 WAS USED THROUGHOUT) 


cal ate is exceeded, the body S then loses its load-carrying capacity and 
oe (t) operates as well in the by - Aj element. The strain required to trans- 
formt the property of to bring into action the additional part by - is assumed 
to be an infinitesmal of higher order so that the element S does not contribute 
to the compression of the whole model. When S is rigid, the by - A, element 
is in effect uncoupled and the consolidation proceeds according to the earlier 
theory. For stresses below the preconsolidation load, the fact that secondary - 
consolidation is unimportant can be accounted for by letting the viscosity 1/)_ 
approach infinity, and in these circumstances only the primary compressibility 
a contributes to the settlement. On the other hand, when structural breakdown 
occurs, the whole model is fully active asa result of the transformation of the | 


the strese- -strain relationship described by this model, ‘the 


theory p previously can be extended with little and the entire 
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process ss of can be accounted for. For the of the 
"additional parameters by anc 1 Ay use ben made of the asymptotic sc solution 


Ad 
rA 


te 


curve occurs. For t, ta< t<t,, Eq. 3 holds. 
Using this procedure, the soil parameters for undisturbed Lilla Edet et Clay 
have been computed. As an example, Fig. 17 shows a typical plot. The plot 
on the left hand side gives the soil parameters before structural breakdown 
occurs and the one onthe right hand side yields the additional soil —— 
required to account forthe accelerated settlement a after breakdown at the time - 
roughly 4.5 by 104 min after application of loading. As ; predicted by Eq. 4, each 


in which te is the at which an abrupt in in the settl nal 


_ of these plots was expected to be a straight line within the appropriate limits. 3 
The soil parameters are of the same order of magnitude in both stages. Fig. — 
18(b) shows a comparison between the with (a) ) Terzaghi’: Ss 
theory, (b) Gibson and Lo’s theory and (c) the modified theory for undisturbed | 
Lilla Edet Clay (test series I, third increment, o' = 1.2 kg per sqcm. Points 
aa from Gibson and Lo’s theory and the modified theory are identical - 


OF SOIL PARAMETERS 
The’ characteristic behavior of soils during ‘consolidation can according to 
Gibson and Lo’s theory be described by two factors, the compressibility ratio — 
b/a and the rate factor A The ratio of b/a expresses the relative 


content for remolded London Clay, Grangemouth Clay and Fornebu Clay are 
shown in Figs. 19(b), 20(b), and 21(b), respectively. All three remolded clays — 
show a slight decrease in the ratio of secondary to primary compressibility _ 
of _ with increasing water content. _ Both for Grangemouth and Fornebu Clay the 
compressibility factor is of the order of 0,1 and for London Clay 0, 05. For 
Lilla Edet Clay the variation of this parameter is more erratic, but the ratio - 
la Ea order of 0,15. Throughout the ranges of water content and applied | 
pressures used inthe experiments, , it may be concluded that the fluctuations 
magnitude of this factor are not significant for remolded —_— #3  ~— 
_ The behavior of undisturbed samples are, however, distinctly different. For 
La Fornebu Clay Fig. 22(b) shows a rapidrise in the b/a ratiofor water contents — 
_ approaching that corresponding to the -preconsolidation load which is a result — 
s sampling for this normally consolidated clay. At pressures above the pre- — 
_ consolidation load, the compressibility factor can be ragarded as a constant © 
while for pressures below the preconsolidation load, secondary consolidation — 
iss so small that the secondary compressibility cannot be evaluated. A com- - ; 
parison of the compressibility factor inthe undisturbed and remolded state in- 
dicates clearly that secondary consolidation is much more pronounced in the 
undisturbed state; the compressibility factor neing twice or more andmay rise 
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loads. Thus in this clay the e secondary § settlement constitutes more than 60% 
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of the etttante compression. For pressures | ‘beyond th this, ‘the compressibility 
reduces to approximately 0 0,2 to 0, 3. 
The combination of parameters which governs the rate of secondary c consoli- 
dation is the factor A/b, as indicated in Eq. 3. The variation of /b with water 
content for London Clay, Grangemouth Clay and remolded and undisturbed 
_Fornebu Clay are shown in Fig. 19(a), 20(a), 21(a), and 22(a). It is both in- 
teresting and surprising (and indeed that there is little change inthe 
ratio for each remolded and clay. The value of this ratio for 
undisturbed Fornebu Clay is 1, 4 by 10-4 1 per min, while that forthe remolded 
samples is 1,9 by 10-41 per min indicating that the remoldedsamples willtake 
time to reach the equilibrium state. 
‘The peculiar behavior of Lilla Edet (Clay requires two rate factors to de- 
he the process of secondary compression at pressures slightly larger than 
7 the preconsolidation load. No general conclusion can be drawn from the data © 
at present available, nevertheless some idea of the magnitude of the rate factor 
Gy be obtained in ‘Fig. 17. Inthe remolded state, the rate factor of Lilla _— - 
Gly varies in the range 1,6 - 2,8 by 10-41 per min. ni 
A knowledge of the compressibility factor and rate factor will not only en- 
i“ us to have a general appreciation of the behavior of the clay during con 
‘solidation, but is also essential in the prediction of thetime rate of settlement _ 
infield cases. A brief analysis on aia of the theory has been given 
‘The variation of the viscosity of the soil structure of the clays 2 are oo ll 
in Figs. 19(c), 20(c),21(c) and (d), and 22(c) and (d). The tendency for 


a to . increase with increasing load and decreasing water content may 


CONCLUSIONS 


that secondary compression curves can be classified into three types — 


(a) I curve shows a gradual decrease of the rate of ‘secondary con con- 
7 solidation until it finally reaches ultimate settlement. = © 
(bd) Type I curve is characterized bythe proportionality of tian com- 


pression withthe logarithm of time for an appreciable range of time before the | 
finalsettlement is reached. 
tej T ype TI curve shows an acceleration of the rate of secondary consoli. - 
dation. This can be either gradual or abrupt. This phenomenon is believed to 
be due to the breakage of bond between the particles and occurs in loosely _ 
posited soils. The rate, however, decreases at at large times and zero 


With the environmental factors well controlled, it has. by 
experiments described that final settlement can be reached within practically - 4 
2 attainable time for all three types of secondary compression. Analysis of 
7 curves Type I and II led to the conclusion that the clays behave satisfactorily 
-— according toa theory of one-dimensional consolidation mn considering secondary 


‘Consolidation. 5 An oneneion of this theory is made, : and the behavior of soils 


to their characteristic shapes” ina of time 4 pe 
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-accordingt to Type III curve canbe an described b by the 1e modified theory. 
Therefore, it is possible to analyze all three ‘different types of curves of sec- 
ondary compression. Additional parameters required to describe the soil be- | 
havior can be directly and easily determined by the methods developed. It has _ 
been found that most of the tests on remolded clays showed a secondary com- 
pression curve of Type I. The undisturbed samples may fall into any one of 
the three categories and exhibit | marked variation in secondary compression - 
both with respect to magnitude and rate depending on the applied pressure and : 
the stress These phenomena can best be described by the two } ratios 


‘The factor b/a is the ratio of secondary to primary compressibility and is — 
a measure of the relative importance of compression due to pore water dissi- : 
: pation and that due to creep or breakdown of soil structure under sustained — 
effective stress. A slight decrease of b/a with increasing water content has — 
been found for remolded London, Grangemouth, and Fornebu Clay. 
‘Tests on undisturbed samples of Fornebu and Lilla Edet Clay showed, how- | 
_ ever, that secondary compression is much more important in the in-situ con-— 
dition. _ For pressures belowthe preconsolidation load, secondary compression 7 
is practically negligible. Within a small range of stresses exceeding the pre- 
consolidation pressure, secondary compression becomes important. For 
Fornebu Clay, it amounts to more than 30% of the primary compression and 
for Lilla Edet, secondary compression constitutes approximately 60% of the 
; ‘total settlement. At higher pressures, the b/a ratio becomes almost constant 


and is. ‘somewhat smaller than that obtained at pressures near the preconsoli- 


dation load. An impression of the 1e variation « of b/a with water content may be 


obtained from Eq. 22 for Fornebu Clay. 
‘Lilla Edet Clay in particular, a 2 rather ‘abrupt 1 increase in the rate of 


may be explained by the breaking of bonds between the soil particles (a struc- 
tural breakdown), | and this appears to confirm a hypothesis put forward by 
_ Bjerrum and Wul , who attributed the sudden change in true cohesion near the - 
preconsolidation pressure tothe destruction of bonding. 
ratio of governs the rate of secondary compression implicitly 
a parameter in an exponential function of time. For remolded London, Grange-_ 
-mouth, and Fornebu Clay as wellas undisturbed Fornebu Clay, this ratio is : 
essentially a constant for each clay withinthe range of water contents andpres- 
sures used in the tests. Data from other clays are insufficient to draw any © 
conclusion 1 regarding this point. This is indeeda fortunate situationas itavoids : 
the mathematical difficulty involved with nonlinearity. A comparison of undis- ; 
_ turbed and remolded Fornebu Clay shows, however, that the r/b ratiois smaller — 
in the undisturbed ‘state, indicating that the ultimate s requires longer 
time to complete in the undisturbed state 
‘The preceding c conclusions are valid within the frame work of one dimensional 
consolidation. applying the test results to practical cases, due consideration 
must be paid to field conditions . Special attention is drawn to the fact that, | 
> presented elsewhere, 5 the time rate of settlement involving secondary com-— 
pression is not ‘related by the ratio. of squares of thickness and for 
a the entire consolidation process, the complete solution has to be used. © Pa : 
_* _ ‘Finally, it must be emphasized that variation in ntemperature | has a predom- 
inant influence on boththe rate and magnitude ¢ of secondary consolidation. * For | 
investigations on secondary compression, therefore, this environmental factor | 
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ed. ‘Further of: the effect also probably 
lead toa better understanding of the mechanism of secondary consolidation. 
Reliable measurements of pore pressure during the entire process of compres- 
sion are — to give a still more stringent test of the aaa par 
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ee APPENDIX. {NOTATION | 


for reference: 
ae = primary compressibility, compressibility of of spring a; 


compressibility of spring b; 


(t) 2 strain at any time t; 
e (co) = ultimate strain; 


54 occurs; 
at which settlement ts occurs; 


~ ; 

upper time bound of Eq. 

upper time b bound of of E 4; 


= density of water; 


6 settlement at time tis 


settlement at time ty 
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= pressure increment for each loading step; 


= viscosity of f dashpot A aS 


= viscosity of dashpot 
variable of of integration. ion. 
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By Ardis White,! M. ASCE, James A. Cheney,? 2M. ASCE, 


- oA detailed experimental analysis is presented of the sturctural behavior of 
a segment of a steel sheet pile cellular bulkhead in Long Beach Harbor, Calif. 
‘The bulkhead consists of 62- * diam circular cells with connecting arcs and 

Measurements were 
‘made, during construction n and extending» more than 4 yr, of sheet 
pile strains and deflections, pore water pressures, and soil and bulkhead er 
The provide both empirical knowledge and a clearer understanding» 
ot cellular bulkhead performance. Of particular interest were the hoop ten- — 

_ sions and associated lateral earth pressures, the consolidation phenomena, 

_ _ the distortions of sheet piles and soil. Useful information was obtained about - 
field instrumentation and the importance of certain construction procedures, — 


ae) ‘The t test results are provided in as full detail as feasible. ors 


=F 


__Note,—Discussion open until January 1, 1962, To extend the aontans date one — 
= a written request must be filed with the Executive Secretary, ASCE, This paper is part | 7 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American 
of Civil Engineers, Vol. 87, No.SM 4, August, 1961, 
Assoc, Prof, of Civ, Engrg., Univ. of Houston, Houston, Texas, 
Research Engr., Lockheed Aircraft Corp., San Jose, Calif. | . 
3 Prof. of Engrg., Univ. of California, Los Angeles, ‘Calif, | 
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"INTRODUCTION 


_— sheet p pile cellular bulkhead is still young ung phenomenon, while 
design methods have been developed4,° there have been to date (1961) only few 
attempts to obtain quantitative data on the performance of such 
the field.6 Such field studies have been undertaken in Japan. 
_ The behavior of soil itself as a structural material is not adequately un ll 
stood, and the complex interactions between soil and steel in cellular bulk- 7 
add greatly to the difficulty of attempting a rational design of 
— Laboratory and theoretical research can clarify some of the ques- 
= Field observations, however, are necessary in order to gain empirical 


complex interactions which both sampling type field 

studies, involving simple measurements of the gross performance parameters o 
on many bulkheads of differing design, and detailed field studies of particular 
cells in particular bulkheads, The latter study is presented herein. 
Construction in 1952-1953 of the cellular bulkhead at Pier E in Long Beach — 


opportunity for this | ‘study, ‘anda cooperative field ‘research. program es- 
tablished between the University of California and the Port of Long Beach. A 
_ pilot test was conducted at one cell. Another cell and adjacent arc of the bulk - 
_ head were fully instrumented for measurement of strains in the ‘steel, pore 
__ pressures and settlements of the fill, and movements of the cell and arc. ~Ob- 
servation was made of gross horizontal movement « of all cells. The field mea- - 
“surements ‘covered the period from July 1952 to September 1956. Previous re- 
| appeared in the form of theses by White and Cheney.2,19 
Pier E consists of 33 circular cells with connecting arcs on water side, as" 
illustrated in Fig. i - Original water - depth varied from 30 ft to 50 ft at differ- 
ent locations along the bulkhead, as shown by the contour lines on the figure. : 
_ The normal tide fluctuation has two highs and two lows per day, with a maxi- 
annual range of approximately 8 ft. All elevations herein are referred 
™ mean lower low water (MLLW). 
Qualitative Description of Bulkhead Performance.—It will be useful herein 
to give a qualitative description of the structural behavior of a circular cell 
= founded in sand, When the steel work has been completed, the cell may be | 
Gauge of as acylindrical assemblage of loosely attached sheet piles with ‘their 


ey 4 «Stability and Stiffness of Cellular Cofferdams, ” by Karl ate Transactions, os 
ASCE, 
5 “Soil Mechanics, ” by D. P. Krynine, McGraw- -Hill Bock Co., New York, 1947, om 
with Cellular Walls on a Marl Foundation, by Miller and 
7 Ishii, Transportation Technical Research inst., 
‘hama, Japan (personal communication). 
8 «Steel Sheetpile Cells Form Bulkhead for Long Beach Harbor Expansion,” by R.R. a 
Civil Engineering, September, 1950, 
lisee 9 “Field Study of a Cellular Sheet Pile Bulkhead” by Ardis White, thesis presented to 7 
‘the Univ. of California, at Los Angeles, in June, 1953, as partial fulfilment of the re- 
quire ments for the degree of Doctor of Philosophy. 
_ 10 Soil Behavior in a Cellular Bulkhead,” by James A, Cheney, thesis presented 
the University of California, at Los Angeles, in June, , 1953, as partial | fulfilment of t the 
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CELLULAR BULKHEAD 


lower en ends. embedded in in sand. As pu pumping ng of hydraulic fill into the cell pro- 


“4 Lateral pr pressure of on fill deposited, and any on tetenatiite excess due to 
water held in the cell above sea level, cause the cell to expand inelastically — 


—25‘PILOT TEST CELL 


1.=PLAN 0 OF PIER E CELLU LAR BULKHEAD - 

=. 2. Weight of the fill cau: causes an an increase. in n outward pressure of the sand be- 
low the original bottom and the cell tends to expand a at t that: level also v until the 
interlock slack is filled over the full cell height. 
3. As soon as the interlocks are yh the cell ceases to be a looseas assem- a 


— 
— 
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cally there will then be an mene of oo tension over most of the full | 
height for each increment of fill placed or hydrostatic excess spplied. While 
_ the hoop tension at any time will increase downward to a level in in the vicinity 
of the natural ground, this effect probably will not be reflected in gross changes © 
in shape of the cell because the corresponding eetormations will then be large- 


ly elastic and therefore small, 


4 In normal construction procedure, the cell becomes filled with water 
soon after pumping starts and remains full until the fill level reaches the top 
7 of the cell, _ Subsequently the water level falls slowly to sea level. _ Each in- 
‘crement. of fill consolidates under its. own weight and of the fill placed 
above it, so that settlement of the lower fill is advanced by the time the last 
- fill is placed. Essentially full settlement of a sand fill will have been accom- | 


plished before any ba backfill is placed behind the cell. ~ a) 


The cell may now. be considered as sand- -steel structure 


strength to resist horizontal loading. When n fill is placed and on the 


le -4 ‘Hoop tensions on the back and sides will be diminished slightly due to. 
_ the inward pressures exerted by the backfill. CS — > a 
cell as a whole will rotate and translate bys amounts on its 


he The c cell will deform, with important shear stresses and strains set up 

‘- horizontal planes } and on vertical planes parallel to the bulkhead centerline. 

Some sliding will occur at sheet pile interlocks in the side walls, a 
& 4. Asurcharge placed behind the cell will aggravate the previous s tendencies. 

‘Soil surcharge placed directly over the cell will aid it in resisting rotation, 


‘sliding, and distortion, a will also increase the hoop tensions re 


a some of the p previous ; relationships for the edeiding bulkhead tested. 


STRUCTURE AND MATERIALS | 


The construction procedure used at Pier E is illustrated approximately in 7 

Fig. 2, whichis a . photograph of Pier A in the same harbor. At the lower right 

_ of the figure is the king-pile support for the steel template, in the next posi- 

_ tion the template is in place, . and succeeding positions show the subsequent 

_ phases of construction. ‘The ‘procedure at Pier E was different from that il- 
- lustrated in Fig. 2 in that arcs were placed on the water side only ‘and were 


placed before filling of adjacent cells. 


__‘The sheet piles were typed SP- Ta, 1/2- in. web, 16 ft long, and type SP- 6a, 
ta 3/8-in. web, 69 ft long, manufactured | by t the Bethlehem Steel Company. _ The 


properties: given the manufacturer v were: — 


Modulus of elasticity 29,500,000 
Ultimate tensile ‘strength ®, 000 psi ‘minimum 

‘Fig. ‘method of fabrication | of the T- pile between cell and arc. 
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FIG, 3 .—GENERAL STRUCTURAL FEATURES AND CELL 16 FILL 
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‘The test cells x were ere built essentially in accordance with Fig. 3, but due to 


‘shape fill was Details of this ir will be given with 


ar Soil Characteristics. —The chronology of filling of the main test cell, to- 
gether with a summary the soil stratification and simplified soil 
_ (3-in. diam specimens) are given in Fig. 5. The dates on the fill profiles _ 
- are for 1952. The soil below el -80, down to approximately el -110, wasofes- 
"sentially: same constitution as that at the higher le levels. The be- 


TABLE 1,—FILL STAGES, CELL 16 


Cell 4/10 full 
9/27 1/52 Cell Lfull 
Ares ful 


Backfill up 
on 


L 
—cROSS ‘SEC TION | OF TYPIC AL T-PILE 


; hind the c cells was obtained by hydraulic dredging from the harbor south of the 
; _ bulkhead. . The fill stages shown in Fig. 3 and used as an independent variable 
| presentation of test results are defined in Tablel, | 
; Contrary to experience with the other cells, the water level in Cell 16 did 
not rise to the top and thus stabilize the cell during the early period of pump- 
Z ing, but remained essentially at sea level until September 19, 1952 (fill pro- 
file, Fig. 5). This is attributed to a rather free flow wot soil and water under the 


— 
| 
— 
— 
its _the jetting used in driving the test piles. On September 19, the ground level a 
dae 
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August, (1961 
outside: the northeast corner had risen several feet (corresponding w with an an ef- 
fective penetration depth of approximately 30 ft) and the flow ceased, the water 
level came up in the cell, and the cell stabilized. It was normal procedure to 
_ pump into the cells in 1 such a way that the fill piled up near the south face. At 
Cell 16, the unsymmetrical soil pressures accompanying this procedure fur-— 
ther distorted the poorly alined cell before the became 


‘TEST ‘PROCEDURES 


‘The principal observational results a obtained at Cell 16, but supplemen- 
tary information was gained at Cell 25 and with the bulkhead as a whole. ii 
At Cell 25a pilot test was. conducted, with horizontal and vertical strain- - 
meters, a driftmeter tube, and a pore-pressure settlement device. Workat _ 
Cell 25 covered the period June 1952 to March 1953, during which the cell was 
and filled and backfilled to el +4, The ‘strainmeters were at Level 
* 6) which, in the case of Cell 25, was 7 ft below natural ground. _ _ 
_ Observations were ‘made of the lateral movements of all cells. on 


- entation. The driftmeter tubes were used to obtain sheet pile deflections. The 
pore-pressure | settlement devices were located at the four levels at points O 

U. In addition to the measurements with the indicated instruments, 
observations of vertical and lateral movements were made with surveying in-- 

struments at the top of the cell. 

Installation and observation at Cell 16 covered the period August 1952 to 
September 1956, during which the cell was ; driven, filled, and backfilled,a sur-_ 
7 charge was applied, and, finally, work was done to realine the cell. Complete 
sets of observations were taken at approximately 10-ft intervals of fill in the 
cell and every 2 days or 3 days for approximately 3 weeks after completion of 
‘filling « of the cell. On two occasions readings were taken at extremes of tide 


over a 24-hr period, Partial sets of were taken at decreasing fre- 


_ Sheet pile strains and deflections were attained on adjacent piles, one pile | : 
strainmeters andthe other holding the square tubing for the driftmeter. 


_ mounted in pairs, one inside and one outside the cell at each location and ori- 7 


_ Straight piles were selected for the test piles and the immediately mee 
‘piles t to minimize irregularities i in load transfer across t the test piles. 

Regional Subsidence. —During the test period, general was 
_ progressing in the harbor area. The rate at Cell 16 was approximately 0.9 ft 
per yr and that at the reference bench mark on shore was approximately 1.1 to 
yr. The rotational effects of this “subsidence > were. determined to have an 
_ unimportant influence on the test, and the citorential settlements over the 

-_ period of filling of the test cells were not significant. | However, corrections 

_ forthis effect were required for determining cell settlements sein time | 


- obtain fill surface configurations - (Figs. 3 vind 5 at the fill stages and at other 


times when instruments were read. ea Water levels inside and outside the cell 
_were also recorded. Elevations were referred to a bench mark on a sheet 


7 pile at Cell 18, with periodic checks against the nearby shore station. ‘Spawn 


Settlement and Pore Pressure.—The divide shown in Fig 7 was used 
| ‘to determine both the | pore water T pressure and the settlement of fill at the the ‘sev 
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2 2FEET 


4 
HOLES DRILLED 
f 


__KNIT COTTON SACK 


- FILLED WITH PEA GRAVEL 


Proves 


DEVICE 


= 

| 


at the a of with an engineer ing level, utilizing: the ‘previ ious beach 


“device 
same devices were installed : as the fill reached the desired levels and were 
“held in position temporarily by fastening them to timber piles laid across the 
top of the | cell, When the fill had risen around the device sufficiently to sup- — 
= its weight, the hegre was changed to provide horizontal restraint only. ; 


‘cana movements of the fill accompanying the early construction procedure, _ 
—_ _ Cell Deformation. —Relative horizontal deflections of points at 6-ft inter- 
vals along the length of seven n piles were c obtained by means of a a . driftmeter - 
_ This instrument is an oil well surveying device which | gives a photographic 
- record of the slope at any desired point within a pipe. In use, a photo-sensitive 
disc was placed in a . special seat in the base of f the driftmeter (bottom of Fig. 
8), and the three small parts v were assembled and screwed as a unit onto the 
= base. A timing device was set and the three-part unit was placed - the hollow 
weer part of the instrument, to which the base was then screwed, ‘The drift - - 
- meter was horizontally oriented, dropped down inside a square tube, and held 
si at the proper depth with the steel tape. After an interval, the timing device _ 
activated a flashlight, which directed light down a hollow plumb bob onto van 
_ sensitive disc. The light burned for 1 min, causing a dark spot on the disc, 


the position of which iothented the two components of inclination of the drift- 
meter to0.1°., | 
Mounting of the square tubes on the piles consisted of making up the 
Pp lengths of 1¢-in. (outside dimension) sq galvanized iron tubing, using wa- 
terproof soldered joints, and fastening this tubing with welded clips to the pile. _ ; 
_ ‘The tube was kept free of direct bond to the pile. The lower end of the tube 
- was sealed and then protected by a short « oblique section of $. in. ‘pipe welded 
to the pile. The piles with driftmeter tubes could be handled safely with a 3- - 
r point pickup, and id were assembled in position in the cell and then c carried iad 
a Marks were established at the tops of all instrumented piles and also at a 
pile in the northeast quadrant between D and (Fig. 6). Horizontal move- 
v ments of these marks were observed with surveying instruments to 0.03 ft, 
and vertical movements to 0.01 ft. Control of the horizontal movements was a 
4 point at the > center of Cell 16, which was reestablished periodically by inter- oo 
7 section of a line from shore and the: theoretical centerline of the bulkhead, 
Control for the vertical movements was the previous bench mark. 7 
The lateral movements of all cells were determined by surveys made after — 
completion of backfilling behind the cell (March 18, 1953) and after applying © 
"surcharge (June 23,1954), 
Strains in Sheet Piles. —Carlson strainmeters were selected for the Pier E 


_ research work because of success with them in a previous research project in 


~~ Long Beach Harbor. il — were mounted in J pairs on opposite cnee of a 


$$$ 
~ “Field Study of a Sheet Pile Bulkhead, ” by Cc. Martin Duke, Transactions, ASCE, 
Vol, 118, 1953, pp. 1131-1196. 
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_ FIG, 8.—DRIFTME 


> FIG, 9.—TYPICAL STRAINMETER INSTALLATION 


7 a U-shaped meter lugs (Fig. 10) \ were welded in saeco using a iia to establish 

"accurately the proper gage length of 10 in. _ After positioning the conduit and 
welding it in place, , strainmeters were set in the lugs, adjusted, and the set 
screws tightened. _ The three- -conductor neoprene- -covered lead wires were 
then pulled up through the conduit for | collection at the top ofa pile and were 
labled. Ends of the strainmeters were coated with anti-rust compound, and 
individual sheet metal channel covers tog 10) were | tack- welded at one e end 


—— TRA NM TER 

— cone POINT, ‘ALLEN 

HEAD SCREW WITH 


X28 THREADS 


} WELDED ALL 


AROUN 
¥ [CENTER ONLY 


| SYMMETRICAL 
ABOUT ¢ OF 
WEB 


FIG. 10,—DETAIL OF ST ETER INSTALLATION 
To prevent damage to the strainmeters 3 while moving the piles, a timber 
- pile was lashed to each strainmeter test pile and was kept in - until the — 
test pile was in a vertical position ready to be placed in the cell, The timber 
7 pile was then removed, and the test pile was interlocked with the previously 
_ driven adjacent pile and was jetted to grade. The lead wires were grouped at 
_the piles and were later extended to a central point. © 
Only one strainmeter or its lead wire, at Pile E, Level 3, became defective 
up to the time of completion of backfilling (August 22, 1953), but subsequently | 
several lead wires were lost or destroyed due to construction operations ac- 
companying surcharging and the later realinement of the cell, Of the meters” 
whose lead wires were still intact in September 1956, ey 5% were 
7 _ Zero strain readings used in the computations we were those » obtained on Sep- 
ra tember 8, 1952, when ‘fillin be ae at that tim time the interlocks were 


nding strains. Fig. 9 illustrates a typical in- 
7 4 a box over the entire space between conduit ends (Fig. 9). The third meter in ' 
Pe hi the rosettes at Pile C was mounted at a 45° angle between horizontal and ver - in 
— - GAGE G.I. 
! 
| 
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loose, readings at Levels 1, 2, and 3 could be made only-to the nearest 40 by 
10-° 6 of strain. Zero readings at Level 4, - subsequent readings at all a towel 


nearest 8 by 10-6 
DATA PRESENTATION 
soit and Water Behavior. —Tables 2 and 3 present the data on soil settlement. 7 
and pore pressure in Cell 16 for the various instrument locations. _ Also given - 
- d in Table 3 are the elevations of the tide and of the tree water level in the cell. 
are illustrated in Figs. 11 li, 12, C 


TABLE 2,—SUMMARY oF SOIL, SETTLEMENTS, CELL 168 


9-9 | 9-18 | 9-22 | 9-26 | 9-29] 10-4 2 | 2 
1952 1952 1952 1952_ 1952 1952 1954 


3 
7 
1 


2 
95 
4 
4 
2 
3 


1- 
19 
0. 
0, 
0. 


nw 


Figs. 11 and 12 show the soil settlements at the four levels during an and > 
shortly after the filling of the cell. Fig. 11 illustrates the geometrical dis - 
tribution of final settlement onSeptember 22 and October 15, 1952, and Fig. 


_ shows the average ae rates at Levels 1, 2, and 3 in relation to depth of fill. 


and water flowed under the cell the early filling period. 
‘Figs. 13 and 14 are typical of the pore pressure distribution at the time 

_ when filling in Cell 16 was just completed (Stage III). Fig. 13 gives the ex- 7 
perimentally determined equipotentials at at one instant, and Pig. 14 portrays — . 
the manner in which pore pressures at three levels responded to transient in- | 


side and water levels. 
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FIG, 12,—TIME-SETTLEMENT CURVE 
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Deformation.—The data on deformation ar and “movement o 162 are. 
given in Figs. 
Fig. 15 shows the two components of horizontal deflection of each instru- 
mented | pile at each fill stage. The bottom of the pile was taken as reference, he an 


_ and the curves were constructed by successive upward projections of the an- 


cell were observed a few times during filling. The deformations were sion 
to those for Pile E, Cell 16, for the corresponding loading. 
Fig. ‘16 ripen the changes in elevation and on of the top of the cell. : 


™ 


cell c on n the plane x X - 


. 
4 
ig BIG, 18,-AVERAGE RADI = 
x, which is the plane of resultant horizontal move- 
P 


‘Fig, 17 shows the changes in shape and position of the cell for the various» 
stages, based on surveying data. Fig. 18 combines the radius measurements 
at the top with the driftmeter readings to obtain an average deformation of the 


In Table 4 are shown the necieeaes deflections of the tops of Cells 11. 
ts 31 (Fig. 1), based on surveying data a corresponding with application of . 
surcharge. 


, _ Stresses in Sheet Piles. —In Table 5 are e given the horizontal hoop tensions 


were also computed but are not included herein because they were ‘e quite small | 
: and were more or less randomly distributed about a zero mean, —_ DAPI 
_ ‘The method of computation of the stresses will be explained. The strain at 

"each location was first computed, for each date, as the average of the strains 


Southward Deflection, in Feet 
8/18/53 to 6/23/54 


0 
0 

0 

2. 

0 
1, 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0, 


_ just before filling was started, as the zero —_—- value. The strains were then 
converted to stress, using the steel properties previously given, and utilizing 


stress theory to obtain principal stresses, and final- 


| 


In the application of biaxial stress theory, three procedures were ere necessary — 
1 At Pile C, since there v were strain rosettes, it was s possible to obtain the 
arse normal stresses uniquely from the strain measurements.  ~ 
cea 2. At Piles A y 3 _E, and the the pilot test, there were horizontal and vertical 


strainmeters at each location, The ‘normal stresses were obtained with the 
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assumption ‘that | the measured strains were the principal s strains, ‘This 
= is rational for these four piles, because, due to symmetry, the inter- 
lock friction theoretically is zero for each of them at all times. _ 
3. At Pile D, only horizontal meters were used, and the horizontal stresses _ 
o were obtained by assuming that the vertical stresses were zero and increasing © 
the computed horizontal stresses | by 10%. This procedure w was based on a de- 
- tailed study of strains and stresses from the two- and three- element rosette 


TABLE 6 .—AVER AGE HOOP TENSIONS, 


ore Co 
ANON QDI 


In Figs. 19, and 2 21, the average” tensions in Cell 16 are sho 
graphically i in their relation to some of the other variables. These figures were 
_ plotted from the data in Table 6, which was prepared by averaging stage data 
_ from Table 5, — Fig. 19° shows how the tensions varied around the cell at each 
: _ fill stage for each meter level. Fig. 20 gives the chronological picture of de- 
_——- of tensions at the four levels, averaging the data on all piles at each 
level. The same information in Fig. 20 is plotted in Fig. 21 against level for 
each fill state, 
The few ten tension saan from th the 7s test are enon as follows: 


Hoop yop Tension Fill Conditio 
(kips per lin in.) 


52 Cell empty 
1/26/ 52 Fill at el 0 
7/30/ : Cell full 


~4 
— 
| 
mm 
| 
iva 
— 


CELLULAR BULKHEAD 


100 


10-30-52 


-27-52 


StageI 

9-18 52 9 


one 
4 


Nn 


NOISN31 dOOH 


‘ 


qv 
|; | | &¢§ 
— 
4 
| 
— 


Hoop Tension 
“(kips per lin in.) 
— 


‘Arcs full 


3/1/53 ‘1. Backfill at el 4” 

| PERFORMANCE DURING FILLING OF CELL 


Settlement. —Soil settlement entered the of cell t in three 
ways, ‘Firs st, the soil within the cell consolidated under its own weight. . Sec- 
ond, the weight of the filled cell caused consolidation of the original amen, 
‘Third, the lateral pressure of the fill depostied unsymmetrically in the cell 
(Fig. 3), and later of the backfill, produced a differential settlement resulting 
7 in a southward tipping. Superimposed on these effects was the regional sub-. 


sidence in the Long Beach-Wilmington area, 
ee. maximum settlements of the pore pressure- -settlement devices were 
fou 


und at depth rather than at the surface. Points some distance down in the - 

cell settled as much as 1. 2 ft (Fig. 11), — the estimated effects of the | 
_ This is | because the fill far - 

= down in n the cell was in n place longer, was more heavily loaded , and = 
-perienced some settlement before the top material was placed. 

The maximum soil settlement shown in the north-south | porfile (Fig. 11) 
occurred in the center of the cell, with smaller settlements at the edges. This _ 
could be caused by the effects of wall friction. The apparent inconsistency at 
_ the lower levels of the east- -west profile is because the early soil flow beneath 


east face of the cell near point U. 
v7 Some sliding and shifting of the fill within the cell occurred while the pore 
settlement devices were still in shallow fill, This caused erratic 
movements of some of the lower settlement points and made difficult the ac 
curate determination of consolidation settlements. 
Pressures.- —The changes in water level in the pore pressure pipes 
were due to three pressure . effects such as hydrostatic excess due to consoli-— 
dation, rising and falling tides, and seepage of water out of the cell. The con- 
_ solidation pressure can be determined | analytically from the measurements a at 
a given time if the sea level and the flow net through the cell are known. | we 
_ Unfortunately, data from the borings as well as the features of the profile 
plots of pore pressure data (Figs. 12 and 13) indicated that a theoretical pre- 


from: the measurements if the “assumption - is made ‘that the consolidation hy - 
drostatic excess does not change between ‘successive extremes of tide. or 
level of water in the standpipe is ages en by 


_in which, C is the water level in the cell; S denotes the sea level; h represents | 


the excess hydrostatic pressure due to consolidation; and 


a number of potential drops from fill eines to a given p¢ oint. 
total number of potential drops 
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1: and 2 Sent high and low tide, respectively 


~ 
es defines the flow net empirically. a the flow net is determined, h ma may 
= 
_ This computation was carried out for point “S” (Fig. 6) for September 27, 
1952 (Fig. 14), and it was determined that h =1, 2 ft. However, because the re- 
sult is obtained from the difference between relatively large numbers, a small 
aa error in pore pressure reading leads to a larger probable error in| 7 
the computed hydrostatic excess. Therefore, the previous value, while 
‘Early in the filling, when the sea level and the water level in the cell were 
< he same, the existence of hydrostatic excess due to consolidation was obvious. 
Between September 10 and September 12, h ranged from 0 ft to 1.2ft. we 
_ Movement of Top of Cell 16.—Fig. 17 ‘shows that the top of Cell 16 was far 
out of its theoretical position before any fill was placed. — _ Except on the east — 
= it was initially too far east and south, and the cell continued to move in 
_ these directions during filling. This movemeat probably was —" due aie 
"uneven deposition of fill and to the early flow of soil beneath the cell. 


‘re by one or more of the following: (ia damaged template was - used in 2 
_ forming this cell, (2) jetting was done to place the test piles, (3) relatively 
‘shallow penetration king piles were used to support the template, and (4) re- 
_ latively shallow penetration of sheet piles resulted from the low natural ground ~ 
level, Great difficulty \ was in closing this cell, the most irregular 
Fig. 16 shows that the maximum sheet pile settlement, at point A, on com- 

: pletion of filling of the cell was approximately 0.2 ft. ‘The greater settlements 
on the south side of the cell are consistent with the data on southward d tilting» 
"provided in Figs. 15and 17,0 

During filling, the driftmeter data showed rotations at the lower ends of the 
square tubes (3 ft above the end of pile) of eee © 0. i for all test il 
except at D, in which the value was 2.0°. 
 ~— Radius Versus Depth Curves. —Correlation of radius measurements at ate 

top: of the cell (Fig. 1%) with the driftmeter data (Fig. 15) permitted the con- 
struction of Fig. 18, representing average radial sheet pile deformation. The 
: general shape of the curves clearly indicates outward movement of the bottom — 


of the cell, increasing with height of fill. The initially large bottom radius 
(Stage I) suggests: the possibility that an interlock may have been injured dur - 


ing pile driving or perhaps that an insufficient number r of points was used to - 
define accurately the shape of the cell. The lateral distance between the curves © 

7 for Stage I and Stage IV represent both an elastic and an inelastic _ (interlock 

takeup) deformation of the cell wall. ‘The elastic portion of this 
is so small, theoretically, that the deflection measurements could oe detect ‘it: -_ 
with useful paccuracy. 
_ Hoop Tensions. —Perhaps the most interesting data are those arising from 

te measurement of strains in the sheet piles. Fig, 19 facilitates comparison 

of the hoop tensions for the different fill levels and at different points | around © q 


4 the cell. The lower tensions for Pile E may be because the natural ‘ground lev-_ 
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el there (Fig. , 3) was higher than it was on the south side of the cell. For the the 
cell full (Stage IIT) condition, the tensions were otherwise ‘comparatively uni 
form around the cell, except for Level 4. 
_ The highest hoop tension for Stage II (cell 4/10 full) \ was 2.8 kips per per lin in. 
: - Level 4, Pile D, as shown in Table 5. 4 oo only one set of readings was | 
obtained for this ‘condition, the same tension appears in Table 6, When the cell 
had been filled, several sets of readings were taken, and the highest individual © 
tension obtained was «6. 1  kips per lin: in. (at Pile D), with a corresponding av- 
- large local radius of curvature (the cell wall was almost straight), or to the 
high rate of filling at that point which occurred when the soil flow under the 
|e 20 reveals a systematic pattern of development of horizontal tensions 


‘The rosette installation at Pile C were made w with the thought of studying the — 
complex structural situation near the T- pile connection between cell and arc. . 
_ The data were examined carefully in this: regard, but the results were incon- 


, Sets of strain readings were taken at short intervals in 2 ‘days to “study the 
affects of tidal changes on hoop tensions. . The a av verage horizontal stresses in- 


a Only a few values of hoop tension were botained for the p pilot cell, Cell: “ 
_ These were of low magnitude, , because the strainmeters were approximately 
_ 7 ft below the natural ground, where the cell wall was not as free to expandas _ 
it was at points intrumented onthe maintest cell = = 
Hoop Tension as a Function of Fill Height.—Fig. 21 indicates for Levels 2, 
: 8, and 4 an almost linear distribution of the average hoop tension in Cell 16 
_ with regard to elevation for the Stage II condition (cell 4/10 full), and a simi- 


7 lar ar trend is s is shown for Stage III (cell full). , The stresses obtained for Level 1, 


“however, are re much higher than would be obtained by linear extrapolation trom 
‘Lev els 2, 3, and 4. The high average stresses at Level 1, Stage III, are in- 
fluenced by" the high | individual stresses: at Cc and A. The fill was ‘slightly 


It is 


_ Fig. 3 shows that the top of the fill { for toe Stags II was between _—" 2 and 3, 
nowever, Fig. 21 clearly indicates a positive average hoop tension at Level 4 
_ which was 13 ft above Level 2. The curve for Stage III in Fig. 21 seems to im- : 
= ply an ‘appreciable | hoop tension at the top of the cell for the cell-full condition. > 
“3 To study this point further, Fig. 22 was drawn by correlating fill- -height 
above meter level with strain measurements at Levels 3 and 4, In Fig. 22, | 
are included for fill levels intermediate those of of Table 5, The 
_ heights above meter level were insufficiently large at Levels 1 and 2 to justi- 
__ fy graphing the data at these levels, It is confirmed by these curves that ten-— 
_ sions exist at a given level before the fill surface has risen to that level, ‘The — 
i dashed line in Fig. 22 represents the trend of the relationship for the greater 
_ ‘fill heights. The line passes through the origin, susgesting that the relation- 


ship becomes linear after the fill surface has risen n appreciably above the me- 


ter level 


tet 
3 
rat uP 
2 
F 
— 
— 
é 
Gall | — 
— 
> 
= 


This phenomenon may be examined aaiiiiutt if it is assumed that the © 
lateral earth pressure increases in proportion to fill depth and that the cell — 
behaves elastically. The pressure distribution is as indicated in Fig. 23 <a 
If a horizontal cut in the cellwall were made at the fill surface, the lower por- 
tion of the cell would expand as in ‘Fig. 23 (b) and develop the —_—— 


in which @ is eo rotation 1 of each cell element in in 1 the vertical plane; K K denotes 
the ratio of horizontal to vertical earth pressure; y represents the soil unit | 
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1 HOOP IN KIPS PER LIN. INCH 


—AVERAGE HOOP TENSION VERSUS 


‘weight; R refers to the cell radius; nisa coefficient depending on eheaiiniie 


ty of sheet pile interlocks; -E represents me modulus of elasticity | of the — 


alent to a modulus of subgrade solution is corrected for the ef- 
fect of the cut 7 pr oviding edge moments (Fig. , 23 - and edge sh shears (if re- . 
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CELLULAR BULKHEAD 
quired) to eliminate the physical discontinuity. To solve for Mo, the redundant 
edge moment, one needs the influence coefficient 12 for edge rotation due to a ; 
unit — edge moment. The edge moment er to eliminate the discon- 
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tinuity in 1 slope is ‘obtained by dividing by the influence coefficient. Fig. 24 
‘shows the theoretical | hoop tension | curves: thus: obtained for two values of the 
while that for k = 7 includesa vinnie allowance for interlock extensibility. 13° 
This analysis verifies qualitatively the implication of Figs. 21 and 22, but | 
the observed effect considerably exceeds that shown in Fig. 24. The difference 
might be reconciled by considering the the analysis the transverse panei of wi 


Coefficient of Lateral Earth Pressure.—From Table 5 , the hoop anaes 7 
“hight ee at Levels 1 A 2, , 3, and 4 for the cell- full condition (Stage II) are found 


assuming in accordance with the analysis leading to Fig. 24 that the horizontal 7 

loading discontinuity effect is zero when the cell is full. The lateralpressures _ 

shown in Fig. 25 (a) are thus obtained. _ ‘The v; value at top of cell is shown as_ 

_ In accordance with the data of Fig. 14 and Table 3, the high . cell water level 

_ and the lateral water pressures and seepage forces from the flow net must be _ 
_ considered. The computation is made for the case where the water outside the 

cell is at high tide, 5.5 ft, because the water levels are fairly stationary at 
time. net t lateral water pressure, which is is from the total 


aa 12 “Formulas for Stress and ‘Strain, *by R. J. Roark, McGraw- Hill Book Co., » Thir 
13 “Cellular Sheet Piling Structures—A Detailed Study of the Stressing in Circular 
walls,” by Descans, L’Ossature Brussels February, 1954, 
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60 Ib per cu ft above and below wa- 

ter level, respectively a are used in . the ‘computation of vertical soil pressure ~ 

Po 25 (c)). _ The seepage force of the downward flowing cell water (Fig. 25 7 


‘(d)) is added. 


4 


& 


_ VALUES ARE FOR 
SEPT. 27,1952, 3PM. 


A) MEAN ma-)) UNBALANCED HORI- VERTICAL SOIL PRES- VERTICAL 


FROM HOOP ‘ZONTAL WATER SURE FROM SOIL UNIT SEEPAGE 


AND VERTICAL PRESSURES, STAGE III, HIGH 


The coefficient of lateral earth pressure is obtained as as eon average of the 
quotients of the net horizontal and vertical soil pressures at the four datum 
levels. this method, it is found, for Ill, that 


- Approximately t the same value was obtained from a similar computation using 
the Stage III average hoop tensions with outside water level taken at mean sea 
level. For design, the engineer must consider the total lateral pressure, not 
just the portion due to soil as determined herein. 
_ During the early stages of filling the cell, the water levels inside a and out- ; 
7) side t the cell were the same. Thus it is probably acceptable herein to neglect | ; 


14 A Soils Study for E Test Cell No, 16,” L. Evans, Long Beach Har- 
bor Dept., , January 3, 1955, 
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‘the seepage and consolidation hydrostatic in the coefficient 


of lateral earth pressure . Fig. 22 may be utilized for this purpose, taking | a 
as inversely proportional to the slope of curve. The proportionality con- 


are shown in Figs. 15, 16, and Be, ‘and in . Tables 2 and 4. . The indicated verti- 
7 es movements over the period from late 1952 to 1954 are subject to correc- 
tion for the effects of regional subsidence , but the horizontal mov movements and 
rotations were not significantly influenced by these effects. 

_ The movement of the cell as a whole, as illustrated in Figs. 15, 16, and 17, 
showed nothing significant through Stage V, where the backfill was raised to 

el 12. _ ‘The surcharge load, however, carried to el 25 or el 26 a 


"effect was a southerly tipping of the sheet piles of more than 3°, which was 
~ contributed to by a 1° southerly rotation of the cellas a whole. ‘The net cell 
_ racking, more than 2° , created relatively large shearing distortions of the soil — 
_ within the cell. The reason why the movement stopped at this point is the con-_ 
4 sequent reduction, , toward the active value, of the horizontal pressure of the 
backfill against the cell. Note from Fig. ‘7 ‘and Table 4 that other deep-water 
cells also experienced exceptionally large moments away from the backfill. 
_ The soil in the cell also settled with respect to the steel, between Stages IV 
and VI, by 0.4 ft at the cell top and 0.3 ft at Level 2 (from Table 2 and Fig. 16). 
_ Hoop Tensions .—Figs. 19 through 21 and Tables 5 and 6 showthe effects on 
hoop tensions of the events subsequent to filling of the cell. Computation of 


' coefficient of lateral earth pressure after filling behind the cell and in the 
arcs, and after adding the surcharge, is difficult to do with confidence and pre- _ 
-_ cision. _ The asymmetry of fill- loading and the forces applied to the cell cell att the 


7 Soil unit weights!’ of 116 and 60 lb per « cu ft above and below water level, 
respectively, were used, and it was assumed that both cell water level and sea 


level water at el 4.3 (approximate | mean sea level adjusted for regional ‘sub 


: face, however, the surcharge should be more directly reflected in hoop a? 

: sion. Because of the wnpgvtngrerd as to shear flow between steel and soil around 

7 the cell circumference, the lateral earth pressure is computed by using the 
south face hoop tensions at Piles A and C only. This leads to an approximate — 
value of 


a 
am og using the average (dashed) line in Fig. 22, The value in Eq. 5 would be cor- _ 
rected downward if pore pressure due to consolidation were large. 
ie = _ Settlements and Deflections.—Movements of the cell and the soil under the g 
| 
— 
— 


which is essentially the same as in Eq. 4 Apparently then, , subject - still to the Z 
_ theoretical limitations previously mentioned, the large cell distortions did not — 
greatly affect the hoop tensions, and the hoop tension increases between Stages © a 
‘Il and VI may be assignable primarily to the increase of fill depth associated - - 
withthe surcharge, 
The largest hoop tension observed in the test, approximately 7 kips per lin 
in., occurred at Level 3, Pile A, in Stage VI (Fig. 19 and Table 6). Tensions 
at Piles A, C, and D generally exceeded those at Pile E, especially during the 
_ latter phases ‘of construction. This pattern can be explained in terms of — 
existence of soil pressure on both faces at Pile E and of the comparatively _ 
large water depthat PileA, 
rr The comparatively small hoop tensions at Pile B (Table 5) apparently fol- 
low from the facts that the connecting arc radius (Fig. 3) is much smaller 


the cell radius, the in the arc is confined ifined in at bin-like 


SE SELECTED FINDINGS 


The t test data have been presented in as complete detail as was feasible, in 
order that others may use them appropriately. . The following summary em- 


phasizes | some of on significant findings as developed herein. 
; ; The coefficient of lateral earth pr essure was ( btained das 0. 66 di the 
pressure was obtained as uring 
filling of the cell, 0.54 shortly after the cell was filled, and 0.53 2 yr later af- 
ter application of a 14-ft surcharge. | The apparent independence of the latter 
values with respect to large celldeformation is an important point. Unbalanced — 
pore pressures must also be included in computation of hoop tension. Bs, 
2, +The large movement accompanying application of the surcharge con-— 
sisted of a 1° rotation of the cell as a whole and an added « shearing distor - 
tion of the sheet piles and the contained soil. 
3. The maximum hoop tensions occurred at Level 4, approximately at the 
original ground level. The largest individual hoop ‘tension determined was — 
7 kips per lin in. ,at Level 3, PileA,StageVI. 
4, The soil in the cell required 10 days after completion of filling to reach 
nities 90% consolidation. During and before this time, pore pressures 
due to consolidation reached values of the order of 1.2 ftof seawater, —__ 
kh The filling procedure used in construction is a major factor in cell 
-alinement. _ Specifications should require adequate outside soil depth before > 
‘filling; the symmetrical depositing of fill in the cell, and the arising of cellwa- 
ter level to the top before the fill rises more than + to + of the cell height. The 
latter measure tightens the interlocks and stiffens the cell during filling. 
6. Hoop tensions where necessary were satisfactorily determined from 
"measurements with horizontal strainmeters a alone, with an appropriate /correc-— 


. This implies p possibilities for simplifica - 
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‘Upward movement or “pile heave” commonly occurs whenever large 
wan of closely spaced piles are driven. Pile heave is generally attributed | 
to Captacement of the foundation soils during subsequent driving of adjacent. 
piles. The majority of engineers appear to be of the opinion that all heaved | 
piles must be redriven. Very few documented case histories dealing with this 

Problem are available in the engineering literature. 

This paper presents detailed field data with stone pertaining to pile- 

a ‘heave problems that developed during construction of end-bearing pile founda- 
tions for a large industrial plant. Two individual areas within the plantsite are 
selected for detailed study. Pile heave occurred at both areas, However, _ 
whereas in one area the heaving constituted a serious foundation problem 
necessitating the redriving of all heaved piles, in the other area it presented 
: -no problem and redriving was not necessary. The successful performance of 

the completed pile foundations at both areas is confirmed by settlement ob- 

_ Servation data maintained on the buildings since construction, = 
‘The paper shows that each pile-heave problem must be evaluated on its © 
Own merits. In order tocarry out such an evaluation, adequate data concerning _ 
_ - items as subsoil profile, magnitude of pile heave, magnitude of ground © 

“heave, pile penetration diagrams, and pile loading tests must be gathered and - 


Note.—Discussion open until January 1, 1962, To extend the closing date one month, | 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the American 
Society of Civil Engineers, , Vol. 87, No. SM 4, August, 1961, 
4 1 Partner, Ripley, Klohn and Leonoff Ltd., Consulting Engineers, Vancouver, B, c. 
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- May be ¢ determined and consequently its probable ¢ effects on the performance of 
the structures res evaluated. 


‘The’ paper is based on data obtained during construction of a large expan- : 


sion to an existing pulp mill located at tidewater on Canada’s West Coast. A 
plan of the development showing both the original and expanded mill layouts is. 
presented on Fig. 1. _ The original pulp mill, built in 1946, was supported on 
timber end-bearing piles driven into underlying glacial outwash deposits - Ap- 


--90--- Contours indicate average tip elevation 
att which timber piles met refusal _ 


-- Location of simplified section through 
"foundation soils-see Fig. 6 
Log pond Reference El. zero (0) is 100 
below mean low tide 
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proximately 6,000; piles, most of which were 35 ft or less in length, were driven 
for support of the original mill. Settlement observations maintained on the 


- had performed satisfactorily and that no foundation settlement problems i 
_ developed. (Settlement records indicated average building settlements sranging — 
0.02 ft and 0.03 ft for the period 1946 to 1955.) 
The 1955 mill expansion program involved the construction of several anes | 
and | heavily loaded | buildings. , Support of these structures required the driving | 
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1 Building observation points A, B, Cc 
measurements show no movement 
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a “FIG. 2. —GROUND SURFACE MOVEMENTS DURING G DRIVING | OF 


STEEL TEST PILE GROUP 


of inna e groups of piles of greater length and at closer spacing than had — 
viously been necessary for the original mill. Approximately 26,000 timber end- 
_ bearing piles” ranging in - length from 10 ft to 125 ft and having an average 
length of about 50 ft were required. In many instances the new structures joined 
‘directly onto existing mill buildings containing equipment sensitive to move- 

ment. As this a a had to be kept 1 in continuous Operation throughout the 
vented, 
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Construction of pile foundations for the new ‘buildings therefore had to satisty a 


achieved for the new mill structures. __ 


first requirement, prevention o of damage to ex existing — 


ay) Th ‘They permitted use of a higher allowable load per pile (50 tons), there- 
by reducing the number of piles in a group from 16 timber to 6 steel. me On 


(2) Being low displacement piles a group of them displaced a relatively 


ae 
small volume of of soil, 
‘The aamiienis - use steel H piles was partially based on data obtained from 
test group of 5 steel piles driven adjacent to an existing mill building. Brief- 
ly, these data (presented in Fig. 2) showed that the existing building was not _ 
affected | by driving the test-pile group and that subsidence of the ground sur- — 
: te did not exceed 0.01 ft beyond a distance of 10 ft from the centre of the 


exterior, walls” of existing buildings and on n sensitive equipment foundations. 
within the buildings, Level observations were maintained on the bench marks 
_ throughout the driving of the steel piles. These observations indicated that 


a distance of approximately 10 ft from existing buildings. These piles offer offered 


the: steel piles had no measurable effect on the existing buildings 


cept at one local area where maximum settlements of 0.01 ft were reported. 

The steel piles were therefore considered a success, having performed their — 
design function of limiting the soildesplacementtoa value that would not affect 
the existing adjacent pile foundations. A photograph showing steel piles driven 
next to the pulp machine room is presented in Fig. 3. All steel piles were 
A driven before driving of timber piles was started. This figure also illustrates 


the difference in size between steel and timber pile groups. Fig. 4 shows - 
ground heave for a large timber pile group and indicates the kind of soil as- 
placement problem that steel piles may have prevented alongside the buildings. — 
_ ‘During the initial stages of the driving, several steel H piles — 
very high penetration resistances at relatively shallow depths and, before 
the problem | could be assessed, the field forces attempted to achieve — 
penetration by continuing to hammer thep piles. The attempt appeared success-_ 
oe as the piles were driven several feet deeper. It was noted, however, that 


the upper ends of the piles nae both rotated and moved laterally. To basened on 


damage. Fig. 5 presents photographs of two such piles, 
_ The extensive damage caused by overdriving is quite apparent, as is the 
corkscrew shape of one of the pile tips, a shape that explains the pile rotation 
that developed during driving. As a result of these findings, all overdriven -_ 
in this local area were pulled and replaced. 


free ‘from harmful ‘settlements, an evaluation of the effects of pile heave ons 
future foundation behavior was required, Carrying out such an evaluation con-— 
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PILE HEAVE- 
ee the major foundation consideration at this site and the following sec- 


tions. of the paper | deal exclusively with this item. as 
ae 
GENERAL LAYOUT AND SOIL ¢ CONDITIONS 
A general layout of the pt pulp and paper mill showing both the original and 7 
the new mill buildings is presented in Fig. 1, Also shown in this figure are > 
contour lines depicting the average tip elevations at which the timber end- 


Kilns Shops and 


stores 


Ablation till —7 large inclusions of 
dense hard glacial till 


(a) SECTION A-A Screen and bleach 


Mill SNF Pile cut-off El. 110 ft 


inclusions of dense 


Kiln Mill El. 119 ft 


Sands and 
an san gravels with large ey 


inclusions of dense, hard © 


SECTION C-C_ 


pearing piles met refusal. These contours indicate that the mill straddles a 
buried spur of hard material that slopes downwards in a westerly direction, 
Geological Opinion is that the spur was carved by erosion from a thick sheet 
of glacial outwash material. deposit, “into which all end-bearing piles 
were driven, is extremely variable, as proved byt both the foundation explora- — 
tion and construction programs, It consists mainly of medium-to-dense sands 
and gravels. However, it does contain large inclusions of dense hard glacial 
‘till and an erratic scattering of boulders, In one area extensive pockets of 
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clean loose sand were noted. In some areas the outwash material is ditty 
overlain by recent deposits of silt andclay which extend to the present ground 
surface. In other areas a soil deposit, referred to in this paper as “ablation . 

till,” occurs between the top surface of the outwash material and the overlying © nq 

silts and clays. The ablation-till deposit consists of a well-graded mixture of _ 
particle sizes ranging from clay to fine gravel. The material has the appear- - : 
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- During the last retreat of the Pleistocene ice sheet, sea levels rose and E. 
flooded the giacial deposits, burying them beneath fine grained marine sedi- 
ments. With increasing depth below the present ocean floor their grain size 7 
“decreases” from silt to clay size. The silts - non-plastic and have natural 7 
water contents. ranging between 35% and 45%. The clays have natural water 
- contents ranging between 30% and 50%, liquid limits of 40% to 50%, and plastic 
limits of 18% to 22%. The natural water content of the clay is usually sli slightly 

_ Typical sections through the foundation soils at this s site are shown in. Fig. 7 
? ‘More detailed subsoil data relating to particular pile problems are nted 
in those sections of the paper in which th the are 


a ‘The term “pile neene? ’ as as used inthis pap paper refers to the upward mot movement 

of a previously driven end-bearing pile caused by the subsequent driving of ae 
additional piles in the immediate vicinity. Apile is heaved by the upward com- 

ponent of movement of the soil displaced by these subsequent piles. The mag- — 
nitude of pile heave to be expected at this site was not known beforehand, A 
search was made of the engineering literature in an attempt to find data that —_ 
would provide some ‘guidance in evaluating this problem, Several references 

(1) (2) were located; although these did not treat pile heave in extensive 
detail they did indicate that heaved piles must be redriven or settlements ap- 

proximately equal to” the “original heave » might occur under the completed 
building loads. This opinion was also reflected by the requirement of most 
building codes (4) (5) (6) that heaved piles must be redriven, As the mill — 
structures are extremely sensitive to settlement the designers considered that 


maximum allowable pile heave of 0.03 ft was tentatively selected. For piles 


_ settlements in excess of 0.03 ft would be intolerable. Using this criterion, = } 
that heave in excess of this figure redriving would be considered er , 


Pile heaves far in excess of the maximum tolerable settlement 
©. 03 ft) occurred as soon as pile driving got under way. Obviously, unless 
under: future applied building loads, -redriving of most piles ‘would be necessary. 
‘The added costs and the additional time involved if redriving was — 


‘would be most critical for this project. ‘Therefore procedures had to be evolved 


Effects of Pile Heave on End-Bearing Piles.—For end-bearing piles that 
derive a large percentage of their allowable bearing capacity from point re- 
sistance (for example, ‘those driven through soft highly compressible soil to 
rock), pile heave during construction can be the cause of excessive ace 
= of the pile foundations at some later date when they are loaded. This problem 


2 Numerals in refer to corresponding items in the Bibliography. 
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bearing ype ta _ When such piles are loaded, they will settle until their tips 
are once again in intimate contact with the bearing stratum. Thus, settlements 
approximately equal in magnitude to the original pile heave may occur, 
the other hand, for end-bearing piles that appreciably penetrate the 
bearing stratum and derive a large percentage of their allowable bearing 
: capacity from skin friction developed within this stratum, pile heave may not 
pose a future settlement problem, For example, if the frictional resistance 
developed within the bearing stratum by the embedded lower portion of the 
pile exceeds the upward drag exerted on ‘the pile by the soft overlying sedi-— 
ments displaced during the driving of adjacent piles, then these upper soils 
will heave without moving the pile. But soil displacement will also occur within 
BS the bearing stratum itself, due to the pile penetrations. Such displacement may 
result in appreciable pile heave; however, as this heave is caused by a — 
upward movement of the bearing stratum ‘soils it ist unlikely to result in a 
y decrease in the resistance that the pile candevelop within the bearing stratum. 
Therefore future foundation settlement due to pile heave is unlikely to be a 
- problem provided the heave is caused by displacement within the bearing stra- 
tum as opposed to displacement within compressible soils overlying the bear-— 
—ingstratum, 
Both types” of pile heave were encountered at this site, and an example of 
each has been selected for analysis. At Area A, pile heave is attributed to 
soil the bearing stratum. At Area B, it is 


_ Area A.—The boundaries of Area A are outlined on Fig. 1, a _—— plan 
of the millsite, Pile driving, pile-load testing, and subsoil data pertaining to — 
this area are presented in Figs. 7 to 10 inclusive. A photograph showing pile- 
driving operations in this area is presented in Fig. 7. Typical pile penetra- 
tions diagrams are presented in Fig. 8. These diagrams illustrate the general - - 
pattern of pile resistance during driving, and show how this resistance builds ~ 
up as the pile penetrates the bearing stratum. All piles were of Douglas Fir) 
and were driven using a No. 1 vulcan hammer, In Fig. 9 are presented the ~ = 
_ results of 4 pile-loading tests made on representative end-bearing foundation 
7 piles, Pertinent information relating to. subsoil profile, pile dimensions, meth- 
od of test, etc., are noted in the figure for each individual pile test. In Fig. 9 / | 
all test piles were of Douglas Fir, The piles were driven using a No. 1 vulcan 
oe ‘single acting steam hammer. The rate of loading for pile tests equals one load 
- increment every 15 min, Permanent pile movement was measured by unloading 
_ the test pile to zero load og recording the permanent settlement the previous 
load had produced, In Fig. 10 are presented the results of a series of pile- 3 
loading -* carried out on a timber friction pile driven into the soft marine 
deposits. © The rate of loading for the test was one increment of load every 
hour. Also shown in Fig. 10 are field and laboratory test data obtained from 7 
a drill hole made approximately 3 ft from the test pile. The test pile and drill 7 
hole data are considered representative for the marine deposits encountered 
a review of the data obtained on each individual pile during driving indi- 
cated that the majority of foundation piles had heaved as adjacent Piles were 


> “most piles had heaved between 0. 10 ft and 0.20 ft. Settlements of this “magnitude 
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in the completed building were intolerable and unless assurance could be given 
that the heaved piles would not settle, redriving would be necessary. Such an 7 
operation would be both costly and | time | and \ to be 


ablation the sand and gravel before 
refusal, The marine deposits that overlie the ablation till are very soft and 
a offered little resistance to penetration of the piles, The frictional resistance 


that the piles would develop within the was judged to be greater 
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‘FIG, SUBS PICAL PILE PENETRATION DIAGRAMS WITH SIMPLIFIED 


parently caused by displacements within the ablation till and underlying sand 
and gravel soils. If this were the case, then intolerable ble settlement of the 
heaved piles at some future date would not occur. 
To check this appraisal of the heave problem several pile- loading tests” 
were run, All test piles produced load-settlement curves similar to those _ 
shown in Fig. 9. ‘Approximate. ‘computations were carried out to determine 
- what portion of the total load applied to the test pile was carried by skin 


friction developed in the soft marine deposits and what portion by the lower 
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end of the pile in the stratum. ‘Ultimate skin friction values 
“for the marine deposits were obtained from the friction pile-loading test (Fig. . 
10). In every case the computations indicated that the load carried by the sec- 
tion of pile embedded in the bearing stratum was from 2 to 6 times the design 7 
pile load of 20 tons, Admittedly such computations are approximate only, as 

pooh conditions at the test pile location may vary somewhat from those 
assumed for computation purposes. However, when viewed in n conjunction 
all the available data these computations were considered to confirm the p pre- 
heey appraisal that pile heave would not presenta future settlement sl 
lem, Foundation piles were therefore cut off and capped without redriving, 
Settlement observations maintained since construction of the mill expansion A " 
— (1956) confirm that pile heave did not present a potential settlement a 
in Area A, The average recorded settlement to date (July 1960) for the ware-_ 
a building is 0.01 ft. _ The initial settlement readings were made whenthe _ 


: 


building was partially constructed and while an estimated 30% of the final live- 
plus-dead loads was acting on the pile foundations, For the completed structure _ 

_ the average live-plus-dead load acting over the entire plan area of the building 
is approximately 2000 psf. As would be expected, the warehouse structure has > 
readily absorbed these very small settlements and is in excellent condition. 

Area -B.—The boundaries: of Area B are outlined in Fig. a general 1 plan 
of the millsite. Pile driving, pile testing, and subsoil data ‘pertaining to this 
area are presented in Figs. 11 to 14 inclusive. A photograph showing pile - 

driving operations in this area is presented in Fig. 11. Typical pile diving, are 


: diagrams illustrating the e general pattern of pile resistance » during driving, are 
: “presented on Fig. 12, All piles were of Douglas Fir and were driven using a No. 
— 1 vulcan hammer, In Fig. 13 are presented the results of 3 pile-loading tests 
‘made on representative end- bearing foundation 1 piles. Pertinent information 
_ relating to subsoil profile, pile dimensions, method of test, etc., are noted on 
& 4 the figure for each individual pile test. Alltest piles were of Douglas Fir. The 
Pr piles were driven using a vulcan No, 1 single acting steam hammer, The rate 
of loading for pile tests equals one load increment every 15 min permanent a 
- pile movement was measured by unloading the test pile to zero load and re- q 
cording the permanent settlement the previous load had produced. In Fig. 14 
are presented a summary of field and laboratory data for 3 test holes located 7 
Pile heave observation data obtained | on each individual pile in Area B 
clearly indicated that large” , heaves occ occurred during the driving operation, _ 


J Most piles apparently heaved between 0.3 ft and 0.4 ft. However, large varia- 


- tions in heave were recorded for individual piles within any one pile group and 
‘maximum m pile heaves ¢ of 1 ft or greater were noted in several instances. The 

j large pile heaves were also accompanied by appreciable heaving of the ground 

_ surface in the vicinity of each pile group (Fig. 4.00 2 
~~ The pile driving data for Area B, particularly the pi pile-penetration diagrams, 
“definitely showed that the behavior during driving ofa large percentage of the © 


bearing in this area appreciably from that of the piles pre- 


‘The overlying marine deposits of! offered appreciable 1 resistance to driving. 
7 -@) The piles did not slowly build up t resistance to further penetration; in- 
stead they refused 


: 
: 
if 
= 
| 
= 


a 


4ad SMOjq Ol 0 


ul 


AVE 


@ 


PILE HE 


| 
=. 


o asueg 


ul 


WS 


$de} May wapun 
payesjauad aouejsig 


| 

| 


These « observations clearly indicated that the subsoil profile at Area B was 
dq appreciably different from that at AreaA. The available test-hole data implied 
that the overlying marine deposits in this area were medium-to-stiff clays 
rather than silts. ‘The | bearing stratum appeared to be a dense glacial till. iz: © 

To check these items, 3 additional test holes were immediately drilled. These 
EE confirmed that the marine deposits at Area B were medium- to- 
_ stiff clays They also confirmed that t the bearing stratum was a dense hard 
_ glacial till into which the timber piles ; could not penetrate more than one or 
~ two feet. In the area that exhibited large pile heaves, the overlying clay ap-— 

peared to rest directly on the dense, hard glacial-till bearing stratum, — 
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_ eueen foundation problem, It was reasoned that those piles that had met sud- | 
den refusal, without appreciable penetration into the bearing stratum, had 
offered little resistance to the upward drag exerted on them by the displaced 
clay and that consequently they had heaved, lifting their tips out of contact with 


the bearing redriven until their tips v were again in sound 
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August, 1961 _ 
To check this appraisal of the pile heave problem and confirm the necessity 
for redriving the heaved piles, , several representative piles were selected ll 
pile-loading tests. Each pile was loaded until it reached what appeared to be _ 
_ its ultimate bearing capacity. At this constant loadthe pile moved downwards. 
7 ‘However, after having ‘Settled a a distance that ranged f from 0.1 ft to 0.4 ft for 
different test piles, the pile stopped moving ‘under the constant load. 
of the test pile was aspen continued until it again moved downwards under con- 
Stant load. In every this final or ultimate load that caused 
movement of the pile appreciably ¢ exceeded the one that caused the initial pile 


movement. A suggested explanation for this observed behavior of the test — 
—_ 


7 ai Initially, a a ai was driven through the overlying marine clay and met 
practical refusal on contacting the dense glacial till, As a result, the pile 
t tip did not penetrate more than one or two feet into the tilllayer,. | 
(2) _ When subsequent piles were driven in the same pile group, they dis- 
pk aced large volumes of the” overlying clay soil so that the ground heaved © 

The displaced clay exerted large upwai upward forces’ on the. previously 
‘driven pile. These forces exceeded the resistance to upward movement that 
the pile could develop with its small penetration into 1 the bearing —_— 
Consequently the pile moved upwards. 

_ (4) In moving upwards the pile tip lost contact with the bearing wonene: - 
and the pile in effect became a friction pile supported by the sii =) 

ae (5) When load tested, the pile behaved initially as a friction pile suppo es 
in the overlying clay, and when the ultimate skin friction value of this mate- | oF 
rial was reached, t the ‘pile failed, moving suddenly downwards until the tip 
again came into contact with the bearing stratum. With the pile tip incon- _ 


6) “The ultimate failure load the ultimate total capac- 
ow! of the pile, equal to skin friction in the the overlying clay plus point resist- 


: ‘The pile- loading test data were interpreted to confirm the original ap- 
i - that pile heave at Area B posed a major foundation problem at this — 


site. Displacement of the overlying marine clay during pile driving had ap- 
parently caused the pile tips to lift off the bearing stratum, In their heaved > 
- condition, the foundation piles were, in effect, friction piles supported by the _ 
_ overtiyi ing clay. As such they would undergo withtime large settlements under 
the applied building loads. Settlements of this kind would be intolerable for 
propose d structures. _Redriving of all in this area was therefore 


mee and close supervision, On the one han hand, heaved piles had to be re- 
driven sufficiently hard that they moved downwards a distance at least equal 
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to the observed heaves. On the other hand, excessively hard that 
might cause brooming or breaking of the piles could not be tolerated. 
| _ The shorter the waiting period between driving and redriving, the a _ 
_ would be the skin friction forces that must be overcome during the redriving 
operation. The available loading-test data indicated that 12 hr after driving _ 
the piles had developed 75% of the ‘Skin friction they would 
“develop three weeks after driving. Beyond 3 weeks, the rate of increase of 
ultimate bearing capacity with time was estimated to be slow. Obviously, 
unless the heaved piles could be redriven within a few hours after the initial — 
driving they would have to overcome large frictional resistances during the 
‘redriving operation. Unfortunately on this project, “obtaining the required pile 
heave observations, and redriving the piles, required a minimum time interval 
2 approximately 48 hr. This delay was caused by several factors, including & 
_ ment in use of the site, and the construction procedures that had been de- 
veloped for driving piles over the site area. This time lag between driving 
redriving made redriving operation difficult than 
‘In order to develop suitable criteria for redriving heaved piles, several 
representative piles in the area were selected for testing. The piles selected 7 


_ for this purpose had originally been driven several days previously. Briefly, 
_ the driving tests consisted of applying 10 blows to the pile with a No. 1 vulcan 
_ hammer (5000 Ib ram, 15,000 ft lb per blow) and then measuring the penetra~ 
tion achieved. All test piles were redriven to a minimum penetration resist- 
= of 1 in. for 10 blows and some of the piles were redriven until at least — 


three consecutive sets of 10 blows had each produced less than 1 in, of = 


tration, Several of the test piles were then pulled and examined for damage. 
= In every case those piles that had moved 1 in, or more under the first 10 
blows of the hammer and had subsequently been redriven w a final penetration 
resistance of 1 in, per 10 blows had undamaged tips. ‘Even those piles that an 
had been subjected to exceptionally hard redriving forces showed only moder- 
ate tip damage. Brooming and splitting of the pile tops appeared to be the 
more critical problem. ‘This behavior was considered normal since a large 
‘percentage. energy applied to the pile heads during driving never 
_ reached the pile tip, but instead was dissipated by the frictional 
™ pile had developed in the overlying clay soils. | ‘ee: eye 
- After reviewing the results of the pile redriving and pulling tests om 


a (1) Unless the e pile ‘moved 1 in. or or more under the first 10 blows of of the 
pile hammer, the second or third set of 10 blows still did not produce move- 
ments in excess of 1 in. per set of 10 blows. In addition, the exceptionally 
hard driving was likely to damage or break the Pile 
(2) Piles which moved a distance of 3. in. or more under the first 10° 
blows of the hammer normally could be « driven several inches before tighten- — 
7 ing up again. This tightening up was interpreted to indicate that the he pile - 
once again in contact with the bearing stratum. 


oa (3) Piles which moved 1 in. or more under the first 10 blows of the hammer 


and were subsequently driven to a final penetration resistance of lin. per | 
10 blows would not have their — damaged d by th the seeeteanes developed in the 
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“tort the basis 0: of these conclusions, the: following procedures were developed — 
for redriving the heaved piles: 

(1) The pile was given 10 blows with the No. 1 vulcan hammer. it the ‘pile - 
Pa.) 1 in. or more, it was given an additional 10 blows. > 


 . This procedure was repeated until a penetration of less than 1 in. for - 


ss The boundaries of the critical areas where pile -redriving 1g was considered 
on the basis of a study of both the pile-heave 
observation data and pile-penetration diagrams. Generally, pile -re-driving 
ye 
was extended over those areas where the pile-penetration diagrams indicated — 
_ that the piles had met sudden refusal without appreciable penetration into the a 
bearing stratum, All piles in the areas selected for redriving were redriven, _ 
regardless of whether or not the observation data indicated excessive heave. 7 
_ The redriving areas were purposefully made to overlap adjacent areas where — 
4 pile heaving be yt not | appear to be a problem, so as to. ensure > that n no 0 piles” 


penetration was achieved for the first 10 blows in areas where pile heaving 
was considered critical , whereas less than 1 in. of penetration for the first 
10° blows was obtained in those areas bordering the critical area that — 
purposefully overlapped during redriving, 7 
At most locations a direct correlation between observed pile heave ~ 
_ observed pile penetration. during redriving (pile redrive) was not evident for — 
— each individual pile. However, a review of the heave data in terms of average 
values per pile _ group, indicated that in areas considered critical for = 
heave the average redrive value approximately equal to the ll 
{ heave value for the same pile group. These data were interpreted to indicate 7 
that the redriving operation was successful and had accomplished its Objective — 
7 of driving back down into the bearing stratum those piles that i been un- 
tes Settlement observation data maintained since construction of the mill | a 
expansion (1956) that the redriven pile foundations have performed 
_ satisfactorily. _ The average settlement to date (July 1960) for the screen 
_ 7 and bleach plant, which is located in Area B, is 0.02 ft. The initial settle- 
ment readings were made when the building was partially constructed and 
_ while an estimated 40% of the final live- plus- dead loads was acting on the 
La pile foundations. For the completed structure the average live-plus-dead | 
load acting over the entire plan area of the building is approximately 3,000 
Ib per sq ft. As would be expected, screen and bleach plant has 
absorbed these very small all settlements and is in excellent condition, ee 


|. 
CONC LUSIONS" 
wei are many aspects of foundation behavior that are not ‘subject to 
exact pre-determination by theoretical analyses. Pile heaving, andthe effects 
_ Of pile driving on existing adjacent buildings, are two of these. Evaluation of _ 
such problems requires the gathering and studying of detailed field — 
tion data as well as a knowledge of the theoretical principals involved. 
‘The effects of pile heave on the ultimate performance of a pile foundation — 
a depend on the mechanics of the heaving action . Each pile heave problem 
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must be individually evaluated. To carry out such an evaluation, adequate : 
information - concerning such items as subsoil profile, magnitude of _ 

_ heave, magnitude of ground heave, pile penetration diagrams, and pile loading» 
= must be gathered and studied. The data presented show that whereas in 
~ gome instances pile heave may constitute a serious foundation n problem neces- 

: _sitating the redriving of all heaved piles, in other instances it may present 
no problem and redriving is not required. 

Where are driven adjacent foundations, move- 


| isplacement to values that did not affect the existing foundations, — 


Owner of the ae and pa paper r mill project. is MacMillan a: and Bloedel Ltd. 
The project was designed and constructed under the supervision of H. A. 

_ Simons Ltd., Consulting Engineers, Vancouver, B.C. Senior soil mechanics | 
consultant was Terzaghi. The writer’s firm acted as general soil 

chanics ‘consultants to the designers and in this capacity worked in close 
cooperation with Mr. Terzaghi. The writer is indebted to: Mr. Terzaghi for 
valuable ‘guidance and recommendations during the site investigation 


and construction stages of the wor of the work and for his encouragement in the prep- 


aration of this paper. 
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dimensional compression ‘waves in dry , Ottawa. sand a are examined. The theory, 
which contains an explicit strain-rate effect, has previously been found to 
give good agreement with experiment. Results from it are compared with > 
elastic and plastic theories of wave propagation in soils. Problems associated 
with: the use and experimental determination of “dynamic” stress-strain 
curves are considered. Laws of similitude required by the present theory are ; 
Tentative data on the ‘response of of impact pressure gages are 


developed. 


__ Note.—Discussion open until January 1, 1962. _ Separate discussions should be sub- © 
mitted for the individual papers in this symposium. To extend the closing date one 
month, a written request must be filed with the Executive Secretary, ASCE. This paper — 
is part of the copyrighted Journal of the Soil Mechanics and Foundations Division, Pro- 

_ ceedings of the American Society of Civil Engineers, Vol. 87, No. SM 4, August, 1961. 


1 Engineer, RAND Corp., Santa Monica, —— BS 
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IMPACT WAVES IN SAND: IMPLICATIONS 
~ AN ELEMENTARY THEORY ~— 
Gama imnlinatinne af oan ploctin thane: fox the naman catinn 
‘letter’ symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, inthe Appendix, 
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ate detailed comparison between a theory of elastic-y plastic v wave e propagation 
in a strain-rate sensitive material and experimental results2,3,4 for impact — 

wave propagation in columns dry Ottawa and Fort Peck sands have been 

previously presented.° The fact that the theory has the striking ability to re- 

_ produce a number of important experimental trends with considerable pre-— 

‘cision would seem to give it es as a mathematical model for these 


Of course, there is a difference between Saying that a particular phenome- 

_ nological theory is credible, and saying that atheory for the dynamic response © 
of even a relatively simple soil such as dry sand has been established. Much 
experimental and theoretical work remains to be done before it is possible to | 
_ say that any theory of dynamic soil behavior has been established. However, | a 
the: very encouraging comparisons obtained with the drastically ‘simplified 

_ theory suggests that it may have some implications of practical value. For 
this reason results from the present theory willbe compared with some other 
_ theories of wave propagation in soils. Also, problems connected with the deri- 
vation of “dynamic” stress-strain curves from experimental stress histories — 
at fixed points in a sand will be examined. The laws of similitude implied by — 
the present theory will be examined and the present theory will be used to_ 
investigate the response of drum-type pressure gages to impact loading. 
~~ Because these topics are considered mainly within the framework of the | 

_ experiments examined previously. 5 The theoretical basis of that analysis will 
be summarized for - completeness. Thus the sand is treated as an elastic- 

last continuum which, for one-dimensional motions, can be described by 


w 


at eee 


(k[p-£e)],p 


The quantities t and x x denote the time and the Lagrange cniieate along the : 
_sand column respectively. The quantity €(x, t) denotes the compressive longi- 
tudinal strain, p(x, t) denotes the compressive longitudinal stress, and the 
stress — f(€) is a function which defines the static stress-strain curve. The - 
_ Young’s 8 s modulus E E is defined as the ‘Slope of the static stress-s strain curve 


as 33 2 «The Behavior of Soils Under Transient Loadings,” by R, V, Whitman, ‘Messen. ; 
< Internatl. Conf, on Soil Mechanics | and Foundation Engrg. Butterworths Scientific 
3 “The Behavior of Soils Under Dynamic Loadings, 3,” by RV. Whitman, et al., 
» Report on Laboratory Studies, M.I.T., Dept. of Civ. and San. Engrg. Soil Mechanics Lab., 
on ER 1954, AFSWP-118, Contract DA-49-129 Eng-227, Office of the Chf. Engr., 


D. Ww. Taylor and R. V. Whitman, M.I.T., Dept. of 
Civ. San. ‘Soil Mechanics Lab., July, 1953, AFSWP-117, Contract DA- 49-129 
Eng-227, Office of Chf. Engr., ASTIA No. AD-28- 
_ 53 “Impact Waves in Sand: Theory Compared with Experiments on Sand Columns, ° - 
by Blaine R. Parkin, ASCE, Vol. 87, No. SM3, June, 1961 
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The initial density of the sand is assumed to be given by the constant p. The 
strain rate constant is k. In the idealized material obeying Eqs. 1, 2 and 3 
disturbances propagate with the sound c =yE/. Eqs. 1, 2 and 
have been given previously by Malvern.9,f 
We shall apply the preceding equations to the "mechanical system 
schematically in Fig. In m denotes the n mass of the gage; is 


ane the sand seutiete velocity, at the reaction end, In this illustration the 
‘sand occupies the interval between x x=0 and x= =i, The stress gage at each end 
_ of the sand column consists of two lumped parameter elements; a mass m and 
a S epring of stiffness S. The alii is started from res rest by giving impulsively 


‘FIG, 1 _-IDEALIZ IZED EXPERIMENTAL AR ARR ANGEMENT 


Vo to the mass of. the right hand gage. Eqs. 1, 2 3 are 


subject. to the initial conditions 


in which is the initial impact The ‘conditions 
which apply on the right hand or the impact en end of the ies = on 


7 t) = 


6 “Plastic Wave Propagation in a Bar of Material Exhibiting a Strain-Rate Effect,” 
by L. E. Malvern, arterly of Applied Mathematics, Vol. 8, pp. 405,1951, | a 
_ 7 The Propagation of Longitudinal Waves of Plastic Deformation in a Bar of Mate- 
rial Exhibiting a Strain-Rate Effect,” by L., E. Malvern, Transactions, ASME, Vol. 18, 


1951, p. — 


| 
— 
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ig 7" and as a result o ie: 


At the left or end conditions are 


= 


Eqs. 1 through 9 determine a functions p (x, t), €(c, t), v(x, t) and V(t). In 
- some cases near the start of the motion, — function V(t) in Eqs. 5 and 6 will 
hot vary greatly from its initial value Vj. These two equations can then be 
approximated with the single dh condition 

It will be necessary herein to refer to experimental conditions discussed 

_previously.® Values of the parameters k and § will be used which are 
rived from the Ottawa sand data with» no | corrections being applied to them, 


= x 10 Ib per per cu in, 


This value of S corresponds to the impact gage cor responding 
_ value for the stiffness of the reaction gage used in the tests - uncertain but 
it appears that value 


= 18x x 103 Ib per cu ‘in, 


ates. ‘The: gage m had the dimensions of mass per 


unit area. In gravitational units the value of m obtained from the data given 
by R. V. Whitman, M. ASCE? is 


initial void ratio e of about 0.53 are ae "a Veen 
E= 40, 000 i 
= 00016 Ib sec? per i 
c= 15, 600 in, per sec 
COMPARISON OF T N OF THEORIES | 


In this sec section we we ‘shall not attempt to exhaust the ‘topic, For examy example if the — 
the present results with those obtained from 


with those of R. H. ‘Smith and Nathan M. Newmark, 8 F ASCE. Moreover, di- 

_ 8 Numerical Integrations for One-Dimensional Stress Waves, by R. H. Smith and 
ON. M. Newmark, Civ. Engrg. Structural Research Series, No. 162, A Technical Report _ 

“" to the Office of Naval Research Dept, of the Navy under Contract Nonr _ we q 
183, Dept. of Civ, Engrg., Univ. of Llinois, Urbana, fll., August, 


150 
’ 
— Me, 
— 
15 > 
— 
(| 
| 
— 
— 
— 
: = 
| | 


IMPACT WAVES 


rect calculations for an elastic- -plastic > material will not be: found in this 
we shall give some attention to the determination of “dynamic” 
‘Wave Propagation. - —In many concerning the dynamic 
response of soil it is customary to consider the soil as an elastic substance. f 
The elastic parameters for t the soil are derived from measurements of seis- 
mic wave speeds in the medium. Such a procedure might be useful for co- 
hesionless masses of sand as long as the amplitudes of the stress waves do — 
~ not differ greatly from zero. Otherwise it would seem that the epeeinnl of 
the e medium must be considered. 


elocity, v (in/sec) 


a. 


In support of this view results are compared from the present theory with 
those of wave propagation in an elastic medium having the same wave speed 
and density as the idealized Ottawa sand. It is imagined that 14-in. columns 

_ 3 both substances are subjected to identical impacts at V, = 60 in. per sec. 
impact-end and reaction-end stress gages of both have stiff- 
“nesses as specified previously. Similarly the strain- rate parameter k for the © 
idealized sand is that specified in ‘the ‘ “Introduction.” An example from these 
calculations is shown in Seine in which the curve labeled ‘ ‘plastic” rome 


. 
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“a It en been a the peak stress at the impact end of the 
column is significantly overestimated by elastic-wave calculations. As 
_ shown,5 the degree | of this overestimate depends upon the values assigned to 
the parameters k and S in the plastic theory. Of course this finding is ob- 


- tained from preliminary estimates of the behavior of the elastic column which 4 


surfaces. The highest compressive stress occurs for the present configuration 
at the reaction end of an elastic column as the result of the first reflection, 
- On the other hand, the present theory, which includes | plasticity, shows that 
: the | reaction-end stress is essentially limited by the static strength of the 
anil, even though | some strain rate and reflection effects may be present. 
These factors combine to produce ‘dramatic differences in reaction end re- 
"sponse of the two media in terms of stresses, particle accelerations, particle 
- velocities and particle displacements. It was | found that peak vi values of these © 
quantities differ by an order of magnitude. 
_ Finally it is observed that the striking - differences in the response ¢ of the 
plastic and the elastic media are in large part due to the method by which the © 
wave motions are excited in them, andthat these waves are longitudinal waves 
in bars. Changes in either of these factors would result in quantitative vari- 
_ ations in the response of both elastic and plastic materials, although in many 
situations the qualitative sense of the present results will be preserved. 


Simple Waves and Dynamic Stress -Strain Curves. this paper a possible 
constitutive relationship for a homogeneous granular medium has been 
_ studied which relates the properties of the substance under dynamic conditions 
those measured under ‘Static conditions. To find such a constitutive re- 
lationship which is valid for a1 any solid substance poses many difficult prob- 
; lems. It is rather remarkable that the simple law used here represents dry 
gand as well asitdoes, 
Of course, the present _ viewpoint is not the only one which might be con- 
7 "sidered for studies of wave propagation in homogeneous solids. Indeed, a 
_view which certainly antedates that of the present study concerns anyon 
in which the constitutive relationship, or stress-strain law, depends only on 


the instantaneous strain and not at all on the rate of straining. In such cases” 
one imagines that the stress-strain law appropriate to the medium ina a 


_ Namic situation may or may not resemble the static stress-strain curve of 


‘the substance. As already noted, the essential idea is that whatever the de- 
tailed nature of the dynamic stress- strain curve, it must be dependent only 
upon the strain. There is only an incidental interest in any relationship be- 
tween the static and the dynamic stress- -strain curves. The problem then, is 

to determine the form of this dynamic stress-strain curve from observations | 
of wave propagation in the medium. If the medium wees truly — 
independent for the range of strain rates of interest, the dynamic stress- 7 
strain curve would then form the basis for predictions of the response of a 
medium to various dynamic loads. There has never been doubt that two theo- 
retical models, one including strain- rate effects: and the other excluding yang 
factor, will generally lead to differing results. The ‘purpose here is to ex 
amine whether o: or. ‘not the additional concept of the dynamic stress- — 

- curve, when used with a theory which excludes strain rate effects, can _ 

satisfactory estimates concerning process of wave propagation in 
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z wave propagation in soils appears to have been G. I. Taylor.9 Taylor — 
‘to determine the plastic stress-strain relationship that the material must 
_ Obey in order that the velocity- time relationship observed at a fixed station 
within a plastic bar may be dynamically possible. Later C. W. -Lampson10° 
used similar ideas to determine dynamic stress-strain relationships from 
pressure observations at points in soil close to the point of high explosive © 
-detonations. Lampson’s results are derived from experiments in which spheri- 


cal or at least axial symmetry prevailed, and he used a heuristic procedure — : 
to correct his pressure-time data for the effects of spherical : spreading. ’ 


space of only one dimension defined by the longitudinal extent of the sand 


complication ii is not considered herein, Rather, like Taylor, a physical 


_ The possibility that one might be able tc todeduce dynamic strains and parti- 
cle velocities from stress (or strain) measurements at fixed stations ina 
soil is grounded in the theory of simple waves. 11 This theory is valid for 
Mea which exhibit no strain-rate sensitivity. (Strictly speaking the equations — 
of motion in the material must be “reducible. »11 For example, simple wave 
solutions cannot be found for spherically symmetric motions even if there _ 
&£ no strain-rate effect.) One manifestation that this theory might be applicable 
in a given situation is that travel-time curves of constant stress level in the 
medium should plot straight 1 lines. In ‘Spite of the fact that the present 


= contours of constant strain from ‘calculations resemble closely 
similar contours deduced from simple wave theory. However, Ottawa sand — 7 
appears to be far. more strain-rate sensitive than the metal considered by 

- Malvern, so that travel time curves of constant stress from the present theory ; 
a modest curvature, The particular example chosen for study corre- 
sponds an impact at t 115 in. per secina semi-infinite column of sand, 5 

Travel-time curves were plotted for a short interval following the impact. q 
_ The resulting plot does not resemble a simple wave solution until disturbances ~ 
— traveled at least 6 in. into the sample. At later times when stress re- 


laxation effects occur the resemblance is lost. Therefore for the present case 
it is plausible to use simple wave theory for the calculation of dynamic stress- au 
strain curves only during the first part of the loading when the stress is 
rising. Clearly, the results obtained depend upon the position in in the medium ~ 
as well as the intensity of the impact. 


g to experimental data in the following way. Suppose for example, that ae 

perimental travel-time “curves of constant stress ina strain-rate sensitive 
material have been obtained. Suppose further that this material is accurately 

> _ described by the present theoretical model. Can simple-wave analysis be 
— used to obtain at least approximate histories of the particle velocities = 


9 “Propagation of Earth Waves from an Explosion, ” by G. I. Taylor, Paper written 
_ for Ministry of Home Security, 1940. See The Scientific Papers of G. I. Taylor, Vol. I, 
“Mechanics of Solids,” Cambridge Univ. Press, 1958, pp. 456. = 
a 10 Final Report on Effects of Underground Explosions, by C. W. Lampson, NDRC Re- 
_ port No. A479, OSRD Report No. 6645, 1946. See Also Appendix B, pp. 410 in The Ef- 
fects of Atomic Weapons, U. S. Govt. t. Printing (Office and the Cc Combat Forces Press, 


11 Supersonic Flow and Shock Waves, by R. Courant and K. O. Friedrichs, Inter- 
science Inc., New York, 1945, Sec. 99, pp. 240. 
rt 
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strains in the e medium? Figs. Figs. 3 and 4 dynamic stress- ow curves and _ 
a particle velocity versus time curves for two fixed stations in the sand are- 

‘presented in answer to this question. ‘These positions are located at x - =2 in. 

“a and at x = 12 in. from the impact end of the sample. At both stations : it will 

ber noted that, at a given st stress, the strain calculated from the ) simple wave 


— Present. theory (strain rate effect included) 
bs, 


im act velocity,Vo IS in/sec 


“10 12 14 «16 18 
Longitudinal strain (percent) 


“Fic. _—-DYNAMIC STRESS-STRAIN C CURVES OBTA OBTAINED ‘FROM 
PRESENT THEORY COMPARED WITH DYNAMIC STRESS- pin 
CURVES OBTAINED BY USE OF SIMPLE-WAVE 


_ assumption exceeds the str strain elias from _ present theory by a factor | 
of approximately two. ‘Similarly, at a given instant of time, and therefore at 
a given stress, the particle velocity obtained from simple wave computations 


exceeds that obtained from the present theory by nearly 50% at x = 12 in _ 


= 2 in, thi in particle vel 
When x = 2 n, s excess particle cle ve locity amounts to about 20%. Evidently — _ 
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war 
simple wave ‘approximation | is chiefly of qualitative value in soil dynamics 
problems where appreciable strain-rate effects are involved, On the other | 
comparison of these results with the elastic calculations /examined 

_ previously indicates that the plastic behavior of the soil is an essential in- 
-gredient which, even more than strain-rate be for 
_ discussions of intense wave propagation in soils. 


_ Dynamic tests on sand bring into prominence certain properties of the 
) which are not of crucial importance in static tests, and the question 
arises as to the proper method of extrapolating results from dynamic _ 
: tests in the laboratory in order to predict the dynamic behavior of a sand > _ 
formation in the field, or as to the proper method of correlating existing 
¥ oratory results. Of course, the information now available does not permit ~ 
to this question. However, the success obtained with the 
/ Bresent theory leads to the hope that that its use | to indicate a a ‘partial answer wil answer will an 
The aims and methods of the theory of similitude certainly do not ae 
4 illucidation here. It is sufficient to note that characteristic physical quantities | 
' for the present idealized substance are defined by the parameters p, E, k and Pt 
te). If we call the characteristic length 1, we find from the usual sort of = 
_ analysis that we may relate the characteristic time to this length by means of _ 
the characteristic wave speed c given in the Introduction. In particular it is | 
| found that if quantities (kl/c) andf(e)/Ehave the same values in the laboratory 
as they do in the field, then scaling laws of the sort proposed by Lampson, 10 
- _ or in a slightly more general form by Parkin, 12 are valid. (G. Birkhoff13 has | 
‘grouped such scaling laws under the heading of Mach similitude.) The first : 
these two quantities (—= ye is similar in form to the quantity | which must be i. 
_held constant in order that these scaling laws “may be used for linear viseo- 
e elastic substances. This result could have been foreseen from the existence - 
of a certain likeness of the differential Eqs. 1, 2 and 3 to the equations of — 


depends upon the | value f(e€)/E. And as ; also found for the elastic- plastic ma- “— 
— terial,12 this requirement indicates that motions in the two materials will be ~ 
similar only if every 0 ordinate: of their stress-strain curves are properly 


dition is obvious. General knowledge of the precision with which it must be 
_ ‘satisfied would be desirable. From the particular cases studied here it was _ 
i - found that for some purposes it may be possible to relax this requirement 
7 _ whereas for others it must be rigidly satisfied. In this connection it is note- 

that R. L. Bjork!4 has found that a condition similar to the writer’s 


- 12 A Review of Similitude Theory in Ground Shock Problems, by B. R. Parkin, The 
Corp., Santa Monica, Calif., Research Memorandum RM-2173, April, 1958 (ASTIA 


_ 13 Hydrodynamics, by G. Birkhoff, Princeton Univ. Press, 1950. 
14 Numerical Solutions of the Axiall 5 mmetric Hypervelocity Impact Process, by 


every strain within the range of interest. The formal stringency of this con- J % 


— 

g 
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4 plasticity requirement ‘must be satisfied in order that experimental results. 
from hypervelocity impact of metals may be properly scaled with impact 
velocity. In this case the pressures involved are so high that the metals be- 
come fluid and it becomes necessary to speak of an equation of state rather 7 
than a stress-strain curve. In the present case the existence of two para- — 
meters which must be ‘Specified, in order that the customary scaling laws - 7 
may be validly applied, places rather severe restrictions upon the scaling 
of laboratory results, This difficult situation may be severely aggravated if 
- the nature of the external agency which produces the motion in the medium 
must _be taken into account. For example, if one wished to scale the 


=— as = and | = are kept invari- | 


 LINEARIZED ‘THEORY 

In order to linearize the equations of motion, Eqs. 1, 2 and 3, it is neces- 7 
; sary only to replace the nonlinear static stress-strain relationship in Eq.3 | 
by a linear function. In terms of the dimensionless ‘stress such 


Al (e) = - 

in which 6 determines. the of the line and (E give the ‘coordi- 

nates of an arbitrarily chosen point upon it, For example if Eg = Oo= 0 and 

6 = 1 the medium will be taken to be elastic, whereas if 6 = 0 for some fe) 
- and for some range in €, the substance is perfectly plastic within this all 


order to: systematize the calculations it is helpful to introduce dimen-_ 
sionless variables similar to those used by Malvern. 6 Thus, 


xk 


.(12b) 


> 
and 


‘ Introducing oy new variables and Ea. 11 into the equations of motion oldie: 7 


‘4 
7 
that the additional parameters U : 
tn the transformation, althou De seen itt 
~ ph : of the linearized theory, there are some conditions of limited applicability in a _ 
og 4 which the last of these added constraints may be ignored. - 
4 7 | 
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of conditions, However, the present examination will be restricted to con- 7 * 
sideration of the first part of the motion of particles near the impact end of a_ 

sand column, In this case the simplified boundary | condition of Eq. 10 may be —— 
used on on the surface x = 0. Since no reflections at the surface x = 1 will be 7 
considered conditions at that point may be replaced by the condition that all 4 
_ stresses, strains and particle velocities must vanish at large distances from q 

0. Thus, the modified boundary conditions may be written in terms of 


to the imp act. -veloct Vv 
impact Vo and the 


ont 


(y, 0) = 0) = 


25 


: mated by a straight secant. " Gaveiens €o and go are made equal to zero <a 
Eqs. 11 and 13. The Laplace transform of r), uly, or of 

_ denoted by superscribing a bar over the appropriate symbol as follows ~~ 

> a s is not to be- confused with | with the displacement s(x, t) previously — 


For an elastic solid 41, the linearized becomes exact, and Eq. 


Maw 


T 


he inverse of this transform is known to bel5 © 


_ 15 Modern Operational Mathematics in Engineering, by R. V. Churchill, McGraw- | 


Inc., New York, 1944, 


may _ a re 12 norreecnnn no wa? 
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1 Eq. 20 is expressed in terms of real time and space instead of the dimen- 
_sionless variables T and y it willbe seen that this result is independent of the 

strain-rate parameter, _ k. This well-known relationship has been rederived | 
: here because it will be needed later and, because it emphasizes the fact that 


when. a spring ‘suddenly pushes against the endof an elastic rod the « elementary 


Impact Gage Response in ina Strain -Rate Sensitive the 


primary interest is in events which occur at the impact surface y = 0, we 


‘Shall - y = 0 in Eq. 18, In this special case the transform becomes ~ 


and its inverse may be found with the aid of a convolution, The precise nature 
of the inverse depends upon the order and the affix of the | poles of (0, : s) in 7 
the complex s-plane. Simple poles which are located on the real axis in the | 
Be -plane lead to real exponential functions in the inverse, and poles which are 7 
located off this axis imply an inverse transform involving trigonometric 
functions. In order to determine when we will have one or the other of these 


may be that if if the = 


is all roots = 0 will be real. As it happens, q(s) = 0 has one 
root which is real and positive for all @ in 0< 6 < 1, but if the condition of 
; Eq. 23 is not satisfied the remaining roots will be complex. If only the ‘equality 
in Eq. 23 is considered this equation defines a locus of points ina 6 -a plane 
which divides the plane into two r regions, In one area the solution o(0, t) will 
involve circular functions. In the egion of the @ - a plane exterior to this 
: o(0, t) will involve hyperbolic functions . Fig. 5 shows that portion of the a 
- a plane which is of physical interest. ‘Points inside the curve correspond 
to solutions involving complex roots so that an initial spike in the impact 
stress records is to be expected. If a point (a, @) lies on the curve or inthe | 
area outside of it, the fact that the linearized solution does not contain com- - 
.~ exponentials is not sufficient to exclude an initial spike from the impact | 
stress record. Whether or not the linearized result predicts a weagenereell 


stress rise without an initial peak depends upon the precise location of the — 
particular point (a, 6) with» respect to the locus. In passing it is noted that 
- aecording to Eq. 20 all points on the line 6 = 1, that correspond to an elastic © 
‘medium, are not in a domain of the a - @ plane where initial stress spikes — 
occur, When the point (a, 6) does not lie on this line the question of whether 
or not a stress spike will occur is not simply resolved. aoe ee natne ae 
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; yy y ~ stress immediately after the impact to a maximum limiting value. When Up pee 
2 is constant as assumed, an initial stress spike or peak does not occur. The _ 
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‘In order to clarify this point and to test the accuracy of the eeeteed ; 
_ theory the inverse transform o(0,7T) was determined and the resulting form- 
“ulas were evaluated numerically for values of a and 6 corresponding to the 
points A, B and C of Fig. 5. The numerical results from the linearized theory 


im then compared to corresponding results from the nonlinearized theory. 4 

The linearized theory was found to give results which were in satisfactory — 
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FIG, 5. GAGE | DIAGRAM 


A was computed from known values of E and c for Ottawa sand together with _ 
gage stiffness and the strain-rate ‘parameter k from the introduction. 
‘The selection of a ‘proper value of @ is somewhat bothersome because we 
' must know in advance the range inthe strain, €(0, T), which will occur shortly 
after impact. For a given” ‘material the amo amount of strain will depend largely © , 


| 4 
4 | 
. 
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In the numerical evaluation of the linearized theory the value of a at point 4 aS 


4 


P 
‘upon the impact velocity This ‘property is used the present study by 
3 selecting the slope of the secant to the static stress-strain curve 6E, so 
this line intersects the static curve f(e) at € = 0 and at that value of € corre- 
"sponding to the peak stress immediately after an impact. For Ottawa sand — 
the necessary data may be taken from the calculations that lead to Fig. 8 of — 
the writer’s previous work. 5 The results for € and f(€) give values of @ for a 
range of impact velocities Vo as shown in Fig. 6. Thus it is seen that the 
point A of Fig. 5 falls within the experimental range for Ottawa sand. Fig. 6, 
although constructed in a rather arbitrary manner, is useful because it — 
cates, at least qualitatively, how one might account for ‘ the value of Vo: in his © q 
choice of @. It seems likely that for impact velocities that are significantly 


_ greater than those considered here, it woul would ld introduce 1 little error ‘if if one were 


velocity 


; EFFECT O OF IMPACT VELOCITY UPON | CHOICE OF LINEAR 


carry out the present. development o of the linearized theory with (This 
additional simplification has been found useful for the calculations of the initial ” 
response of the impact gage to a second impact. 
_ The three examples from the linearized theory at points A, B and C in Fig. 

and the result given by Eq. 20 for the elastic medium permit a full interpre-_ 
tation of Fig. 5. It seems clear from both the linearized and the nonlinear = 
theories that if a and 6 are chosen so that the point in the a - @ plane lies” 

« @ below the lower branch of the curve of Fig. 5 then the impact stress-time a 4 
7 16 Impact Wave Propagation in Columns of Sand, by B. R. Parkin, The RAND Corp. ; 
Santa Monica, Research Memorandum RM-2486, November, 1959, ASTIA 

AD234961. 
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_ record should show no peak. If the (a, 6) is inside this locus. the first 
a part of the stress-time record should show a peak. Moreover, for this case _ 
the linearized theory indicates that the solution is oscillatory in form and that _ 
these oscillations are damped. Finally, when the point (a, @) is above the 
upper branch of the curve a stress peak may be expected in the impact records. - 
os a This peak is not associated with oscillatory solutions; and as 6—~1, that is, = 


as the medium approaches perfect elasticity, the stress spike should vanish, 


The consequences of choosing alternative theoretical assumptions regarding 
_the propagation of longitudinal stress waves in dry Ottawa Sand have been 
‘compared. For example, suppose calculations are matched to the experimental — 


Ls seismic wave speed but the study of large amplitude wave propagation in the 2 | 


-_ medium are based upon the theory of elasticity. Then compared to the — 
theory, or any other theory that accounts for the plasticity in the medium, — 


- elastic calculations can lead to significant overestimates of stresses, particle 
displacement, particle velocities and particle accelerations in the onl 
_ On the other hand one might wish to include the assumption of plasticity ~_ 
‘his calculations but use a “dynamic” stress-strain curve in order to avoid 
taking account of strain-rate effects. If strain-rate effects are important then, e. 
depending upon the magnitude of the strain-rate parameter, it appears a 
this procedure may lead to results which are no more accurate than would — 
be obtained by a straightforward neglect of all strain- rate effects = iin oe 
_ In the present simplified theory the occurrence of the impact stress peak 
of the experiments is ascribed entirely to an interaction of the — 


1G derive design criteria for impact stress gage response under a variety “4 


conditions, The possible usefulness of these results in experimental assess- 


ii ments of the strain- of soils may be nates. 


| 


"The preceding results have been made available through “ne efforts o 
Mrs. Susan Belcher who carried out the numerical work, of Mrs. onal 


Lopez who prepared the manuscript, _andofG. I . Margadonna and his staff who 


APPENDIX. -NOTATION 


| The terms used in this ener are listed here for convenience of of nie 


sensitivity and static stress-strain properties of the medium with the » elastic” 
‘as of the impacting stress gage. On this basis it has been seamed 


 &§ 

: 

tm 
3 
as 
= 


strain- rate parameter; 


mass of gage per wit areas 


= gage stiffness (pressure/ deflection) ; 


mensionless particle velocity; 

= = dimensionless impact gage velocity; _ 


= di mensionless initial impact gage 


= particle velocity; 


= gage 


= compressive longitudinal 


= initial density of sand; > 


= dimensionless stress; 3 and 
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DYNAMIC TESTING OF PAVEMENTS?» 


to by W. and C. R. Foster 


the aaa on page 153, inthe April 1961 Journal of the Soil Mechan- 
ies and Foundations Division (Proc. Paper 2804) should be deleted and replaced 


W. HEUKELOM®5 and C, R. “FOSTER, 36 F, ASCE.—The final paragraph of 
by the following: | 
3 The numerical evaluation of the damping force, , considered by ALA A. Eremin, 


Mz. ASCE, is not based on Eq. 15 but on a combination of Eqs. 15 ‘and 16 fror 


which it can | be be determined tl | 


The value of A is determined by measuring S and ¢ at a known frequency w. 
T 


he writers are grateful toO. M. Christof Ehrlerfor his remarks and coe | 4 

tions. _ The relationship between the elastic stiffness R and the dynamic deflec- 

of the medium is 3s basically given in Eq. 8. A 
A a 
ac 
in which Gp isthe shear modulus of the unten p represents the factor which - 


depends on the distribution of the load over the loaded area with radius a a: 


= 
I-p 
3 according to J. ‘Boussinesq37 for 


ibut ion, 


Ras 
February, 1960, by W, Heukelom and C. R. Foster (Proc, Paper 2368), 
35 Konin Klijkel Shell- ‘Laboratorium, _Ameterdam (N. Vv. Der Petroleum 
36 Coordinator of Research, Natl. Bituminous Concrete Assn., Texas Agric. and 4 
- Mech, College, College Station, Tex, Formerly, Asst. Chf,, Soils Div., Corps of Engrs., a 
U.S. Army Engrs, Waterways Experiment Sta., Vicksburg, Miss, __ 
37 “Application des Potentiels,” by J, Boussinesq, Paris, 1885, a? 2 
38 “Druckverteilung im by O. Vienna, 1934, 
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For the s gives the best fit. By comparing tl the ob- 


_ served values of R and Gm (from v wave velocities), the following values of the ~ 


Poisson ratio were derived according to Frohlich: = 


Soft wet a 0. 42. 


Soft loam = 0. 39 


‘is 0.4 rather than the value 0.5 which is used in cases, 


| "eae of the R-value during vibration have been observed, es: especially 

: on clay soils and have, so far, been attributed to thixotropy. In the same way, 

the amplitude of the vibrations affectsthe stiffness, but this was found to affect 

7, ¢ primarily the observed value of M and to a much lesser extent the value of R. 


far as the observations enable conclusions to be drawn, R can be 


to. 5 ‘enn: for bare soils, 


= tons flexible ‘pavements, ont 


The mans ‘of the ection. machine is 2 the of the soil- mass 

was considered by H. 
The value of 9 has been found to be small in the range of practical mates 


times t.. _ For various reasons the actual loading time in — vibrations 


8 ‘cps. It is doubtful whether resonance at higher frequencies is real. 7 In the — 
opinion of ith vi 
observed mass, which can also give rise to a maximum inthe amplitude- -versus- 
frequency curve, thus true resonance where, the 


would be much greater. 


39 “Fourier Transforms, by I. N. Sneddon, London, 1951, 
40 Grundbau-Dynamik, by I H. Lorenz, Springer, 1960, 
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E, RICHART, 13 F, ASCE. —The writer wishes to thank all the dis- 
cussers for their generous contributions of ideas and information on the ge gen- fs 
- eral subject of vibrations of foundations. It is evident that much work remains 
to be done on the problem of developing satisfactory design procedures = 
foundations subjected to impulsive or repeated loads. 
_ Coates feels that the repeated static load type of test has practical i 
_ bilities for evaluating the behavior of soils under dynamic loads. The writer 
agrees t that a series of plate loading tests at several depths should give | useful 
data, provided that the range of deformations corresponds to that encountered 
under loading conditions. order to obtain reliable load- 
~ must te used that are several orders o: of ‘magnitude more sensitiv e than those 
_ customarily used for plate bearing tests. This increased difficulty of meas- 
urement also requires § skilled test personnel. Consequently, tt the choice of meth-— 
_ odused for local evaluations of a building site, by seismic methods, by steady - 
7 state vibration methods, or by repeated static loadings depends on the budget | 
_availablefor site exploration. If the static method is to be used as the primary 
_ exploration procedure, the writer wouldstill recommend that one of the dynamic 
-methods be usedasacheck, 
= With regard to the problem of starting or stopping, the — problem — 
tobe resolved is the probable maximum amplitude of vibration to be encountered 
7 asthe resonant frequency is passed. If the system is highly damped (low value | 
ot b) a low magnification of amplitude occurs over a rather wide rangeof fre- 
- quency, but intolerable amplitudes of vibration should not occur, The amplitude 
of oscillation under the steady state condition represents the maximum value 
7 to be developed. If this theoretical maximum amplitude is larger thancan be 
tolerated indefinitely, an estimate can be made of the largest amplitude to be 
developed during starting and stopping by the graphical phase- -plane method 


drop hammers can be estimated by use eo Fig. 8 if the soil parameters are 
carefully evaluated for the entire range of stress involved, ree eee 
Ih compliance withthe suggestion by Coates and Eremin that Fig. 8 be pre- 
sented with E as a parameter instead of G, Fig. 20 has been prepared z as an al- 
_ ternate: for Fig. 8(a). . Byt using g these coordinates the curves for the constant 
- force oscillator nearly coincide, and those for the rotating mass oscillator ap- 


pear t to move. closer oser together. ‘However, the use of E instead of G as a a para- 


4 August, 1960, by F. E, Richart, Jr, (Proc, Paper 2564), = 
Prof, of Civ. Engrg., Univ. of Florida, Gainesville, Fla, 
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= meter ordinarily will not eliminate the need for ¢ tinting an effective Poisson’ S- 
4 ratio of the soil. ‘The shear wave in a saturated soil is a much more reliable | a 
- measure of the stiffness of the soil structure than is the compression wave, a 
~ because the measured value of the latter | is influenced by the presence of water. 7 


to be evaluated. = 
ae _ The use of E as a parameter in curves of amplitude factor versus frequency 

actually causes increased divergence of the curves that can be prepared from — 

the information shown on Fig. 8(b). Consequently, this information is not shown 4 
38 Fig. 20. It should be noted that Fig. 8(b) is plotted using a different pair _ 

of parameters as ordinate and abscissa than are used in Fig. 9(a) and Fig. 


10(a). _ Because of this arrangement, a numerical value of amplitude of vibra- 
tion forthe rotating mass oscillator can be obti obtained without hout requiring any e eval- a 


FOR ROTATING MASS OSCLLATOR 


FIG, 20. ALTERNATE PLOT OF MASS RATIO VERSUS. FREQUENCY» 
s,s FACTOR RELATIONS FOR VERTICAL OSCILLATION 


Fe 


by Whitman, M. ASCE, of MIT. 
a Coates and Perri are correct in pointing out that the radius of gyration of 


the rectangular prism for an axis through a center line in the base should be 


we + a2), ‘This ot changes esthe computed amplitude a and frequency i in Example oo 
0, 000176 radius and 218 rpm, respectively, 
Alan has raised objections tothe possibility of resonance in a semi-infinite, 
- homogeneous, isotropic elastic body anduses Polz (39) and Christopherson (40) 
as referencesto support his contention, Alpan has not carefully identified his 
_referencedinformation with the problems to which they apply. Inthe first pa para- 
: = graph of the section under Vertical Impulses Acting at the Surface, the writer 
= noted that Lamb (17) didnot find resonance resulting from a single concen- 
This is load type of impulse which acted on the semi-infinite elastic medium. 
'‘This is the basis for Christopherson’s criticism of the method for determining © 


\ 
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SM 4 ‘DISCUSSION 
_ the soil constants by a single impulsive loadas proposed by Crockett and Ham- 


mond, Polz also_ noted that resonance was not possible for the condition of a 
concentrated pulsating load, which is again consistent with Lamb’s findings. © 
_ Resonant conditions may be obtained in semi- infinite, homogeneous, iso- 7 
tropic, elastic solids whenthe periodic loads are distributed over a finite dl 
tact area at the surface of the solid. The study by Polz (39) also showed this -_ 
and by use of the curves given on Table 1 of his paper withthe data for Exam- _ 
ple A, the ‘computed amplitude of oscillation is very close to that computed 
from Fig. 8(b). (It should be noted that Polz considers onlythat caseofacon- _ 
stant force type oscillator). It appears that Alpan has not considered the radius 
= the contact area, rp, which appears inthe expressionfor resonant frequency 
(Eq. 11) to represent the characteristic distance required to establisha natu- — 
frequency. The expression for frequency (Eq. 11) obtained for the circular 
contact area and the expression tor frequency obtained by Polz (39) which con- 
tains a a length parameter, h, are all dimensionally c consistent with Eq. (21) as — 
well as with the expression for frequency obtained by Lundgren (65) on the 


_basisofdimensionalanalysis, 
4 Alpan notes that the methods discussed inthe paper didnot include the prob- 
lem of damping. The mass ratio, b, represents a measure of the dissapation 

of energy by elastic waves in the elastic body, but no internal damping of the 
soil is considered. In most cases, the amplitude of oscillation in the operating 
frequency range is the design criteria. Consequently, the internal damping in 
the soil will act to reduce the amplitude below that computed from the elastic 7 
theory, and, as reliable information becomes available onthe damping ‘capacity 


@ soils, _ this secondary correction can be made if desired. 


of vibration of a soil- system, The amplitude obtained by 
i= method depends on the damping ratio chosen, and typical values have re been — 
accumulated for years for different test installations on different soils, 
a This value of damping ratio computed from test results lumps together the 
7 results of the internal damping of the soil and the dissipation of energy —— 


dar - damping r: ratio for conditions where the geometry and masses involved r may 
be appreciably different from the test conditions. This is where the elastic 

_ half space theory has an advantage, because the results are presentedin a man-_ 
: ner that permits extrapolation of model test informationto prototype conditions. 

_ Relatively small mechanical oscillators have been used extensively to pro-_ 

wide test information on the dynamic behavior of soils. Many of these tests 

_ have indicated a non-linear behavior of the soil under increasing load ampli- 
"tes . The writer pointed out that the amplitudes of motion developed by the 
: oscillating mass often were large enough that the computed acceleration be- 
came greater than 1.0 g. This indicates that the oscillator jumped free of the 

| round during the upward motion. Asa result, the contact pressure at the os-_ 
cillator base varied between zero and a maximum during each oscillation, and 
_ the “hysteresis modulus of elasticity” had a lower slope than it would have had | 
_ if the minimum pressure had been higher (Fig. 21(a)). This effect of range es : 
stress on the effective modulus of deformation, or hysteresis modulus, was” 
pointed out by Terzaghi | (66) for repeated static loads on sands, and by“ Tschebo- 
= (discussion of reference 33) for dynamic loads on foundations. Because 
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August, 1961 
this effect, in the conclusions to the paper the author recommended that os- 
cillator tests utilize the same vibration parameters (which in this discussion 
refers to range of stress) as are to be usedinthe prototype. 
‘The decrease in hysteresis modulus with increasing range of stress corre- 


= to the idea of a soft spring characteristic as discussed by Lorenz (45). . 
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\STATIC 


to 


THEORETICAL | PRESSURE 


(b) P PRESSURE ON CIRCULAR BASE 
FIG. aa, |.—FACTORS CONTRIBUTING TO NON-LINEAR 


must be. a in the pressure dis distribution at the footing base, Fig. 
21(b) illustrates one concept of the change in ‘maximum pr pressure distribution 
4 from curve A under static load alone, to curve B for static plus smalldynamic | 
loads, to curve C for larger dynamic loads. The maximum pressure cannot 
exceed some which is by the curve 


— 
wi 


M4 


changes from one wherethe majority of the loadis carried near the perimeter 
5 of the footing to one where the pressure is concentratednearer the center, the _ 
resonant frequency will be reduced (Sung, (3) and Richart (24)) iad 
‘Finally, Alpan suggests that because the amplitude of oscillation computed _ 
from the elastic half-space theory | was appreciably less thanthat he computed 
_ fromthe damped mass -point theory using an assumed typical value of the damp- - 
ing c constant, that the former method was unsafe. The applicability of any theory 
for design” purposes must be proved by field data obtained from operation of. - 
‘ the actual installation. There is stilla scarcity of reliable published informa- — 
ff tion concerning the dynamic behavior of large footings under periodic loading. 4] 
_ However, at the present time . (1961) a comprehensive series of tests is being | 
conducted at the United States Army Engineer Waterways Experiment Station, © 
for the Military Construction Division of the Office, Chief of © 
Engineers, Washington, D. Vibratory tests are being carried out on four’ 
circular reinforced concrete bases of 62 in., 88 in., 108 in., and 124 in. indi- 
ameter using vertical, torsional, and rocking modes of oscillations, It is the _ 
writer’s opinionthat the resultsfor vertical amplitude of oscillationfrom about — 
4 the first forty tests agree well with the values computed using Fig. 8(b). 
Perri has pointed out the expressions that may be used to evaluate ampli- 
“tudes of oscillation and resonant frequencies of the coupled sliding and rocking = 
- motions for an undamped mass subjected to a periodic horizontal force. For | a 
a sucha system, the amplitudes of motion become infinite at the respective values ’ ; 
of coupledr resonant frequencies, Consequently, this particulartype of analysis 
is of little use in evaluating the maximum anticipated amplitudes of oscillation 


resonance, 
However, the wri writer is grateful | to Perri f for emphasizing the very incom- 7 


> 4 on Figs. 9 and 10 represent the parameters involved in uncoupled atiding | and a 
rocking that were > obtained mathematically be Byeroft, and 


_ the oscillator contact zone and introduced into the mathematical solution, Physi- 
- cally, this requires that the center of gravity of the oscillating mass must coin- 
cide with the line of action of the resisting elastic forces that, for the horizon- 
: tal sliding case, requires that the mass be concentratedin a disk of zero thick- | 
ness. In order to get some idea of the order of magnitude of the vibrations, 
: <4 the writer suggested that the amplitudes: computed from the two uncoupled modes 
qt of vibration be superposed. This procedure is correct from the standpoint of 
, mechanics, but may not represent more than a crude approximation to the ac- 
4 tual dynamic behavior of afoundation, The differences between the frequencies 


; ‘tation, or rotation without translation, the boundary conditions < are ‘set up for ’ 


determined for the uncoupled modes and the coupled frequencies may be esti- 

mated from the curves given by Jacobsen and Ayre (64) p. 337, or the lower 

of the two couple frequencies may be evaluated by the Southwell-Dunkerly pro- 


It should be noted that t there are mistakes in the equations. presented by 
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oF Fig. 12 was included for use in estimating the maximum stiffness that could — 
we obtained bythe use of point bearing piles that rest on an infinitely rigid stra- 
tum, This is an weer limit that can be reached only ‘theoretically. Actually, © 
the pile tip will move, , and Eremin notes that the effect of surrounding piles “e 
- should help stiffen the bearing layer that would, in turn, reduce the pile tip mo- 


| 
> (5) -[n+n-5+1 n= (39) 
n{1+ +1] 4 n[{1+ +1| -4n -(40a) 
= 9, | +3 eal + 
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—* has presented some interesting data onthe relations between confin- 
‘ing pressure and the shear modulus of Ottawa sand, as influenced by changes © 
in void ratio. The effect of change in void ratio is quite significant. He also 
observed that the modulus varied with about the 1/2 power of the confining pr pres- f 
‘sure, which means that for this material the shear wave velocity varied with . 
about the 1/4 power of the confining pressure. This isquite reasonable inview | 
7 of the deviations from the theoretical slope of 1/6: found by Duffy and Mindlin 
(12) for longitudinal wave propagation in bars made upof uniform steel spheres. | 
In the range between 1000 psf and 2000 psf, and for the first mode of vibration — 
with the (1,0,0) orientation, their test data indicates that the longitudinal wave 
velocity varies with confining pressure according. ‘to the 1/5 power for neh 
tolerance balls (1/8 + 0.000010 in, diameter) and the 1/3.9 power for low tol- 
erance balls (1/8 + 0. 000050 in, diameter), . Actually it is remarkable that 
a the test results for Ottawa sand correspond s so well to those obtained from a 
7 granular medium composed of but slightly imperfect spheres. | 
_ Smoots, Stickel, and Fischer have brought up several interesting points. Of - 
course, the use of a homogeneous, isotropic, elastic medium to represent a 
mass of soil represents a drastic simplification, and the approximations in- 
volved in this substitution must be re-evaluated for each different use of the - 
theoretical results. Inc connection with the vibrations of foundations, the elas- . 7 


‘tic waves representa non-destructuve method of evaluating the stiffness ofthe 
soil in and adjacent toa constuction site, 
They mention the observed effects | of ground water onthe compression 
note that the shear wave seems to be little affected bythe water table. ‘The 
_ theoretical studies by Zwikker and Kosten (67) and Biot (68) onthe elastic __ 
"wave propagation in saturated porous materials have shown the influence of aq 
saturation on the shear and compression waves, The shear wave is ; uncoupled — 
from the compression waves and is influenced only by the added mass of the 
fluid, whereas the compression waves in the fluid and in the soil structure — 
are coupledand dependon boththe stiffness and mass of each component. Har- 
din (69) has combined the characteristics of a granular medium withthe theo- — 
ries of Zwikker and Kosten and of Biot and found that for saturated materials © 
the shear wave velocity varied with the 1/6 power of the confining pressure as 
it does for the granular materialin the dry condition. Consequently, the velo- 
city of the shear wave as determined from seismic tests, or from steady state 
vibrations, should permit « evaluations of the shear modulus of the soil that can — 
be used in the theories for the dynamic behavior of foundations. = 
-Smoots, Stickel, and Fischer have made computations for the resonant fre- 
quency of the oscillator described i in Example Aby several methods, © Tschebo- — 
tarioff’s method utilizes the single degree of freedom t undamped mass theory 
and a typical value for sandy clay soils. No determination of the amplitude of 
oscillation can be made by this method. Toriumi’ Ss approach tothe problem is_ 
_ essentially the s same as Reissner’s (2) in that he assumes a uniform distribu- 
= of pressure over the circular loading area and a constant force type of 
~ oscillator. - Apparently Toriumi has an error in his evaluation of the function | 
f2, that would affect the answers obtained from his curves. Because Toriumi’s 
4 method is based on the assumption of a constant force oscillator instead of a _ 
rotating mass oscillator as used in Example A, it es al myemes that his 
The test conducted at ‘the Kwajalein ‘site field data that 


a is needed for a better evaluation of the theoretical procedures. Use of a shear 
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August, 1961 
modulus of 1075 psi as obtained by Pauw’s method would, of | course, givea 


_ “computed resonant frequency that was lower thanthe test value. A considera- 


er of the curves on Fig. 16 and the average static contact pressure atthe base 


AMPLITUDE FACTOR = g 


(a) AMPLITUDE FACTOR VS FREQUENCY FACTOR RELATIONS 7 


(b) MASS RATIO. FREQUENCY FACTOR RELATIONS 


‘FIG. 22. .—CH ARACTERISTICS OF TORSIONAL OSCILLATION 
FOR AN OSCILLATOR RESTING ON A 


INFINITE ELASTIC MEDIUM 


psi in the region directly below the base. Consequently, the oemetiemee 
mene onthe basis of this shear modulus is unduly low. The measured am-_ 
of 0.0022 in, is lower than that given by use of 8(b). 
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DISCUSSION 
the results — just one test are not sufficient to evaluate the effectiveness 
: The motion of the radar dish can also impart torsional oscillations to ‘i ; 

_— that may develop resonant vibrations. - Solutions for this problem 
have been developed by Reissner (70) for the case of shear stresses cule 
linearly with radius, and by Reissner andSagoic (71) for a rigid circular base. 
Fig. 22 22 gives the curves from which the maximum amplitude and resonant fre- 
quency: canbe computed, ‘Because the torsional oscillation is an uncoupled mode 


of vibration, these motions can be considered 


B study, moth analytical and experimental, should be conducted on each of these 
topics, as well as Or on the subject of Sastrumentation for this particular rai range o of 
_ Wilson and Miller have contributed valuable te test data cc concerning the effect 
_ _ of confining pressure on the compression wave velocity, and on the compres- 
_ sion modulus andshear modulus of elasticity for several granular soils. 'Gen- 
erally, the compression modulus varies with about the 1/3 power of the con-_ 
_ fining pressure. The two materials that deviated markedly from this behavior 
designated (SC) clayey gravelly sand and (SM) silty gravelly sand, 
‘The resonant column procedure for determining the elastic wave velocities 
in soils was used by Hardin and by Wilson and Miller. In this test, a column 7 
of “soil is subjected toa confining pressure then set into resonant reased,the 


: to determine the wave velocity. As the confining pressure is increased, the 
void” ratio of the material decreases in accordance with its compressibility. 7 


of confining { pressures shown on Figs. 6 or r 18, but it is of importance for softer a : 
‘materials . Hardin’s curves indicate thata considerable increase in slope 


ment of confining pressure to result in the corresponding change in effective 
> _ stress in the soil. If the dynamic measurements are taken before consolida- 
tion is complete, the wave velocity corresponds to some pressure which is | 7 
lower than the confining pressure, and a flatter slope will be produced onthe 
- confining pressure versus wave velocity diagram. Tests carried out by J. RO 
‘Hall, dt. As M. ASCE, in the writer’s laboratory have shown this effect for 7 
_ Vicksburg loess. For a confining pressure increase from 360 psf to 1320 psf, 7 
- the computed shear modulus increased from 4100 psi to 5470 psi after 20 min, 7 
and to 6670 psi after 20 hr. corresponding increases inthe exponent c corre- 
- sponding to the relation in Eq. 44 was from 0,22 to 0.37 as a consequence of 
the additional time the confining pressure was applied. There also appears to. 
be > some influence of the maximum pressure applied, on the preconsolidation 
q pressure, as the confining p pressure is reduced. 
Additional experimental workis required to establish the confining p pressure- 
- time- -wave velocity relations for silts and clays and for granular materials 
containing relatively small amounts of these fine-grained particles, It appears 
from the test datafor samples C-1 and NJ-1, given by Wilson and Miller, that 7 
relatively small amounts of the fine-grained soil particles exert a considerable 
influence on the dynamic behavior of a clayey or silty gravelly sand. . 
_ The writer again wishes to express his appreciation to the discussers for 
contributing much valuable information, It is hoped that efforts to obtain field _— 
by Smoots, Stickel, and Fischer, _ 
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at become a recognized part of the design procedure 1 for dynamically loaded 
foundations. It is also expected that testing programs similar to those de- . 
- scribed by Hardin, and by Wilson and Miller will soon provide information on 
the behavior of soils containing a significant amount of fine-grained particles. 
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 Errata.—p. 3, line 1 - “text” instead of “ “test”: 11, - - - ordinate c ot Fig. 

should instead of ”; p. 17, - 8 from 

be “has led” instead of “was led”; p. 27, - lines 1 and 2 should be “1,0 04” in- 

_ stead of “1.07”; line 2 should be «2. 45” instead of “2.52”; lines’ 3 and 6 should — 


instead of “1.19”; p. 27, - bottom line read 


(P+?) —mach (aa)? =: 3. 6 x 10° ft-Ib Ib sec”; 


lines 25, and 15 anne be “13.3” instead of 13. 4”; line 4 

“2,9” instead of “2.3”; line 6 should be “Fig. 9” instead of “Fig. 10” and should — 
be “0. 87” instead of “0.97”; lines 6 and 12 should be “3, 0” instead of “2.5”; : 
lines 8 and 10 should be “919” instead of “240”; and line 8 should be ‘3. 63” it in- 
stead of of 02”, p. 33, under reference 18, JW. Clawson 
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_SEEPAGE REQUIREMENTS OF FILTERS AND PERVIOUS BASES@ 


HARRY R. CEDERGREN, 24 M. J. Turnbull, F. pointed id 
x atendency, presumably more prevalent in past years than at present, for engi- » 
neers to call for drainage blankets, but not to rely on them to function. One 
—~ factor that may be largely responsible for this attitude is the fact — 
that filter design has been almost entirely on a qualitative basis, and designers | 
have had no way of knowing in advance of construction if filters of proposed — 


dimensions and materials would be capable of doing the drainage job expected 


and permeability of filter layers and pervious bases to the permeability of the 
soil from which water is being drained. 
a Possibly the manner in which perunechitity is defined has contributed to the 
_ Willingness of many engineers to allow relatively low permeability aggregates 
. _ to be used for filters. The coefficient of permeability actually relates to flow 
under a hydraulic gradient of unity. The flow potential in a given layer is pro- 


portional to the coefficient of permeability times the available gradient. = 
much for example, the available hydraulic gradient is 


0.05is onlyo one- -twentieth of that under ahydraulic gradient of ‘unity. This point 
is amplified by Fig. 11. Line 1-1 represents the slope of a pavement, which 
also is the maximum gradient that can exist in the drainage layer without a 
head build-up that will produce uplift of the pavement. Line 2-2 represents a_ 
a gradient of unity. For this gradient to exist under a pavement there 
7 _ Would need to be available ahead of 100 ft for every 100 ft of seepagedistance. _ 
A graded aggregate that would ¢ drain readily under this steep gradient might be 
- “rather slow-draining under conditions that actually can exist under pavements, 
The writer expects that much general thinking about the seepage capabilitiesof 
7 fone aggregates is based on a gradient of 100 ft in 100 ft rather than the 
small gradients that can exist in the field. 
* +S. Barber, M. ASCE, has shown by an approximate formula that the ca- 
7 pacity of horizontal drainage layers to transmit water to the sides varies in- 
_ versely with the square of the width of the layer. ‘He shows that a foundation 
_—- at Idlewild Airport with a width of 5,000 ft and an available head of 10 ft 
— ould be capable of transmitting seepage at a rate of less than 2 in, per yr. 
- The soil has a permeability of 100 ft per day. His computations demonstrate — 
a - the need for considering boundary conditions inthe design of drainage systems, 7 
r a Barber also — curves shin otha showing that the kind and quantity of 


—a October 1960, by He arry R, Cederg gren (Proc. Paper 2623). 
24 ‘Senior Materials and Research Engr., Calif. Div, of Highways, Sacramento, Calif, 


of them. The approach suggested inthe paper relates the physical dimensions 


on 
&g 
| 

| 

4 | 

tm 

2 4% 

ad 
4 
4 
— 


to one one- hundredth of the original value, w whereas 10%of 
4 fines reduced the permeability to one ten-thousandth of the original value. Only 7 
2% of silica fines reduced the seepage potential to less than $% and 2%of clay — 
lowered the water- -removing potential to about one-thirtieth! 


- _ Barber’ s curves lend considerable support to the hypothesis that Tealistic 


gre face of this kind of evidence (which substantiates general information 


‘that has been available for many years), itis amazing that clay-bound subbases 
- are stillconsidered by some engineers to be “pervious.” Simple computations 
7 such as those | made by Barber can demonstrate in afew ‘moments caine 


CMs tne in acy 
of such layers to serve effectively for subgrade drainage. 7 


-Hydraulic. Gradient = Unity 
par for computing 
"permeability coefficient 


um available 

“lo 


— 


"Pervious. 


Original 


e 
FIG, 11.—RELATIVELY SMALL ‘GRADIENTS AVAILABLE TO INDUCE 


 SEEPAGE IN PERVIOUS BASES FOR HIGHWAY 

— J. Gibbs, M. ASCE, referred to extensive filter studies by the U. S. Bu- 


4 reau of Reclamation, in which maximum flow is a principle basis for establish- 7 
ing limiting gradations of filter criteria that he are 
« a welcome addition to the paper. 


y 
broad limitations were allowed that an acceptable gradation could be miscon- . 


strued such that flow would be restricted and possibly be no better than the — 
q material protected. ” This statement emphasizes one of the main | points of the 


paper. . Under the low gradients that can exist in pavement drainage layers, ; 
_ the rate of escape of seepage can be unbelievably slow with aggregates that 
would | be entirel 


y adequate for filters fordams. A big difference can be attrib-— 
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narrow dimension of the. filter - (ag in most dam filters) can occur much more | 
readily than flow along the length of the filter (as in many highway filters). 
When this difference is taken into consideration, flow restriction is or ll 
tially ‘worsened in | pavement drainage. As “suggested by Gibbs, ‘determination — 
_ of the permeability of filter materials is an important part of filter tests . When | 
filters are required to function under low gradients, knowledge of the permea- 
bilities of available aggregates is essential to the development of adequate, ; 
economical designs, 
writer wishes to thank ‘Messrs. ‘Turnbull, Barber, and Gi for their 
a Filters and pervious bases are re serving an important function inthe seein 
mance of many engineering works. In the past, most filter | design (in terms 
of water- removing capabilities) has been entirely qualitative. While there have 7 
been a great many successes, there. have been some failures. 3 A rational ap- 
proach to filter design, with consideration being given to water-removing ca- 
pacity in relation to permeability and boundary conditions, can provide some 
assurance that designs will at least be “in | the ight county.” 


—Page > 16, 5, line 24, word at” be replaced v with “ to” page 19, 
7, word “orientation” has tranaposed letters; page 20, line 29, ‘word “had? 
should be. replaced with “has”; page 21, line 26, word “not” should be elim- 
inated, line 29, “14: 1 to 6:1” should read “4: 1 to 6: 321,” footnote 11 word “Me- - 


by “an un underlying page 30, line 9, word “gradients” should 
replaced with “gradings”; page 31, line 8, word “of” should be replaced with ¥. 
“or,” line 25, numeral “140” should be *180” and “0, 15 ft” should be “0.11 ft”; 


page 32, line 7, “0.7 ft per “0.07 ft per min,” 5th line above 
bottom, word “record” should be “second”; and page 33, line 1, “looking” should 
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Discussion by Stafford Happ 


Tuttle Creek 
Dam is an , interesting 1g example of careful embankment use of dredged fill and 
- ; broken shale and limestone. They have also reflected justifiable pride in the — 
- handling of foundation and rock fill problems, but their concern with the prob- 7 
lems may obscure the fact that success infitting design and construction meth- 
ods to the existing conditions also justifies the geologic forecasts which con-_ 


tributed to site selection This analysis is devoted to the geologic 


Favorable geologic features, as the authors have noted, included long Creek = 


limestone at a suitable elevation and location for the outlet works foundation, 
and Sallyards limestone at suitable elevation. These provide a favorable sur- 

face for placement of spillway concrete, with incidental removal from the foun- _ 

dation of the weak Legion shale, in a topographic saddle leading toa ‘relatively — 

> long tributary y valley which pr provides added safety against excessive spillway 
_ erosion, In addition to these obvious advantages, even the presence of Neva 

_ limestone within the spillway cut was considered favorable, relative to — 


7 nate sites in similar geologic terrane, because it insured exposure to cee 


vation and would facilitate any necessary remedial treatment. The authors 
suggest that calyx observation holes should be included at future projects where 
similar conditions are suspected. Calyx holes were judged unnecessary dur- 
“ing the Tuttle Creek investigations, and it seems doubtful whether they would 
have resulted in any appreciable change in plans, or saving in costs. 
_ Weakness of the lower part of the Neva | Limestone was recognized early ; in 
the investigation, and was the subject of a special study by S. S. ‘Philbrick, — 
Pittsburg District geologist who was called to consult onthe Tuttle Creek site __ 


because of his } experience WV with similar foundation materials. It was concluded, | 


_ - however, compensated by favorable features of the site. Later decision to — 


ehh that deficiencies could, if necessary, be made t up by quarrying Cot-_ 
_tonwood limestone along its outcrop upstream fromthe dam, at greater but not 
prohibitive cost. _ Thus a computed risk was accepted, offering some prospect — 
g significant economy, and protected by assurance of anacceptable alternative. 
Perhaps later disappointment with the Neva reflects some natural 
= toward wishful thinkin , one decision had been made to try to use it = as rock 
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a December, 1960, by K. S, Lane and R, G, Fehrman (Proc, Foor 2681), 
U. S, Atomic Energy | Comm,, Grand Office, Gr rand 
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economy wasa major factor inthe for berms- 
on- shale inthe spillway side slopes. There seemed no point in excavating more _ 
Shale than necessary, ¢ especially since laying back the shale slope would ex- y 
- pedite erosion of the overlying limestone. At that time, however, knowledge of 
"similar procedure was limited to a deep railway cut in eastern Ohio where the ~ 
_ berm was ; placed on shale vat base of a vertical face in sandstone, but with a 
_ thin limestone only a few feet beneath the top of shale, There, as Philbrick _ 
_ has reminded the writer, “The berm was preserved by the limestone and the 
_ primary criticism of the berm on Shale, which is its rapid weathering or —_ 
of durability, was overcome.” Hence, unless additional experience came to at- 
tention later, the Tuttle Creek berm design was essentially nt re ap tog 
reports” of future stability and maintenance experience should be helpful indi- 


cations a as to the desirability of similar practices elsewhere. Seepage is to be 7 


_ expected from the base of the Cottonwood limestone, providing a rigorous test _ 


of erosion resistance, 


The outlet conduit foundation involved a problem illustrative of those fre- | 
quently associated with minor geologic features, At the top was the thickest 
(approximately 1 ft) and firmest bed of the Long Creek limestone, essentially 

a the elevation desired for base of the conduit concrete. Preliminary drilling 
_ showeda weathered rind 1 in. or 2 in. thick at the top of this bed, with second-_ 
4 ary carbonate cementation, and also deeper weathering along vertical joints. | 


possible open solution channels, and foundation cleanup specifications stressed 
exposure of firm limestone in order to reveal any open channels along the joints, ; 

Removal of the weathered rind requireduse of jackhammer and rock spade, — 
vith which it was easier r to excavate the entire top bed than to chip off only 1° - 
in, or 2 in. of badly weathered limestone. Furthermore, weathering discolora- 
tion extended downward for several feet along many joints, although there were 
open channels. In these circumstances, there was some tendency simply to 
remove all discolored limestone, but this would have weakened the foundation 
by isolating pillars between the joints, and would have required excessive con- 
crete for replacement. | Thus the easy distinction by color had to be ignored, 
and reliance put on close and discriminating inspectionto assure sufficient ex- 
posure of joints to revealany open channels, yet avoidover-excavation of rock | 


which was firm although slightly weathered and discolored. | yay ey 
_ Reduction in riprap thickness on the flat 1 on 13 berm (as. described in the 
paper) seems reasonable in principle, but what are appropriate criteria for 
 eeeniekg the thickness needed? Presumably | 18 in. may not exceed the di- 
mensions of some of the individual p pieces. "7 There is a . scarcity of gubiliatend 
data on riprap design and experience, and it should be of interest if the authors 
could give a summary of the governing criteria. 
_ The data on density ‘of dredged sand fill seems clouded by the restriction of 
to “reasonably uniform non-stratified material,” and citation of only 
; the “more reliable” results. It would be helpful to know why stratified sam- _ 
ples were excluded, the range of values actually obtained, the proportion of © 
measurements excluded as “less reliable,” , and reasons for such exclusions, | 
Concerning protection of shale against ‘Slaking and 


. was handled. keep-— 


— 
— 
ae 
= 
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“surfaces moist without causing erosion or slumping, and 1 mo e details o on jet 
“how this was accomplished would be of interest. 


- A principle of possible wide significance is suggested by the statement that _ 


‘the | “old channel deposit” clay was | “similar to the terrace clay.” The latter 
is presumably a deposit of Pleistocene age, whereas the tine fill” is much 
nee, egg a depression or swale followed by intermittent tributary | 
_ drainage across the present-day river flood plain, and normally | boggy so that 

_ adequate sampling was difiicult. Presumablythere should have been much more — 
- opportunity, as there certainly was a more time, for consolidation of the 
clay. 

_ that age differences in tens of thousands of years s and soil loads of tens of feet 
are unimportant in consolidation of alluvial clays. — ‘This one case of course 

tohave more details onthe evidence > of similarity between the two clay: deposits. a 

Weathering was the most critical geologic factor in the Tuttle Creek foun- — 
dation conditions, marked by lightening of colors and limonitic stains, solution © 
of limestone, and ‘softening, cracking and crumbling ‘of both limestone and shale. 

- Weathering of limestones seems to be almost unpredictably erratic. Relative 
_ abundance of different clay minerals, or other impurities, isa suspected factor, 

but lacks verification as a general controlling factor. Possibly the relatively 7 
_ rapid and extensive weathering of the lower Neva limestone reflects greater 
abundance or a more susceptible type of clay content, but pre-construction — 
drilling results suggested that the more extreme Neva weathering might b be re- 
- lated to local absence, e, by erosion, of the higher Cottonwood limestone. The > 

- Cottonwood is a principal aquifer, and Neva weathering may have been retarded 

4 in areas where overlying Cottonwood served as a horizontal drain to intercept 

_downward- percolating waters. This hypothesis 1 might be confirmedor disproved | 
. additional information developed during construction, regarding the rela- 

_ tions between areal extent of severe Neva weathering and of overlying Cotton- 

Weathering of shale (or mudstone, claystone, or siltstone) seems to pro- 

gress relatively uniformly, although there were apparent differences among 

- Gifferent kinds of shale at Tuttle Creek. A large amount of laboratory test data 
_were accumulated, and it would be desirable to make this more widely avail- — 

- ghle by publication, especially if it could be correlated with observations on 
7 7 _ behavior in construction, Limited studies of some of the Tuttle Creek shales © 
"were ‘made by Ada Swineford, of the Kansas Geological Survey, indicating a 
predominance of illite which is considered among the more stable clay miner- _ 
als. ‘Hence extrapolation of dataon strength ar and stability of these shalesshould 
be subject to confirmation of similar mineralogy. Rate of disintegration and 
of shales in artificial embankments is obviously in | the 


factual data to support design based on reasonable assumptions. It would seem 
_ desirable to undertake investigation by sampling and testing embankments of 
_ various ages, known to have been built from various shales of various clay-_ 
_ mineral compositions. _ Embankments of canals and related dams should provide © 7 
a examples a century or more old, and there is an abundance of highway as well © 
: as dam embankments of lesser age. Philbrick has called attention to nearby 
Kansas highway fills constructed from shales similar to those at Tuttle Creek. 


| The whole field of weathering effects, relative to construction problems, should — 
be worthy of ‘Sponsored investigations. 
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y William Daniel Finn, BV. da Silveira 

WILLIAM DANIEL FINN, 13 M. ASCE. author is to be 
for his significant contribution to the theory of consolidation. His ingenuity in | 
applying the results of other fields to soil mechanics merits admiration. 
_ The writer has some slight reservations about the establishment of the ini- 

tial dial reading by the requirement that “consistent equal values of the co- — 
efficient of consolidation result from a comparison of compressions at differ- 


ent values of N.” Such a requirement is necessary | to satisfy the assumptions — 


of the theory. The initial dial reading, however, has an independent objective J 
existence; its value is indifferent to any set of thought processes brought to | 7 
bear onthe consolidation problem, The corrected initial reading coupled with | 
“the assumptions of the e consolidation theory may result inthe closest fit to avail- 7 
able consolidation data, but doubts may be reasonably entertained that it is the 
‘real initial dial reading. 
a As an extension of Scott’s work, an analysis of three-dimensional consoli- 7 { 
"dation with axialsymmetry is presented herein. Examples of thistype are (a) 
€ consolidation of a triaxial specimen in which radial and vertical consolidation 
- may take place and (b) ‘consolidation in the field where vertical sand drains ~ 
are used to accelerate the consolidation process. It is assumed that layers 
exist which permit vertical consolidation | 
‘The equation describing such consolidation is 
om” +t r 
r 
in which u is the excess over hydrostatic water pressure, t denotes time, Cc, 
‘represents the co- -efficient of consolidation for drainage in the Z or vertical ; 
direction, r denotes radial distance from the Z-axis and Cy is the 


for radial consolidation. 
< The may be into two. — 


92 


4 February, 1961, by Ronald F. Scott (Proc, Paper 2746), a > 
13 Asst. Prof, Civ. Engrg., Univ, of British Columbia, Vancouver, 
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DISCUSSION: SM SM 4 
Each part may be si solved separately with boundary conditions. So- 
lutions areto be found in the paper. _ The corresponding ng percentage consolida- 
tion for each case may t then be found. ee see 


The partial results” can be combined to the total U% consolidation at 


= (100 - Uy (100 ) (1/100) (23) 


which U%is the combined percentage consolidation, Uy%the percentage con- 


solidation due to vertical — percentage | consolidation due 
radial drainage. 


‘Then, as as in the p paper, 


However, now U(T) = n= 
solidation and Ty the time factor for» vertical consolidation, ‘The compression 
ratio is also = g (Cyy Cy) . Thus two values of the compression ratiotheoreti- 
_ cally suffices for the solution of Gi and Cy since they are assumed constant. —_ 


ee Since the equation wa was solved by parts, it is s preferable to determine C, _and 
~ independently as in the paper and ultimately compute U by Eq. a3. oat 


fetes problem. The consolidation may may, , for example, ‘be assumed to entirely 
_ due to vertical drainage. , If Eq. 24 is then used, an equivalent vertical consoli- _ 
dation co-efficient is obtained. In like manner, an ‘equivalent radial co- -efficient 


nay be determined. Either of these co-efficients may then be used to charac- 
terize the consolidation problem. The latter approach is analogous to trans-_ 
_ forming an anisotropically permeable section to one of uniform permeability ye 

-V. RANGANATHAMI6, _ The valuable ; paper by the author on a new 
od of consolidation-coefficient evaluation represents the only one of its kind 
ies endeavors to probe into what happens withinthe loading interval because — 

nothing but the average value canbe computed by the use of the hitherto known 
methods. This, though a significant contribution to our knowledge, seems to 

a sufferfromthe want of giving uniquely averaged value which is the result needed 


by the practicing engineer. 


This can best be demonstrated with reference to a laboratory time curve 
for atypical load increment (taking car care to choose one with little initial com- 


‘pression to eliminate correction) on a Saturated black cotton soil. The com- 


_ putation of Cy, the coefficient of consolidation, was conducted in accordance 


- with the example worked out in the paper, , but for a number of combinations of 7 


- initial time and N-values. Table 4 and Table 5 were p prepared from such com- 
putedvalues. The coefficient of consolidation is plottedagainst N-values aes 
3 family of curves (Fig. 2(a) ) which are described by the values of initial a 
_ written adjacent to them. From the same results, another family of curves of 


a “Simple Two and Three Dimensional Cases in The Theory of Consolidation of Soils,” 
by Carillo, Journal of Mathmatics and Physics, Vol, 201. 
. 15 “Hawaii’s Experience with Vertical Sand Drains, ad K. B, Hirashima, Bulletin 


in Soi ‘Soil Mechanics, Indian inst. of Science, Bangalore, India, 
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TABLE 4, —COMPRESSION DIAL READING* 


Compression dial reading (in.x 10™4) 


Time, in Black cotton clay Black cotton clay 

‘ 


1 
2 
4, 
6 
9 
12 
16 
25 
36 
49 
64 
81 
100 
121 
44 
69 
96 
5 


2 
0 
.0 
.0 
.0 
.0 
.0 
.0 
0 


TABLE 5.—COMPUTED® VALUES OF Cy 


Initial time, 


| 


Cy, insq in, per minx 10-4, 


- 
299 | 2482 
(b) BI ype | 
309 
4 


-coeffic 


(in*/mins!@ 


lidati 


~conso 


oO 


ation 


ient 


c 


| 


5 


3 4 


-0412.25 


Aw 


| | 


Black cotton s 
tonsa tt 


Ange 


— 
— 
— 
2,—COE FFICIENT OF CONSOLIDATION VERSUS N-VALUES 


> 


| 100 


-Cin/ mi 


- wo 


_N- “ratio | 


ding 


ion dial reading ten thousandths of an 
i 


a 


n di 


Compression. 


q 


: 


| 
pf times 
if 
| 


August, 1961 


pst 


PRESSUREO-s ton/sq.ft. 


= 4 
ly 


Tey 


Le} 


= 

| ii | | 
| 
4 
= 


= 


| 


he 


4 


| 


-= 


ECE 


~~ 


iz: 
| 
g 
— — -+——— — 
ca 
| 


Cyy versus N which are described byt the of final time written adjacent to 
- them is also produced (Fig. 2(b)). From these figures it can easily be dis- 

cerned that the coefficient of consolidation so computedvaries with the choice © 


of values of initial wane, final time and their ratio. i The extent of variation is 


reproduced in n Fig. 3. can be that there is no unique 


variation either. _ Although it is not possible to get the average coefficient af 


- only bring out the effect of the hitherto vaguely smemieates variable, second- 
arycompression, 
“Fig. and (b) show laboratory compression curves for two different 
soils, chosen: to demonstrate the other practical difficulty in having to properly | 
correct the zero reading. For the case reported in Fig. 4(b), uncorrected or 
improperly corrected zero reading gives absurdvalues for « compression rz ratios, 
These two aspects, taken together, will demonstrate the necessity for obtaining 
2 continuous record particularly in the early stages of compression. ae ae —_ 
From a general study of the family of curves produced in Fig. ‘it 
worthwhile examining rather 1 more closely inthe range e for time factor between 
0 and 0.1. Since U is an exponential function of T, it is considered befitting to 
use logarithmic plotting. | The Salient f features are self- explanatory f1 from the 
‘results of this further reproduced in Fig. 5. This among other 


N need not be limited only to values greater 


| ~ 


ASCE. —The method suggested. in the paper is not 


= identification of the ‘degree of consolidation corresponding to certain ‘dial. 7 
reading obtained duringthetest. The complete consolidation curve can be ob- 
tained from three dial readings in in substitution to the normal series of meas- > 
urements as standard procedure of consolidation tests. serve 
Such results are naturally possible of criticisms, because ‘it isa al 
many other advantages" must be recognized in the ideas of 
thor bythe analysisofhis work, 
‘The possibilities of the method for | the early stage of consolidation, where 
the parabolic form, 2v T/n, was “used, can be successfully ov overcome by 
using otherformulas for determining the ratio r = U(T)/U(NT). Itisa matter 
of more accurated numerical values and o a in plotting the e results — 
the graph (Fig. 1). 


that a required of consolidation U(NT corresponds to toa 
factor 0,2 or up and 50% of consolidation. 

; The correction of the zero reading can be obtained also eg! using Eq. 9 if the 
ratio (r) is well- -determined, thus dg > 


17 Con 7 Cons. . Engrg., , Com, . Solos e e Fundacdes, Rio de Janeiro, Brazil. 
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ON, 
This procedure is more general than those using the parabolic form, as- 
sumed as good approach for consolidation law at its early stages. — 
itt corresponds tothe case of radial drainage, mentionedin the paper, which 


Eq. 13 was this type of test = the writer18, 19° and 


Since § some discrepancies between test and numerical values ob- 
tained with this formula were emphasized, the writer re-examined the prob- _ 
= in Northwestern Tech. (1956-57) and checked the formula recently for the 
condition of equal strain with appreciable accuracy. | A consolidation curve was 
= by plotting Uy, against the ratio (T,/Ty - 50%). 
procedure has been successful for comparing different — 


urves in its was , and eile a basic seal time factor (T - 50%), for 


| 
| 
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es 4 e utilidade do Ensaio de Consolidacao de com m drenagem radial ex- 


Vol, 1, Session 1 to 13, Zurich,  .. 
7 20 Numerical Solution of Some Problems in the Consolidation of Clay, by R, E. Gib-| 
son and ‘P. ‘Lumb, Proceedings, British Inst. of Civil Engineers, ‘Paper 5877, London, 
a (21 Discussion by I. da Silveira of“ Procedure for Rapid Consolidation Test, ” by Hsuan- 
“Loh Su, Proceedings, ASCE, Vol. 85, No. SM3, February, (1959. — 
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GROUTING OF GRANULAR 


Discussion by Judson 


intelligent and exhaustive thinking processes prior to or to moving ahead ona foun-— 


dation treatment program by grouting methods, 


_ The statement that “. . . grouting has been handledalmost as an art, as Op-_ 


posed tothe science it should much more nearly be,” no doubt was deliberately - 
_ The principle contribution, that this paper makes toward advancing the “ art, Vd 
or science, of grouting isto emphasize the need of designers for extensive, - 
a telligent ‘and thorough investigation and exploration of proposed foundations, 
The next logical step is to assess and evaluate the data received, not only through 
the facilities of laboratory procedures, but throughthe studies and judgment of 
highly qualified geologists, and soils and foundation « engineers. Inthis special- 
ized type of work it has always seemed that the orderly sequence ofintelligent — 
engineering deductions were most effectively furthered and brought to a suc- : 
‘cessful conclusion by on-the-site testing under actual field conditions | in the 
material or rock itself to betreated. In many of the writer’s experiencesthis © 
was carriedone step further, or without a division line between design and con- 
"struction, into the actual construction work or contract itself and continued to 


make itself felt throughout the work. 
el Soils mechanics is considered to be a science by most of the. engineering 


profession, and yet some of our most famous authorities in this field have | 
adopted a personal and ‘somewhat “artistic,” in their own words, approach to 
extremely difficult problems and with great success, Geology has been con- © 
_ sidered a science for many years, and yet some of our foremost engineering — 
geologists have achieved notable success by. novel, individual, and ingenious 
theories, predictions ; and practices later proven sound as exploration and con- © 
The authors imply in the paragraph Ragga Mio goennt that the only mate- 
rials to be used for grouting clays are“. . . cement, or combinations of cement _ 
and other solids that will form layers or bulbs of hardened material.” 
It is extremely difficult to believe that, either in the laboratory or in the - 
field, clays in place can be grouted (that is penetrated and saturated) by any 
materials, solids or otherwise. It would be indeed interesting to have described a 
the field applications and test cores” or pits in which combinations of undis- 


turbed in-place sands, clays, and silts were consolidated _ — 


‘SgApril, 1961, by J. C. King and E, W, Bush (Proc, Paper 2791). 
Civ, Engr., Rock Island Dist., Corps of Engrs., Rock Island, ml. 
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— The statements are made t that. the way y to a choice of oun materials is 7 


J determined by e! effective pore size andassumption that all sand =“ in nature 
__ The implication seems to be from the statement, In other er words, if the di ; 
ameter of grouting particles is not over 1/5 x 1/3 = 1/15 of the particle di- 
7 ameter of the materials to be -grouted, the grout will pass freely through the 
void system,” that inasmuch as this formula was derived from some labora- 


tory source, presumably, ‘it should apply to similar materials found in — ] 


in nature, As tothe uses of fly -ashversus clays the writer feels, that 

volumes - when used with cements - being approximately the same, the size - 

and variations in size of openings and comparative costs should be the control- a : 

ling criteria assuming, naturally, ‘that the fly ash can meet ASTM standards. 

: It is alsothe writer’s opinionthat “bleeding” is of no particular importance © 7 
in pressure grouting underground as contrasted with the problem it often cre- 

ates in concrete or pre-packing operations with cements and gradedaggregates. 
_ The authors delve rather extensively into “Grouting Materials and their 
Limitations. ” It would be interesting to know the costs and availability in local 
7 markets of Portland blast furnace slag cement in this country. It is assumed 

~ “that they refer entirely to laboratory “ne orto carefully controlled nt 
_ Separated sand sizes where they state,” ', whereas for sand careful choice © 
- of the clay andcement enable one to omni the usefulness of cement asthe so- | 
lidifying agent at least one extra sieve size, as shown in Fig. 4 and Table 1 - 
. authors’ comments and conclusions on bentonites and chemicals will a 
be discussed except to say that they, inthe writer’s S opinion, have an extremely | 
limited use under certain, at times rare, foundation conditions. The assump- 
tion is taken that the sodium silicate grouts a are considered useful materials | 


ence inthe last 25 yr does not include what he would consider 
- of fine sands with sodium silicate or clays, together or separately. As to_ 
compressive ‘Strengths shown in n Fig. 6 and 7, were these derived from labora- 
tory cylinders orfrom 6-in. . cores recovered from field grouted areas? 
The question also arises as to whether the statements contained in the last 7 
‘joe under “Chemicals” has been proven under | in-service field operation 
or is largely a ‘matter of conjecture from laboratory experiments s and long- time 
7 ‘tests. The writer’s work has led him to the conclusion that the most oil 
= and valuable controlcriteria are derived in the field under existing conditions 
immediately prior to commencing the task of formulating and constructing, the — 
_ strengthening foundation or impervious barrier or combination of both, as may 
_ The writer finds himself forced to disagree that soft fine-grained satruated 
wa that cannot be practically grouted by true chemicals such as chrome- 
lignins or AM - 9, can be grouted with cement and cement-sand grouts. ie 
a Such a theory would seem to imply that ‘sands or solids can be grouted by 
other sands or solids with relatively the same void ratio, give or take one or a 
_— two sieve sizes. It is of course quite true that in some instances ‘sufficient 
:— ‘pressure: will permit injection through rupturing causing displacement. _ Whether 
“na “take” is low or high in soils under pressure does not necessarily seem an 
indication that — is not as much in the vertical } direction 3 as in wes 
wr 
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an 
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to obtain a large number of lly records of foundation 
-ance both where vibratory trouble occurred and where no suchtrouble was reg- 
istered,” because “ ... the progress of our r knowledge i in this field is depend- 

ent | upon the ‘systematic collection and analysis of records... . concerning the 


erformance of full scale structures.” 


; opportunity of examining certain aspects of the problem of resonance in mocking 
: - The paper lists five factors that influence the natural frequency of a soil-— 
foundation system. It is to be regretted that in the subsequent description on of | 
the troublesome machines these factors are not dealt with more fully. In par- _ 
ticular, the weight of the machines and their foundations as } well as the shape 
and magnitude of the contact areas would be of interest. 
‘It can be shown’ that, under certain simplifying assumptions, the resonant 
"frequency of a machine- -foundation- system can be byt the follow- 


in which fp, is the resonant frequency in cycles per unit of time; Cc denotes a 
coefficient dependent onthe Poisson ratio of the soiland the shape of the foun- 
dation contact area; VR is the propagation velocity of Rayleigh waves in the ~ 
soil; Y represents the unit weight of the soil; A is the contact area of the foun- aie 


dation; and W, is the weight of machine and foundation. rane 


A typical value of VR for dense sands would be 200 m per sec, Thus, Eq. © 
1 could be checked as to its S accuracy | in predicting resonant frequencies. 


April, 1961, by John P. . Gnaedinger (Proc, Paper 2793). 

Senior Lecturer, Israel Inst, of Tech., Haifa, Israel, 
_ 6 «Performance Records of Engine Foundations, ° by ( G. P. Tschebotarioff, ASTM 
Tech, Publication No, 156,1953, Vl. 11 
“Machine and Resonance,” by Geotechnique, , Vol. 11, 1961 
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| GROUTING TO PREVENT VIBRATION OF MACHINERY FOUNDATIONS 
Discussion by I Alpan and Judson Elston 
ALPAN,5 M. ASCE,—In order to appreciate the value of Gnaedinger’scon- _ 
a hypothesis concerning the operation of the structure under consideration, cul- -_ | 
: 
— 
— 
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itt is , stated inthe paper that « vit is probable that the etait constant will 


increase with increased amplitude of vibration 


ic the prior statement quoted implies a 
7 — experimental findings of H. Lorenz8, 9 who established soils to have - 


t? (sublinear) spring ch as show shown in ig. «6. 
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FOR 
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6, AMPLITUDE VERSUS TRANSMITTED FORCE 


his ¢ characteristic of the “soil spring” ' might explain the smaller than an- 7 
‘iia increase in natural frequency after grouting underneath the Expander _ 
foundation: instead of the computed ratio 558 ‘300 = 86, only” 402/300 = 340 


natural after grouting. This tendency could have been predicted 
Eq. 1 because yand, , consequently, Vp, would have ‘increased because of grout-— 


“Der Baugrund als Federung und D&ampfung schwingender Kérper,” by H, Lorenz, 
Der Bauingenieur, Vol, 25, J., Heft 10,1950, 
9 “Elasticity and Damping Effects Bodies ons Soil, by H. Lorenz, AST! 
Special Tech. Publication No. . 156, 1953, 
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The same curves ‘show the at the 1 new resonance frequen- 
F cies to be smaller than those at the original resonance frequencies. It is the 
writer’s s contention that tend to increase with their corre- 


From the phase diagram at resonance of a simplified dynamic model - 


in which X, isthe amplitude < at resonance; resonant | frequen- - 
CY; Mo is the eccentric rotating mass, r denotes the radius of neice c 4 
% the model the unbalanced force varies with the s square of the ede 
It appears from Eq. 3 that the amplitude may be expected to increase with 
the resonant frequency unless the damping is similarly i increased, rn aa 
od al In the case under consideration, the process of grouting had the effect of 
increasing the density of the foundation soil; existing evidence shows10 that ; 
densification increases the elastic modulus, the amplitude, and the wave 
ofa soil but reduces its damping characteristic. 


JUDSON P, ELSTON ASCE.—The writer is acquainted the large 
‘quantity of foundation work by grouting methods going on inthe Western Hemi- _ 
sphere, as well as in other parts of world. And there isa tremendous 
amount of work with baffling problems : as yet largely unsolved coming up in = 
_ immediate future inthetreatment of soil and soil-rock foundations by grouting 
or other foundation treatment it methods. . Unfortunately, in most cases novel and ; 
_ ingenious job or fieldsolutions are not brought tothe attention of the engineer- — 
ing profession, except through hearsay, but are lost in in prosaic job oo 
and semi-confidential project completion reports, 
It is noted in the Introduction, “Into that clay was pumped 34 cu yd of a 


7 


-ment and cement-sand grout.” Clark is discussing only thethick residual clays 
derived from disintegration and decomposition of the limestone itself. How- 7 
ever, the author gives no mechanical analyses of the thick residual clay. | ‘The 
writer assumes that the author is ; referring to uncompacted material that i 
sluffed from its natural position resulting in large cracks and openings. Per- 7 
Clark would kind enough to clarify this question. It is possible with 
sufficient pressures to rupture, heave, and displace soils or stratified sedi- q 
mentary rocks and nea m poecenn to fill up, by grouting, the zones of weakness, 
Clark does mention a“. . .gravelly pervious zone from 0 ft to 20 ftthick....7 


it is assumed that he may have had this specific material in mind when he re- q 
ferredtothickresidualclays, 

_ The pattern of grout hole andsequence of grouting ; shown i in Fig. 4 is simi- 
; dar to one experimented with by the writer in 1956 on the St. Lawrence re 7 

- ‘Project. A circular pattern was used, grouting from the centerout. The quan- 
tities of grout were controlled by varying mixes and pressures and stopping - 7 


a 10 “Ingenieur-Geologie,” by L. Bendel, Vol. I, p. 457, Springer Verlag, Wien, 1949, 
ll Civ. Engr., Rock Island Dist., Corps of Engrs., Rock Island, ne 
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all grouting in n the pattern after heavy mixes under pressure nent 
-turnedfrom all periphery holes, 
A somewhat related principle was utilized in a circular pattern test grout 
program in 1959 onthe Niagara Power Project to determinethe need for grout-— 
ings hole spacings; depths and grout acceptances. This particular pattern = 
the procedures used permitted the writer to determine grout travelin 4 uncon- 4 


nected and unrelated zone depths containing separately connected and discon- 


= cavities, seams and water-bearing openings. Information that was of 


inestimable value was found as to spacing of holes; packer settings, mixes, 
_ pressures and helped to evaluate the cost of sealing off several miles of foun- — 4 
priortoawardofcontracts, 4 


Th The writer takes exception to cement-sand mixes of; 1 bag cement: 2.22 | 
cu ft ‘sand; 0. 8 cu ft water (commonly known in the trade as a “brick mortar 
or dry pack mix” _ The writer’s faith in funnel viscometers is thus severely 
strainedat such times when required to convert into large scale (volume) field 
mixes some of the results of these laboratory instruments. ‘There is no 


of those ‘intimately acquainted with grouting, that the secret to success of the 
_ construction phase of this work lies withthe abilities and experiences of those a 


spect, as stated before, Clark not only had a a ‘difficult problem, but solved it “a 7 


ave satisfactory manne 


tyne o work is lacking af which comes hack to an oaninion held hv mos 
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"Discussion | by Jt Judson P, ‘Elston 


—Ker ’s paper is anoutstanding contri- 


JUDSON P. ELSTON,4 F. ASCE. — Kennedy’s 
“bution in the field of research with cement grouts and represents painstaking, 
thorough, meticulous effort to compile in one paper this voluminous subject. 
However, the writer f finds in the mortar” are cc that he is 


s fact, the writer feels the word “mortar” should be barred from the vocabu- 
ht commonly used by people in the drilling and grouting business. | 
yy _ This minor criticism does not reflect or imply other than whole hearted ad- 


“neering profession, and especially to the “art,” or science, of grouting. | ‘The 
writer believes it will become a classic description and compilation of the ad -— 
vances made in the field on cement grouts for founda- 


& 


-s April, 1961, by T, B, Kennedy (Proc, Paper 2794), 
4 Civ. Engr., Rock Island Dist., Corps of Engrs., Rock Island, fl. 


RESEARCH IN FOUNDATION GROUTING WITH CEMENT@ 
&g 
| and rescnect for the contrimition that Kennedy has made to ti 10°} 
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Ait, 


aft. 
Discussion by Judson P. Elston 


_ JUDSON P. ELSTON, F. ASCE.- —The writer | is S personally acquai nted with 


oa Polatty is aware of the writer’ 's vintinen on the use of sanded mixes inthe _ 
- field, and his reluctance and hesitation, at times, to consider seriously the use 
: of such mixes except when filling large open cavities and seams. The writer | | 
is of the opinion that all attempts to use sands in clay, silt or sand anes open- : 
: ings are doomed to failure, or only partial success at best. _ oe 
It is the writer’s wish that the Corps of Engineers could wie on a pro- 4 
_ gram of field testing under known conditions of clean cavities and seams versus ~ 
7 partially or completely y filled cavities and seams with sands, , silts, and clays § 


proceeding then to use the methods and so well described by the 


author, 


4 April, 1961, by James M, Polatty (Proc, Paper 279! 5). 
2 Civ, Engr., Rock tol: and Dist., Cor ps of mn, Rock Island, 
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pas year are by numbe below. ‘Technical- division 
aponnorshipis indicated by an abbreviation at the end of each Paper. Number, the symbols referring to: Air 
_ Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
« - draulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (SA), Soil 


- Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and Harbors _ 5 


(Ww), divisions. Papers sponsored by the Department of Conditions of Practice are identified by the sym- 
_ bols (PP). For titles and order coupons, refer to the appropriate issue of “Civil Engineering.” Beginning _ 
= with Volume 82 (January 1956) papers were published in Journals of the various Technical Divisions. To 
locate papers in the Journals, the symbols after the paper number are followed by a numeral designating ie 
the issue - wt particular Journal in which the paper appeared. For example, mane 2703 is identified as ve 


AUGUST: 2564(sM4), 2565(EMA), 2888(PO4), 2569(P04), 2570(HY8), 2571(EM4) 


2597(WW3), 2598(IR3), 2599(WW3), 2600(WW3), 2601(WW3), 2602(WW3), 2603(WW3), 2604(HW3), 
2605(SA5), 2606(WW3), 2607(SA5), 2608(8T9), 2609(SA5)°, 2610(IR3), 2611(WW3)°, 2612(ST9)¢, 2619(IR3)°, 
OCTOBER: 2615(EMS), 2616(E Ms), 2617(sT10), 2618(SM5), 2619(EM5), 2620(EM5), 
2622(EM5), 2623(SM5), 2624(EM5), 2625(SM5), 2626(SM5), 2627(EM5), 2628(EMS5), 2629(8T10), 
9630(sT10), 2631( POS) ¢, 2632(EM5)°, 2633(ST10), 2634(ST10), 2635(ST10)¢, 2636(SM5)°. 
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2697 

JANUARY: 2698(PP1), 2699(PP1), 2700(HY1), 2701(SA1), 2702(SU1), 2703(ST1), 2704(ST1), 2705(SU1), 2706 
(HYI), 2707(HY1), 2708(HY1), 2709(PO1), 2710(HY1), 2711(HY1), 2712(ST1), 2713(HY1), 2714(PO1), 2715 
2716(HY1), 2717(SA1), 2718(SA1), 2719(SU1)°, 2720(SA1)°, 2721(ST1), 2722(PP1)°, 2733(PO1)°, 2724 
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